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Abstract 


Until recently the drilling of a well was much more of an art 
than a science. With the advent of mathematical models of the drilling 
process it has become possible to estimate quantitatively the effect of 
certain of the parameters involved. 

The object of this thesis was to investigate various optimiz- 
ation procedures and how they could be used in conjunction with the mathe- 
matical model selected, in order to reduce the cost of drilling a well. 

The present knowledge of the drilling process restricted the 
number of parameters to be optimized to two, namely, the weight on the 
bit and the rotary speed. Other factors, such as mud properties and 
bit type, were assumed to have been properly selected. Within these 
limits, three methods of varying complexity have been developed. The 
first method seeks to minimize the cost per foot drilled during a bit 
run, the second method minimizes the cost of a selected interval and the 
third method minimizes the cost over a series of intervals. The methods 
are listed in order of increasing complexity. 

It was found that all of the methods gave a worthwhile cost 
saving and that the saving increased as the complexity of the method in- 
creased. The data requirements for the methods increased with increasing 
method complexity. 

Further improvements will depend on greater sophistication of 


the mathematical models used to describe the drilling process. 
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I. Introduction 


The nature aoe oi] industry is such that it is always facing 
uncertainty. It is never certain that any particular well drilled will 
be an economic producer, although the probability of finding a commercial 
well is, of course, much higher in a developed area than it is in a wild- 
cat situation. Because of the number of dry holes that do occur, it is 
necessary that the few successful finds give a large return on invest- 
ment. 

There are several ways that an oi] company can improve its 
chances of achieving a reasonable rate of return on investment. The 
first indication that oil or gas may exist in an area comes from the ap- 
plication of geological and geophysical techniques, so that any improve- 
ment that is made in the interpretation of this type is helpful. To 
this day, however, the only way of confirming the existence of oil or 
gas is to drill a well. Hence it follows that the lower the cost of 
drilling each well the greater the number of wells that can be drilled 
and the greater the probability of finding producing wells. 

In spite of this rather obvious fact, the drilling of wells 
was, until very recently, more of an art than a science. It was not 
until the late fifties and early sixties that the advent of mathematical 
descriptions of the drilling process changed the situation. 


The object of this dissertation is to examine how optimization 
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techniques may be used in conjunction with the published models in order 
that the drilling cost may be minimized. It should be realized that 

the total cost of a well involves other factors such as the amount of 
testing and coring and that the return on money invested in these areas 
Should also be considered. 

Any technique suggested will give poor answers if the model 
describing the system is inadequate. The second chapter of this work 
will therefore examine the progress made to date in forming a mathe- 
matical model of the processes involved in drilling a well. Chapter 
III will give the mathematical background to the various optimization 
techniques suggested and Chapter IV will review the results obtained 
using each method. 

In view of the fact that the models currently available may 
be improved upon in the future, every effort has been made to write the 
computer programs in a general format. This should permit, relatively 


easily, the substitution of any improved model in the programs. 
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II. Review of Available Models 


Although many prior investigations had taken place into the 


various aspects of drilling performance, Speer (28), in 1958, was the 


first to propose a generalized method for determining optimum drilling 


techniques. Speer suggested that the drilling performance depended 


primarily upon five controllable factors: 


is 
ex 
oe 
4. 
oe 


Properties of the circulating medium. 
The hydraulic effects. 

Type of bit. 

Weight on bit. 


Rate of rotation. 


Investigations since 1958 have confirmed and refined Speer's suggestions. 


to affect 
2 


The properties of the circulating fluid that have been found 
drilling are as follows: 

Mud weight: 

It has been shown (21,13) that drilling rate decreases as the 
mud pressure increases, and that the decrease is actually more 
correctly attributed to the excess of the hydrostatic pressure 
over the formation pressure (9,6). The reason for the decrease 
is thought to be due to the fact that compression of the rock 
makes it harder for the bit to break it up. In order to keep 


the penetration rate at a reasonable level the pressure, and 
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hence the mud weight, should be kept as low as possible while 
allowing for the highest formation pressures to be encountered. 
Solids content: 

Eckel (8) and Moore (20) have shown that the larger the solids 
content the slower the penetration rate. The reason for 

this appears to be because the particles interfere with the 
removal of the cuttings although the actual mechanism is un- 
clear. 

Mud viscosity: 

Increased drilling mud viscosity results in a slower rate of 
penetration (20). In the case of normal drilling muds the 
solids content and the viscosity are strongly interdependent 
and hence it might be thought that the viscosity was not an 
independent parameter. Moore (20), however, showed that the 
effect could be achieved using high viscosity glycerine. 

Moore suggested that the high viscosity resulted in lower 
fluid velocity in the vicinity of the cuttings at the bottom 
of the hole and hence the cuttings were not removed efficiently. 
Oi! content: 

Shales, especially, drill much more quickly when oi] is added 
to the drilling mud (20). This effect probably results from 
the oi] preventing "balling-up" of the bit. The economic ad- 
vantage of adding oi] is not always clear. The loss rate must 


be strictly controlled or else the cost of the oi] used may 
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well be higher than the cost saving it produces through in- 


creased penetration rate. 


The necessity of adequate hydraulic horsepower has long been 
recognized (1,14,20,25,28), but quantitative treatments are limited 
and generally hydraulic programs have been designed for a given rig to 
give maximum bit horsepower or impact force with the available pump in- 
put horsepower at maximum surface pressure. The thought is that, since 
actual requirements are unknown, the more horsepower the better (10). 

The function of the drilling fluid is to remove the cuttings from the 
area of the bit in order that fresh material may be drilled. Because 

of this it follows that the level of hydraulic horsepower required is 
related to penetration rate (28,3) and hole diameter (3) in that the 
volume of cuttings generated in unit time is proportional to the drilling 
rate and to the square of the hole diameter. 

Circulation rate must be considered as the upward fluid velocity 
has to exceed the cuttings' slip velocity in order that the cuttings be 
removed from the hole (14). 

Hopkin (14), discusses the factors affecting removal of cuttings 
during drilling,while the choice of mud properties and the hydraulic horse- 
power required,as dictated by differing hole conditions, is well documented 
by Reid (25). 

The bit selection is, naturally, very important and depends 
mainly on the degree of hardness of the formation being drilled. The 


Hughes Tool Company (15) has discussed, in general terms, the bit re- 
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quirements for drilling different formation types. It is not, however, 
proposed in this work to discuss the intricacies of bit selection,as to 
this day it is a matter for a certain amount of trial and error. 

The effects of weight and speed have been extensively investi- 
gated but it would appear from the variety of expressions that have been 
Suggested that the processes involved are complicated and require care- 
ful analysis. Bielstein and Cannon (1) were among the early investigators 
and they noted that drilling rate appeared to vary approximately linearly 
with weight and somewhat less than linearly with rotary speed. 

Similar results were obtained from the analysis of field re- 
sults conducted for the American Association of Oil Well Drilling Con- 
tractors (30,7). 

Moore (20) was the first to suggest an analytical method for 
finding the optimum weight and speed. He suggested that the drilling 


rate equation was: 


ala 
=i 
it 
= 
= 
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(II-1) 


where , F = footage drilled 
T = rotating time, hours 
k = proportionality constant 


W = weight on bit, 1000 Ibs 


N = rotary speed, RPM 


Bit life limited by tooth wear was considered but no formula 
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was given so that optimization under this limiting condition was not 
considered. An equation for bearing life was proposed and so optimization 
was possible for a bearing controlled situation. The equation Moore 


suggested is 
B= k /NW4 pane (11-2) 


where, B = bearing life, hours 
k = proportionality constant 
q = exponent on weight, Moore suggests 2. 
King (16) subsequently confirmed the form of equation (II-2). 
From previous data Moore found the values of the constants, 
and combining equations (II-1) and (II-2) with such factors as the rig 
cost, bit cost and trip time he was able to predict the optimum weight 
and speed at any point. The optimization technique suggested by Moore 
has been used extensively by other investigators on the more sophisticated 
equations developed later. 


Much effort has been concentrated on the drilling rate equation. 


All the early investigations gave results of the form: 
oF = k f(W,N) (11-3) 
qT * : 


Maurer (17,18) suggested that 


f(W.N) = N(W - W,)4 (11-4) 
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for "perfect cleaning", by which he meant instantaneous removal of al] 


the cuttings. The equation he proposed was: 


Z 
N(W - W_) 
ipo are one 
H Sy 
where, Wo = threshold weight before the bit penetrates the rock 
H = hole diameter, inches 
3} = formation drillability strength, constant. 


In this way Maurer related penetration rate to weight per inch of bit 
diameter and made the constant 34 independent of bit size. 

Perhaps the most complicated theoretical analysis done to 
date is by Outmans(24). The equation for drilling rate that Outmans 
derived described the rate in terms of weight, rotary speed and hydraulic 
horsepower at the bit. The resulting equation is complex and contains 
several unknown constants which would have to be estimated from previous 
experience in the same area. The simplifications that arise when weight 
or rotary speed is the only variable give a general insight into field 


results. If weight is the only variable, 


dF = ALC1ot Bu) 7? -7]° (11-6) 
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In this equation A and B are constants which include hydraulic horse- 


power and rotary speed terms. Similarly for variable rotary speed, 
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where A and B_ include hydraulic horsepower and weight terms. 


Examination of the asymptotes of equation (II-6) gives, 


1. °BW <<] then (1+ Bw)l/? = 14 5 BW (11-8) 
2. BN >> 1 then (1 + Bu)!/2 = (py)! /4 (11-9) 
Hence, for case l, 
creme een ; 
& = Ay BW) (11-10) 
and for case 2, 
dF | ; 
aT = ABW (II-11) 


Outmans showed that BW is much less than 1 for most rocks 
drilled by the oi] industry so that the linearity observed in field re- 
sults was a reasonable approximation. 

The equivalent equations in the case of variable rotary speed 


are: 


N (II-12) 
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when BN << 1, and 


jo (or 
1 


= AB (II-13) 


when B N >> 1. 


Thus the response to rotary speed will be at best, linear, 
and at worst, nil. The field results which show 


a no 


(11-14) 
where q <1 
are thus simply an approximation over a small range with q constant. 
This explains the reason for the variation in the exponent found by dif- 
ferent investigators. 

One other important point brought out by Outmans' work is the 
danger in examining one variable at a time as appears to have been done 


by all experimentalists. It should be noted that 


does not imply 
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where g(W) and h(N) are functions of weight and rotary speed respectively. 
This independence of the effects of weight and rotary speed cannot be 
theoretically justified and therefore the equations suggested may be 
Slightly suspect. 

The final step towards a comprehensive model was made by 
members of the Hughes Tool Company. Cunningham (5) first noted that, 
even when a homogeneous formation was being drilled, the proportionality 
constant of the drilling rate equation appeared to vary with the age of 
the bit. Galle and Woods (11) introduced the concept of bit dullness. 


The equations they developed are of the form: 


dre KrF ay 9) (lieie) 
d h, (0) 
where D = normalized tooth wear 


0 if new 


1 if completely worn out. 


dp Kof_(W)95(N) ; 
qi * ar (II-17) 
dB, 

ar = k3f3(W)g,(N) (11518) 


BY = normalized bearing wear 
0 if new 


1 if completely worn out. 
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The assumptions made by Galle and Woods will, of course apply 
to any optimization procedure which is based on their equations. These 
assumptions are as follows. 

i “Drilling cost 1s the summation of bit cost, rotating cost, 
connection cost and hoisting cost and these. costs can be 
expressed equivalently in terms of rig hours. 

2. Diamond bits are excluded from the analysis. 

3. Bit life is limited by bearing failure or tooth wear or 
drilling rate or a combustion of these factors. 

4. Circulating hydraulics are adequate and do not limit 
drilling rate. 

5. The drilling rate is a function of only bit weight, rotary 
speed and degree of tooth dullness; that is for the purposes 
of this investigation the effects of pressure, lithology, 
fluid properties, hydraulics and drill string dynamics are 
assumed to be unalterable at any particular depth. 

6. Within the range of rotary speed specified there are no re- 


strictions brought about by prime mover performance. 


Later investigations by Edwards (10) in 1964 and an on-well 
computer application by Young (31) in 1969 modified the Galle and 
Woods equations slightly but the industry (19) appears to consider 
that Galle and Woods equations are the best empirical set available at 


this time. 
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It cannot be said that the models proposed thus far are 
highly sophisticated, but the development is sufficiently far advanced 


to justify the use of optimization techniques. 
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III. General Theory 


In view of the fact that the Galle and Woods equations are 
the best currently available, the various procedures discussed will 
be formulated in terms of these equations. For the sake of clarity, 
the original presentation of these equations by Galle and Woods is set 


out in the following section. 


A. Summary of the Galle and Woods Equations. 


1. Drilling Rate Equation 


ale. 
ain ha eee (III-A-1) 


where F = distance drilled by bit, feet 


==) 
Il 


rotating time, hours 
Ce = formation drillability factor 


W = equivalent 7.88" bit weight = 7.88W/H 


W = bit weight, 1000 Ib. 
H = bit or hole diameter, inches 
k = weight exponent 


r(N) = exp(-100/N*) yo 428 + 0.2N(1 - exp(-100/N‘)) 
in the case of hard formations, and 


exp(-100/N2)NO-75° + o.5N(1 - exp(-100/N*)) 
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in the case of soft formations 
N = rotary speed, RPM 


a(D) = 0.928125D¢ 


+ 6,0) +. ] 
D = normalized tooth wear 
p = exponent on a(D) 
1.0 for flat-crested wear, 0.5 for A eae 


wear and 0.0 for "button" bits. 


The form of the functions in the above equation is such that 
the drilling rate increases as higher weights and rotary speeds are 
applied, and decreases as the bit dulls. The drilling rate is also pro- 
portional to the formation drillability parameter Ce which includes the 
effects of bit type, hydraulics, drilling fluid and the formation. It 
should be noted that, if it is decided that the bit type or the hydraulics 
or the mud system should be changed, the parameter will change. 

Formations of a very soft nature, for which the penetration 
rate is not a linear function of weight, are covered, to a certain extent, 
by the use of the exponent k. The value of 0.6 for k was found to be 


the best compromise in these cases (11). 


2. Rate of Dulling Equation 
dD _ i(N) ae 
ar (III-A-2) 


where Ae = formation abrasiveness factor 
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i(N) = N+ 4.438 x 1072 N° 


m(W) = 1359.1 - 714.19 logy, W 


The rate of dulling increases, as would be expected, as higher 
weights and rotary speeds are applied. The rate of dulling decreases 
as the dullness increases. This happens because the conical shape of 
the bit tooth results in a larger area being available as the tooth wears. 
As in the drilling rate equation, where Ce was affected by down-hole 
conditions, Ae in this equation includes the effect of bit type, hydraulics, 
drilling fluid and formation, and will be altered if any of these factors 


is adjusted. 


3. Bearing Life Equation 


wee 5 (TiteA=3) 


at 


where B = bearing life, hours 


rT} 


S = drilling fluid factor 


1(W) = 21340./(1 + 0.03W)?"*? 


The bearing life decreases as greater weight or speed is used. 
In this form the bearing life equation is not particularly useful, and 


it is found convenient to use the fraction of bearing life consumed. 


(III-A-4) 
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In this equation BY is the fraction of total life expended. 
Assuming constant weight and speed this equation can be differentiated 


to give the rate of wear: 


ai (III-A-5) 


Again, it should be noted that the parameter S is a function 
of the bit type, mud and hydraulics, and would change if the hole condi- 
tions or bit type were altered. In the case of the sealed-bearing type 
of bit theparameter for a particular bit type should, in theory, always 
be the same. In practice, however, some variation may be expected be- 
cause the bit action, and hence the forces on the bearing, is affected 


by down-hole conditions. 


B. Calculation of the Parameters in the Galle and Woods Equations. 


Before any optimization procedure can be used with the Galle 
and Woods equations it is necessary that the values of the three parameters 
C 


A, and S be known as a function of depth. The only way that the 


faa: 
parameters can be found is from previous wells drilled in the same area. 
From bit runs made in the previous wells the following information is 
required: 

1. Rotary speeds used during bit run - N, 


2. Weights used during bit run - W, 
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3. Footage drilled during each interval in which the speed 
and weight remained constant - AF; 

4. Time that each of the intervals occupied - AT, 

5. Final teeth condition - D 


f 
6. Final bearing condition - Be 
7. Hole size - H 

8. Formation response to rotary speed - soft or hard 
9. Formation response to weight - soft or hard 


10. Wear type - self-sharpening or flat-crested or "button" 


It is also necessary to rearrange the Galle and Woods equations 
in a form which permits solution for the parameters. 

Evaluation of the formation abrasiveness factor Ae is considered 
first . 


From equation (III-A-2) 
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é A(N) ee 
‘| aids | ae (I11-B-1) 


In the practical case the finite approximation to the integral 


on the right hand side of equation (III-B-1) is used. Hence solving for 
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si UP a m(W. ) 


bye Yo 3t02-- basqe yystor oF senagest aad *. a * . 
| iT bash yo S762 - inpiew od s2nogzey nobseerh: ig 
| : 
“Hordud" v0 ¢ batear9-f61 + 10 patnsqustiz-tfee - sqys e0W hist ‘a i 


J 


Lae) 


andtisups ebooW bas isd oft spnsv1s97 oF yiazegasn o2f6 at oe “at vse 
eystemer6q edt vot notfuloe estmagq Wo 


7 ae 7 
barebienos 2t 4A vas267 zesnovtzsrts notte? any to nofipuleva i 


; 
(SAT) aotdsups WOH * 


s 
| 3 
(f-8+J11) bea ; a = xb (k)s Ri 


> 


i alitietss eansH .beeu et (- Ug, aottaupe DS ‘ebhe basd 4 


” 


i 
| “(TAG WE a 
($-8-i11) iT RG tou 24h 


19 


where the weight and speed are constant over the interval aT, and 


n 
Nir ir (TI1=B=3) 


Te - total rotating time 


and 


ae ‘ 
| a(x) dx = U(De) (III-B-4) 
0 


De - dullness at the end of the bit run. 


It should be noted that equation (III-B-2) has little meaning 
when "button" bits are used and there is negligible tooth wear. When 
this is the case the value of As tends to infinity and the rate of dulling 
to zero. If this is the case, it is better to set the rate of dulling 
to zero and ignore any slight dulling that occurs. 

The drilling fluid parameter S is calculated from the integrated 
form of equation (III-A-5): 


iP 
| qa (II1-B-5) 


where Bf is the bearing wear at the end of the bit run. Rearranging 


equation (III-B-5) to solve for S: 
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The parameter C¢ is calculated differently for the several 
wear types. In the case of "button" bits the procedure is as follows. 


From equation (III-A-1) with D = 0, 


dF ok 
dt = Ce WriN) | (III-B-7) 
AF . 

or igo —k 


C,=2 ; Rh (I11-B-9) 


j=l Alj Wr (Nj 


Self-sharpening and flat-crested wear may be considered to- 


gether. Combining equations (III-A-1) and (III-A-2), 


aa, = 
C, Wr(N) Az a(D) m(W 
of Te afss (2 Syad piode ton | AN at gana (III-B-10) 
(a(D))Pi(N) 
—k i-p 
C, Wer(N) m(W) (a(D A 
fo te ent a a feel Sal (III-B-11) 


By exactly the same procedure as shown when finding Ae 
equation (III-B-11) may be integrated and solved for Ce. 
The answers for self-sharpening and flat-crested wear are 


slightly different because of the different value given to p in each 


DDI 88 
eee aa 


dais “1 


hay 


(wh 9° FS 


74 
(im) fF 49 a : 
Lae | 
j eyo VOT pntvioe «mot ennsvettth af hart at 
34 fi i ; 
ps —- & 7) on. 7 
(4-111) (aH ETA Sy nn ¥ vie ne 


iio 


wt beisbtemos sc yen veew beszev2-ssel? bis potneqinne-F fae ” wee! : 


bee ‘) 
(S-A+LET) ne (T-A-ITT) 2nobdeups,porntdea® ary, 
i.e (eugty 
a (in (a)s.4 (hlaAH 9 Be . ) if 
a (A) F9((O}a) if 
7 La 34 9(a)s) @a (MW 59 a 
a 34 patbatt nedw nwofle 26 syubs201q smsz ait s capali hi 
i .- 142 10? dav [oe ‘bap! badsypesat od 


te | P . 16 Aeow betzsvo-t6 17 nn tieine tart pag 


it. figsa nt q of navip sufsv siiaaileliaaa tan 
hie | 


zi 


case. 
If the wear is self-sharpening, p = 0.5 
; 1 : i(N.)AF. 
: sept III-B-12 
FA WE) 521 WE rN) m(Wy) 
where De 1/2 
V(D_) - | Ca(x))'"" dx (III-B-13) 
0 


If the wear is flat-crested, p = 1.0 


i Ne beAF 3 


sled Sab 
2 ON, t pol (IITI-B-14) 


f j=l War(N) m(W;) 


Thus all the equations necessary for calculating the parameters 
have been derived. For convenience the necessary equations for each 
wear type are summarized below. 

1. Self-sharpening wear - Ag from equation (III-B-2), S from 

equation (III-B-6) and Ce from equation (III-B-12). 

2. Flat-crested wear - Ae from equation (III-B-2), S from 


equation (III-B-6) and C, from equation (III-B-14). 


f 
3. "Button" bits - S from equation (III-B-6) and Ce from 


equation (III-B-9). 


The derivation of the equations from which the parameters can 
be found includes one important assumption, namely that the parameters 


are unchanged throughout the bit run. It is impossible to discover the 
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variation of the parameters with depth unless the bit is pulled at various 
times during a run and the wear measured. Obviously this would be much 

* too costly and the only method that can be used is to examine the drilling 
rate profiles of the bit runs for any sudden changes in rate. Combining 
these rate changes with knowledge of the lithology, the values of the 
parameters as a function of depth can be estimated in order that a better 
approximation to the actual conditions can be found. Certain of the 

bit runs will, in fact, be made through one formation type only, and 

these are used for control on the values of the parameters. In this way 
the accuracy of the parameters will be improved as more and more wells 


are drilled in the same area. 


C. The Galle and Woods Optimization Procedure. 


There has already been developed what might be called "point" 
optimization. The solution of the governing equations has been described 
(11,12,2), but in each case the solution obtained is not truly optimal, 
as several important points have been ignored. However, as this solu- 
tion is used as a stepping off point for the other methods to be dis- 
cussed, it will be examined in some detail. 

The object is to find the combination of weight and rotary 
speed which results in the lowest cost, the weight and speed remaining 
constant throughout the bit run. 


There are three possibilities, any of which may control the 


answer. 
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Teeth control bit life. 
Bearings control bit life. 


Drilling rate controls bit life. 


Each of these conditions must be investigated when deciding 


on the optimum weight and speed. 


out is as 
ita 
2% 


The information required before the optimization may be carried 
follows: 

The parameters found by the method described in part B. 

The rig cost - Cae dollars/hour. 

The bit cost - Cc. dollars. 

Trip rate - t. - time to go one foot (both directions) during 
trip. 

The change time - t - time required to change bit at surface 


plus the time required to change the number of drill collars. 


Before the Galle and Woods equations can be used in any optimi- 


zation procedure it is necessary that they be converted into a form 


amenable to solution. In the case of self-sharpening or flat-crested 


wear this 


dullness. 


is done by expressing each equation in terms of changes in 


Hence by combining equation (III-A-1) and the inverse of 


equation (III-A-2) we have: 
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Inverting equation (III-A-2) gives: 


A, a(D) m(W) 
Ss Maryse (III-C-2) 


Combination of equation (III-A-5) with the inverse of equation 


(III-A-2) gives: 


DS (III-C-3) 


It is convenient to simplify the above equations by defining 


the following three functions: 


f, (WN) = AT (III-C-4) 

f,(W.N) = #5 (WN) C Wr (N) rreces) 

Theale I11-C-6) 
f,(W.N) = wane (III-C- 


Substitution of these functions in equations (III-C-1) through 
(III-C-3) gives: 


ap =f) (WN) a(D) (LilsCe 7) 
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= fp(WyN) (a(D))-P (III-C-8) 
dB _ 
ap = f3(W.N) a(o) (III-C-9) 


For any section drilled at constant weight and speed these 


three equations may be integrated as follows: 


Do 
Alert, (W,N) [oa (ITI-C-10) 
D, 
AF_ = f,(W,N) oe, dx (iii Ge 1a) 
S eee 
us 
Do 
AB. = f(W,N) | a(x) dx (III-C-12) 
D 


In these equations AT. is the rotating time, AF. is the footage 
and AB. is the change in normalized bearing wear which occurs when the 
dullness changes from D, to D, and the weight W and speed N are applied 
to the bit. 

As previously stated this derivation applies to self-sharpening 
or flat-crested wear only. For "button" bits the derivation is simpler 
because the rate of dulling equation is omitted. 


The inversion of equation (III-A-5) yields: 
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Combination of equation (III-A-1) (with D = 0) and equation 


(III-C-13) gives: 


dF 
No a) eee ee (III-C-14) 
BY N 


By defining two functions equations (III-C-13) and (III-C-14) 


are much simplified. 


Hence: 


(III-C-15) 


(III-C-16) 


Substitution of these functions into equations (III-C-13) and 


(III-C-14) and integration between the two values of normalized bearing 


wear (during which the weight 


AF. = 


and speed remain constant) gives: 


(III-C-17) 


(III-C-18) 


The object of the derivation is to obtain a formula for the 


cost of the bit run. 


In view of the fact that the rig cost is usually 


described in dollars per hour it is convenient to express the cost in 


terms of rig hours. Thus the 


trip time, rotating time and change time 
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can be used directly as cost items. The bit cost is divided by the rig 
cost per hour to get its equivalent cost in rig hours. Remembering that 
the trip rate is defined as the time per return foot, the simplest way 
of calculating trip time is to multiply the trip rate by the average 
depth of the bit run. With these points in mind, the total cost per 


foot can be expressed as follows: 
Gree (t (ee Fear /2) qet cee Gr teertAtajiite (III-C-19) 


C, = cost, rig hours per foot 

t. = trip rate, hours per return foot 
e. = depth at start of bit run 

Ce bit cost. dollars 

Cc, = rig cost, dollars per hour 


t_ = change time 


Some of the terms in the above equation are functions of W 


and N, while the rest are constant. Separating these two types: 
Ga (c. + t AF. (WN) /2 + AT. (WN) )/OF (WAN) (II1-C-20) 
where CAS ees Cy /c. ale (III-C-21) 


Galle and Woods, by utilizing the calculus of variations, 


solved equation (III-C-20) for the optimum weight and speed. The de- 
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tails are fully discussed in their paper (12). Because of the complexity 
of the derivatives it was necessary to solve the resulting equations 
numerically. In view of this, a direct search technique could equally 
well have been used. 
The step-by-step procedure by which the optimum values of W 
and N are found will now be outlined. 
1. Choose the search area of W and N. If, after the search is 
completed, the best weight and/or speed coincide with the 
lower limits of the search area, then these limits may be 
too high, and the search should be repeated with different 
limits. 
26) Caliculatcer (WM N)stt 


(WN) and f,(W.N) if the wear is self- 


] 2 3 
sharpening or flat-crested; otherwise, calculate f, (WN) and 


f,(W,N) and then proceed to step 5. 
3. Assume final dullness is 1. 


Calculate By from equation (III-C-12). 


tf 


li Bf alk De is indeed 1 - proceed to step 4. 
ih Bie at rs 1 and must be found. From equation (III-C-12): 
Ds 
Be = ]= F(W,N) | a(x) dx (III-C-22) 
0 
ue 
f(W,N) | a(x) dx = = 00 (I1I-C-23) 
0 


This equation can be solved for Dy using the Newton-Raphson 
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technique (4): 


D (II1-C-24) 


Tee Pe a(D 
Equation (III-C-24) is repeated until the convergence 
criterion is satisfied. 

Calculate AT, and AF, from equations (III-C-10) and (III-C-11) 
respectively; go to step 6. 

Calculate NTE and AF. from equations (III-C-17) and (III-C-18) 
respectively. 

Calculate the cost per foot from equation (III-C-20). 

Using a search technique, repeat steps 2 through 6 until the 
best combination of W and N has been found; these values are 
called me and N. 

If a "button" bit is being used, proceed to step 10. Other- 
wise, reduce De by 5% and using the values YW and N found 

in step 7 to start with, recompute the lowest value of cost 
per foot for this new value of De 

If the cost found in step 8 is greater than that from step 7, 
go to step 10. Otherwise, proceed with search on De to find 
the best value and the associated cost and me and N”. 

The procedure is complete and the best operating parameters 
and the resulting time and footage have been found. 


It should be noted that it is not necessary to examine the 
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rate controlled situation for a "button" bit because, with no tooth wear, 
the rate does not alter with time. 

This method, although simple, suffers from several grave dis- 
advantages. The procedure does not allow for variation in A 


Cz and S 


fie oT 
during a bit run. No allowance is made for formation changes. A hypo- 
thetical case will show how this can lead to errors. Suppose a soft 
formation of interval 600 feet is being drilled, followed by a very 
hard formation. The "point" optimum procedure is being applied and gives 
400 feet as the optimum footage for the first bit and 410 feet for the 
second. The second bit is thus pulled approximately half way through 
its suggested run. The question then arises as to whether one bit could 
have drilled the interval more cheaply. Alternatively, two bits, each 
drilling half the distance, might be best. The reason for this uncertainty 
is that the method is truly applicable only where the section being con- 
sidered is very large compared with the distance drilled by each bit. 

One modification that can be made rather easily to the basic 
Galle and Woods method is to consider the intervals in groups. This re- 
stricts errors to the end of the last interval in the group. The pro- 
cedure is essentially the same as previously outlined except that after 
each optimum has been found, the associated footage is checked to see 
whether it extends further than the end of the interval. If it does, 
the dullness and bearing conditions at the end of the interval are calcu- 


lated and these values used as the starting points for optimization of 


the rest of the bit run. 
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In this way, reaching the end of an interval need not imply 


changing a bit, unless of course, it is the last interval of the group. 


D. Modified Galle and Woods Method. 


This modification of the previous method is designed to over- 
come all the objections of the previous method with the exception of vari- 


ation in Ces Ae and S. The well is divided into intervals where A,, C 


feet 
and S are constant and each interval is optimized separately. This method 


can therefore be called "interval" optimization. 
The procedure is outlined as follows: 
1. Select the interval and the constants. 


2. Using Wain and N calculate f, (WSN); f,(W.N) and f,(W,N) or 


min a 3 

f, (WN) and f, (WN) depending on bit type. If a "button" bit 
is being used go to step 4. 

3. Assume De =e: 
Calculate Bf from equation (III-C-22). 


If Boe < 1, De = | as assumed - proceed to step 4. 


ak 
Otherwise calculate D, from equation (III-C-24). 

4. Calculate AF, from equation (III-C-11) or (III-C-18) depending 
on bit type. 


5. Calculate the minimum number of bits required: 


M, = Integer part(Interval/AF, ) + ] (III-D-1) 
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For this number of bits calculate the distance drilled by 


each. 
AF. = Interval length/Number of bits (III-D-2) 


Calculate f,(W,N), etc., as before. If a button bit is being 
used go to step 11. 
From equation (III-C-11): 

D 

= iP ie 
f,(W,N) | (a(x)) Bieditee AF. = 0 (III-D-3) 
0 

and again using the Newton-Raphson root finding technique: 


D 


ie 
-p r S 
| (a(x)) dx #5 (WN) 
Z 0 nts 
D. ahaa i> a ee OA (III-D-4) 
jt j a(D¢ ) 


Equation (iII-D-4) is repeated until the convergence criterion 
is satisfied. 


Tf De >", tne values of W and N are outside the range of 


possibilities; go to step 13. 
Calculate Bie from equation (III-C-22). 


If Bye > 1, the values of W and N are outside the range of 


f 
possibilities, therefore go to step 13. 


10. Calculate AT. from equation (III-C-10); go to step 12. 
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From equation (III-C-18): 
ep ecael Se (III-D-5) 
xe ‘1a 


If Be PoiGO tes Lepw a. 

Calculate ae from equation (III-C-17). 

Calculate the total cost of the interval and hold the values 
of the cost and associated weight and speed if the cost is 
the lowest to date. 

The total cost of drilling the interval is n times the cost 
of one bit drilling at the average depth. This is because 
the symmetry means that the bit is above the mid point of the 
interval for exactly half the time, and below the mid point 
exactly half the time. Thus each trip can be considered to 
have been made from the mid point of the interval. The cost, 


ani rig Bours, 1s stnus: 
C, =n (ae + c/c., pitts (e. + e,)/2) (III-D-6) 


Continue the search for the best values of W and sea a 2 

other grid points remain), by repeating steps 7 through 12. 
If not on the first cycle, compare the best cost with that 
found for the previous number of bits. If the latest cost 


is not smaller, go to step 16. 
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15. Increase the number of bits by one and return to step 6. 
16. The optimum number of bits is one less than the current 


number. 


This procedure thus finds the lowest cost of drilling an 
interval on the assumption that it is best that each bit drills the 
Same distance. 

An examination of this assumption will now be made. 


Suppose an interval y feet in length is to be drilled, then: 
n 
jooee 2 f(z.) (III-D-7) 


and 


n 
) hay y (III-D-8) 


where Z. is the distance drilled by the ae 


bit and 1 is the total 
rotating time for the distance y. Differentiating to find the approximate 


change in total time for small changes in the distances drilled: 


: of(z,) 


ce 
a 
tt 


Az. (III-D-9) 
; 


(Note that double use of the subscript i implies Summation) . 
Now since the sum of the distances is constant, the algebraic 


sum of the changes in the distances must be zero. That is: 


n 
Tez 10 (III-D-10) 
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If the derivative is an increasing function of Z.: 


of (z + 6) ; afl 2 aay) of ( 
z + 6 ees a way - 
mean mean mean 


and hence MS will be positive for all displacements of Z from eens 
In fact, the derivative will be an increasing function of the distance 
drilled in almost all cases. This can easily be seen by considering the 
following example. If a bit can drill 100 feet when weight W and Ny 
are applied, it will be able to drill 110 feet only if the weight and/or 
the speed are reduced. This reduction means that the time to drill the 
first 100 will be longer than before and to drill the full 110 feet will 
be considerably greater than 110% of the first time. Hence the variation 
with footage is greater than linear so that the derivative is an increasing 
function of the distance drilled. The best policy, therefore, as far as 
the rotating time is concerned, is to drill the same distance with each 
OMe 

The total trip time for n bits in an interval can be found by 
considering the following. The first trip, at the top of the interval, 
is one way only, so that the trip time is t e,/2. At the end of the 
first bit run the depth is e. + 2, and as the trip is made both ways the 
trip time is Pee + Z,)- Similarly at the end of the second run the 
depth is e, + z, + Z and the trip time 1st (e. + °Z) + Z5) The trip 


time at the end of the rah 


bit run is tee eda) Fics Zeer ine 
final trip is one way and takes t,, e,/2- The total time spent tripping 
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Zot tie ee oe ee /2) (III-D-11) 


It is convenient to replace e/2 by eas e./2 and to complete 
the series by adding and subtracting en (which is equivalent to e. z Zz, ue 
. + Z,) When this has been done the total trip time is written as 


follows: 


- e./2 + (e. + Z,) +... + (a. + 


CU T= a)ra9) (i i ==12) 


It is now possible to find the total cost of drilling an 


interval y drilled with n bits. 


ca oat 


ye ips 7 Cu cpm ee (I11=D-13) 


In general the form of Vy is not known. If, for example, it 


can be considered as a quadratic function of the distance drilled, the 
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COSE 1S: 


(a, + a.z. +a - 


kal ae aD tnCec(n alee 


(in tal = i)z,) +n c./C,, +n te (III-D-14) 


Applying the method of Lagrangian multipliers to equations 
(III-D-14) and (III-D-8), as shown in Appendix A, the lowest cost is 


found to occur when: 
ee) ecnieeel ) (II1-D-15) 


Thus, if the rotating time can be represented by a quadratic 
function, the lowest cost occurs when the distances drilled by successive 
bits increases by t /2a,- 

That it will always be true the later bits should drill further 
can be seen from the following argument. It has been shown that the 
rotating time is minimized when the distances drilled are equal. Hence 
the only way of lowering the cost by altering the distances drilled from 
the mean is by lowering the tripping time more than the drilling time is 
increased. The only way of lowering the tripping time is by reducing 
the distance to be tripped. This can be achieved by making the distances 
drilled by the bits greater as one proceeds down the hole. 


Hence in general: 
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cy aa cage oe (III-D-16) 
2 3 


In the general case, the problem of allocating a distance to 


be drilled to each bit can be solved by means of dynamic programming. 


minimize 


where 


and 


Using the terminology of Nemhauser (22), it is required to 


n 
C., = by ra(x. sd.) (III-D-17) 
n 
) F =y (III-D-18) 
i= 
Xa ts (x, .d,) ToS a a (III-D-19) 


In these equations rs (x, .d5) is the cost of commencing in 


state x and making decision d., and ts (x, sds) is the transformation 


function which results in state Xe] when decision d, is applied to 


the system in state x 


i: 
In this case the transformation is simply: 


Xe 9g = Xe = d. (III1-D-20) 


Defining f5(x;) as the minimum cost of an i-stage process 
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starting in state Xa» there exist the following recursive relations: 


F(x) = HulD r(x, »d)) (III-D-21) 
=X 
ie 
f.(x.) = ee (ri (x;sdj) + fay (x, - ds) (III-D-22) 


In practice there will be an upper and lower limit on the 


value of d. which will result in limits on state Xa At the jth stage: 


ined 
max 
j Shia <x, nin} i | (TII-D-23) 
There will also be a minimum number of stages: 
We = Integer part(y/d_..) + ] (III-D-24) 


Equation (III-D-16) imposes further limits on d.. It becomes 


(allowing for the fact that the last stage is designated number 1): 
ip kdoe- (deuce od, (III-D-25) 
Since a fel ici Cate Reeaey oe aacile le (III-D-26) 


it. fol lows? that: 
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Hence only values of d. which satisfy the above inequality need 


be examined. 
The optimum number of stages is that which produces the lowest 
value of f.(y). 


Iney.cost function is: 


fal Oe alg Vas + Ch/C,, tet te, Say d./2) (III-D-28) 
This is the minimum cost associated with commencing a bit run 
x. feet from the bottom of the interval, and pulling the bit after drilling 
d. feet. The rotating time "a, is assumed to be already optimized, in 
other words, the best combination of W and N to drill d. feet has pre- 
viously been found. 
The procedure for solution is as follows: 
1. Calculate f, (x) for X, over the range specified in (III-D-23); 
set i = 2. 


* 


Zee tcaculate fF. (x) = y(xe,d.) For, 1(% 


* * A 
- d.) where d. is the 
Teas i7-1°°% j 


value of the decision i which gives the lowest cost and d, 
is found by a search on d.. 
3. lf 1 1S Greater Chan Mos go to step 4; otherwise, increase i 
by 1 and return to step 2. 
eit f.(y) < f._y(y), increase i by 1 and return to step 2. 
Tt F.(y) > F,_yfy). decrease i by 1 and go to step 5. 
5. The optimum value for the number of bits is the current value 


of i; trace back through the stages to get the optimal policy. 
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There are thus three possibilities which have to be tested in 
connection with "interval" optimization. These are: 
1. Each bit drills the same distance. 
2. Each bit drills a certain distance more than the previous one. 
In other words the increase from bit to bit is constant. 
3. The distance drilled by each bit is atadestey by a dynamic 


programming technique. 


E. "Multi-interval" Optimization. 


The methods described thus far have foundered when the para- 
meters in the three controlling equations have changed. In some situ- 
ations this is not important because a change in the parameters signifies 
a change in the rock type which requires a change in bit type. If this 
is the case, the bit has to be pulled and replaced by another type. 

In many cases, however, the same type of bit is to be used 
in the new formation and it is then open to question as to whether, just 
because theparameters have changed, the best point to change bits is in- 
deed where the next formation starts. 

The method about to be discussed involves the problem of 
deciding where the optimum point is for changing the bit. The analogy 
between this problem and, for example, catalyst replacement is obvious. 
Dynamic programming can be used to find the best policy. The discussion 


that follows, closely parallels that of Roberts (26) in his catalyst replace- 


ment problem. 
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Some definitions are necessary: 

B. = the state of the normalized bearing wear at stage i. 

F = the state of the normalized tooth dullness at stage i. 
f<(D, ,B,) = the minimum cost of an i stage process starting 
in state i and following an optimal policy. The quantity i 
refers to the number of stages remaining. In this problem 

a stage is characterized by a certain footage drilled. 
g,(D;.B, .W. oN.) = the cost of stage i for a system originally 
in state D.,B, and where W. and N, are the weight and speed 
chosen for stage i. 

he (Di, 7 Big Wey Na 47) = the transformation function which 
gives state D; .B, at stage i as a result of starting in state 
Dy4,> By4, at stage i+] and applying Way and N.,, during 
Stage itl. 


R= replacement cost at stage i. 


Because each stage is a certain footage drilled, it is obvious 


that, after replacement, the stage footage remains to be drilled. As 


f (0,0) is the cost of drilling the remaining i stages optimally, com- 


mencing with a new bit in the hole ready to start drilling, f, (0,0) + Ry 


is the cost of replacing a worn bit with a new one at the beginning of 


stage i and drilling the remaining i stages optimally. 


The recursive relations are: 


fF. (O,0} 0s aC gs (0,0,Wa Na) (iit-ksi) 
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£10) 58, ) sain sien: | (III-E-2) 
WON, R, + f,(0,0) 
f.(0,0) = _min ew oa (III-E-3) 
alba A (ARIEL 


ain jer a eee eel 
= + 
f.(D.,8.).= min.| "4°"4l €...(h, .(D,.B..M. Ne] . (11T-E-4) 
risen Wen, 1b Meet ks aa 
art R, + #,(0,0) 


The equations are evaluated in order and then equations (III-E-3) 
and (III-E-4) are repeated for successive values of i up to n. 

The recursive equations have now been established in general 
form, but, of course, it is necessary to construct the equations in terms 
of the Galle and Woods equations for actual solution of the problem. 

The transformation equations which apply in the case of self- 
sharpening or flat-crested wear will now be derived. Recalling that the 
final stage is given the number 1, what is required are equations re- 


lating D to D; and B;_, to B,. 
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Rearranging equation (III-C-11) gives: 
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defining D 


] 
Q, = | (a(x))P dx + AF./f,(W,N) (III-E-6) 
) 
and solving for D. using the Newton-Raphson root-finding technique, 
the recursive relationship becomes: 


D 


2. 
D aide be J LSP cei 1-p dey 
mega we (a(x)) " dx - Q /(a(D2 )) (III-E-7) 


Rearranging equation (III-C-12) in order to solve for B.: 

= 2 

WN) | a(x) dx + B (III-E-8) 
D 

Equations (III-E-7) and (III-E-8) can be generalized by in- 


serting D;_, for Dj, 0; for D,, B;_, for B, and B; for B,. The equations 


then become: 


(D5). 
J ie 
fava? ax = Q, 
L Dicer) = (Dym)), = a (III-E-9) 
ID Siw es speaed itev@ tre Shon 
J 
z net 
By_1 = Fa(W,oN,) f a(x) dx + B, (111-E-10) 
D. 
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In the case of a "button" bit there is only one state variable 


and hence only one transformation equation. Using the same procedure by 
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which equation (III-E-10) was calculated, equation (III-C-18) may be 


transformed to give: 


B,_; = AF,/f,(W,.N,) +B, 


; ; (III-E-11) 


This completes the derivation of the transformation equations. 
It is now necessary to express the cost of drilling stage i in terms of 
the rotating time. The equations are different for self-sharpening or 
flat-crested wear, and for "button" bits. 

In the case of self-sharpening or flat-crested wear the cost 


is obtained by generalizing equation (III-C-10) to give: 


to give: 


- B.) (III-E-13) 


The replacement cost is 
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R. = C/c. i es te. (III-E-14) 


is the depth at the start of stage i. 


This completes the derivation of the required equations, and 
possible to outline the method of solution. 
Set up a grid of normalized dullness and bearing wear for 


each stage of size L. by M. such that: 


=) 
— 
—i 
~~ 
iT] 


B. (1) = 0 (III-E-15) 


p.(L;) = B.(M.) = 1 (III-E-16) 


5 * he * 
Calculate g,(0,0,W, .N,) where W, and Ny are the best values 


of W and N as found by a search technique. Thus F, (0,0) has 
been found. 
Calculate the total cost of changing the bit before commencing 


the last stage and then drilling the stage: 


Gee et (080) (III-E-17) 


1 


Evaluate f(D, »B,) for all other values of D, and B, between 


1 and L,-1, and 1 and M)-1, respectively. Obviously, for edge 


1 
points, the bit MUST be changed, and hence 


f,(D,.1) = #,(1,B)) = C, (III-E-18) 


pw ke oe) ’ 
7. i 


(af-3-T11) 0 = (f),8 = (1) (0 
(af-3-111) he Ginga Gaye — Sees ve een 


. * 
esulsv deed ald 975 rit brs ;W sod C(t fle. 0) -@ adsfuoted oe 
. i 7 


26 (0,0); 2udT .oupinasys rovsse & yd bnuot 26 Yons Wo 
pnfonsmmoa voted tid oft patpnsio to S202 Istot edt stefuotsd .€ 


repste2 ort patifinb neat bow spste teof ent 


(Xf-3-111) (0,0), + + (A= 3 


47 


9. Using the recursive relationships of equations (III-E-3) and 
(III-E-4), calculate f;(D; 5B.) for each grid point of each 
Stage until the final point is reached. 


6. Work back through the data to get the optimal policy. 


Several points should be noted: 

1. Only values of W and N which allow transfer from D.,B, to 
De Be without making Ds and/or Bs exceed | are ad- 
missible at any stage. 

2. Minimum values are assigned to W and N in order that the 
search area is cut down. If the minimum values violate 
condition (1) the bit will have to be changed. 

3. In the case of "button" bits there is no dulling of the 
teeth, and hence the state variable D. does not apply. The 
state variable grid is thus variable in B. OnlVe i.e. 1b IS 
a vector in B.. 

4. Close examination of the controlling equations for self- 
sharpening or flat-crested wear shows that the state vari- 
able B. does not affect the optimum procedure except where 
the bearings wear out first. If the teeth control bit life, 
therefore, the state variable grid can be reduced to a vector 


of D. values only. 


This completes the theory behind the methods investigated. 
The results obtained from each method will be discussed in the next 


section. 
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IV. Discussion of Results 
A. Introduction 


The difficulty in obtaining good data has restricted the 
amount of testing possible. Only recently has it become the practice 
of the major oil companies in Canada to keep sufficiently good records 
to make possible application of any of the techniques developed. Where 
records have been kept, sufficient data is still classified as confi- 
dential thus making it impossible to compare the methods with a well 
that was drilled without the benefit of optimization methods. Mobil Oi] 
Canada Limited were prepared, however, to release a set of data which 
consisted of their interpretation of the equation parameters as a function 
of depth for a complete well, plus the necessary cost information, and 
how the formations might be grouped together. This information proved 
to be quite sufficient for testing the comparative worth of the different 
methods suggested. 

Before the Galle and Woods equations were used, it was necessary 


to calculate approximating functions for: 


The function 1(W) was given in tabulated form by Galle and 


Woods (12). It was found that a hyperbolic fit was best. 
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The polynomial approximation of degree four represented the 
integral very well. It was discovered later that an exact solution! 
existed but the polynomial was retained because of its greater simpli- 
city. The final results, together with the appropriate computer pro- 
grams, can be found in Appendices B and C respectively. 

Apart from the papers of Galle and Woods (11,12), and Billing- 
ton and Blenkarn (2), no original data was available for testing the 
program that calculates the parameters. However, the program was 
thoroughly checked by calculating the parameters from the results of 
the optimization programs and comparing with the parameters actually 
used in these programs. The listing of the program and some typical 
results are given in Appendix D. 

The Mobil data used in the optimization is listed in Appendix E. 
The programs and results for "point" optimization, "interval" optimization 
and "multi-interval" optimization are given in Appendices F, G and H re- 
spectively. Each of the methods will be discussed separately and then 


a comparison of the methods will be made. 


B. Point Optimization 


The simple version does not have a separate program written 
for it but can be taken as a special case of the "multi-interval" modi- 
fication with each of the intervals considered separately. The improve- 


ment using the "multi-interval" modification was significant. Certain 
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sections were not changed by the modification but, comparing the sections 
where a change did occur, the times involved were 1161 rig hours and 1132 
rig hours respectively. This represented an improvement of 2.5% as a 
result of the modification whilst there was no change in the amount of 
computing time. Each computer run on the IBM 360/67 took approximately 
O.1 minutes. 

Because this program took so little time to run, it was con- 
venient to use it for testing the different search techniques suggested 
by Wilde (29) and Rosenbrock and Storey (27). It was found that the 
more complex methods suggested such as "Partan" and Rosenbrock and Storey's 
rotating co-ordinates method did not improve on the simple gradient 
method (with weighting factors). The weighting factors which gave the 
best overall results were 1 for the weight and 4 for the rotary speed. 

This method is quite adequate for the situation where neither 
the intervals nor the constants are known to any degree of accuracy. 

If used with a little common sense, under these conditions it should give 


reasonable results. 


C. "Interval" Optimization 


This was the second method tried. The most interesting result 
was that the modifications of the method proved to be wholly ineffective. 
The results were 1303.44 rig hours for the basic method, 1303.10 rig 
hours for the direct search method and 1303.51 for the dynamic programming 


procedure. Close examination of the answers suggested that the only case 
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where any appreciable savings might be made is where the drilling rate 
is extremely fast and the trip rate is extremely slow. Needless to say, 
this is a very unlikely combination, and so for practical purposes it 
can be concluded that making the distances drilled by each bit the same 
is quite satisfactory. 

The dynamic programming technique gave the worst result by a 
small margin. In general, the program correctly assigned greater foot- 
age to the later bits but occasionally this was not true. For example 
in the interval 5070' to 6770' the second bit drills further than the 
rest. The reason for this would appear to be that the large number of 
interpolations result in smal] errors which are enough to influence the 
final answer. 

The computing time for the simple procedure was 0.16 minutes, 
while that for the more complex methods was about 2.3 minutes in each 
case. Clearly, this alone would require a reasonable improvement when 


using the more complex methods. 


D. "Multi-interval" Optimization 


The most complex method is the "multi-interval" technique and, 
as might be expected, the largest number of problems were encountered in 
perfecting this method. Originally it was assumed that the 2-D grid of 
the state variables was necessary at all times, but as explained in 
Section III-E, it is only necessary to have a 2-D grid when the bearing 


wear is the controlling factor and there is appreciable wear on the teeth. 
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As the time required for this method increases exponentially with the 
number of dimensions, the time saved by this observation was considerable. 
In the case of the Mobil well, the time to compute the optimum procedure 
was reduced to less than 5 minutes from 30 minutes, while the accuracy 
was improved by using a larger number of elements in the dullness vector. 

Obviously the starting point for a gradient search should not 
influence the resulting minimum. It was found, however, that at certain 
times it did make a large difference and so it became obvious that the 
function being searched was not unimodal. The reason for this became 
clear when the tables of optimum cost were examined. 

As shown in Figure 1 the previous cost function increases up 
to a certain point and then becomes constant. The change occurs when it 
becomes cheaper to replace the bit before commencing the stage. Obviously 
once the bit is replaced the cost is the same no matter what the condition 
of the old bit was. Adding the cost of the current stage to the previous 
cost gives the total cost to date. The cost of the current stage de- 
creases as the weight and speed increase but, simultaneously, the amount 
of wear during the stage increases, resulting in a higher previous cost. 
The form of the two functions is such that the total cost function has 
two minima, one internal, and one at the maximum allowed values of weight 
and rotary speed. Which of these two minima is the lower depends on the 
two curves which make up the composite curve. 

Thus it is necessary to find both minima, and this can be done 


by starting at the lower limits of W and N and finding the internal minimum, 
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FIG. 1 COST VARIATION WITH WEIGHT AND SPEED 
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and starting at the upper limits of W and N and finding the minimum there 


eee and Lane may be in the forbidden region in which case it is necessary 


to find the edge of the allowed area). The lower of the two minima is then 
taken as the true optimum. 

Because of the time factor involved in increasing the number 
of state variable grid points, extensive investigation was carried out 
into the effect of altering the number of grid points. In conjunction 
with this, variation in the number of intervals (or more correctly stages) 
and the degree of interpolation to be used in computing the previous cost 
were investigated. The results are presented in Figures 2 through 5 and 
Tables 1 through 6. 

Figures 2 and 3 and Tables ! through 3 show the effect of varying 
the number of intervals. The apparently haphazard plot is easily explained. 
The formation which was used for testing was of such a length that it was 
best drilled by 3 bits. Thus when the number of intervals is a multiple 
of 3, the cost is lower because each bit drills the same distance. If 
the number of intervals is not a multiple of 3, the bits cannot drill equal 
distances and the time taken is longer. This illustrates one point that 
must always be kept in mind, namely, that the number of stages must be 
chosen well. 

Figures 4 and 5 and associated Tables 4 through 6 illustrate 
the effect of increasing the number of grid points. The general trend 
is for the results to improve as more grid points are used. It would 


also appear that the matrix size should be larger as the number of intervals 
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FIG. 2 COST VARIATION USING FIRST-ORDER INTERPOLATION 
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FIG. 3 COST VARIATION USING SECOND-ORDER INTERPOLATION 
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Table 1 


Cost Variation with Number of Intervals 


Using a 5 x 5 Matrix 


No. of Cost Using Lee Cost isthe ghd 
Intervals Order Interpolation Order Interpolation 

6 55.64 55.64 

7 Sy fe al 57.44 

8 56.58 56.69 

9 55.46 55.40 

10 56.79 56.64 

11 56.64 56.38 

12 5553 55.44 

ie 55.80 56.16 

14 56.04 55.64 
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Table 2 


Cost Variation with Number of Intervals 


Using a 9 x 9 Matrix 


rc a eS 


No. of Cost Using a Cost Using and 
Intervals Order Interpolation Order Interpolation 
6 55.60 55 997 
7 57.08 57.06 
8 56.03 56 726 
9 bbaa/ 55 35 
10 56.729 56.11 
11 56.25 56.20 
2 55500 55.41 
13 56.31 55.80 
14 56.59 55e93 
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Table 3 


Cost Variation with Number of Intervals 


Using a 13 x 13 Matrix 


No. of Cost Using ee Cost. Using pnd 
Intervals Order Interpolation Order Interpolation 

6 SosD0 55. 50 

7 56.96 56.99 

8 56.30 56.24 

9 Seiscisl 55735 

10 ale 24 56.05 

4 55.91 55.69 

12 55.46 55.28 

13 Soul 7 S508 

14 55.61 55.67 
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Table 4 


Cost Variation with Matrix Size Using 9 Intervals 


an oe eee 
e344 55.58 5557 
oy ets) 55.46 55.40 
6 x 6 55.63 55.47 
1 oO 55.55 Be 43 
8x 8 55.46 55 35 
9x9 55.37 55.35 
10 x 10 55.48 55.36 
Wx 55.43 55.35 
12. xel2 55.35 55.35 
13 xels 55.38 55.35 
14 x 14 55.38 55.36 
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Table 5 


Cost Variation with Matrix Size Using 13 Intervals 


A ee MS Bee 
ie spiel 57.03 
ae 55.80 56.16 
ce Boece 56.06 
7x7 55.90 55.98 
Bx 8 55.89 56.41 
9x9 56.31 55.80 

10 x 10 55.79 56 93 
Wx 55.67 55.88 
12 x 12 55.70 55.65 
13 x 13 55.77 55.63 
14 x 14 55.64 55.72 
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Table 6 


Cost Variation with Matrix Size Using 18 Intervals 


a eee 


Matrix Size cose ian or Gost sing 2 on 
4x4 57.21 37.17 
5 x 5 57.32 96.29 
6 x 6 55.97 96.62 
PSU 55.47 95.27 
8x 8 55.33 95.38 
9x9 55.36 95. 30 

10 x 10 55.26 99.22 
Wx 55.49 95.24 
Tene 55.44 95.27 
130X013 55.45 95.25 
14 x 14 55.28 95.26 
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per bit increases. For example, in Figure 5 an 8 x 8 matrix is required 
before stability occurs with 3 intervals per bit (9 overall), and a 
10 x 10 matrix if 6 intervals per bit are used (18 overall). Also the 
results when using an odd number of intervals (i.e. in this case a non- 
multiple of 3) are more scattered. 

The conclusions may be summarized as follows: 

1. Quadratic interpolation may be considered slightly preferable 
to linear interpolation. The improvement is, however, small. 

2. The size of the state variable matrix needed depends on the 
number of stages per bit. The suggested rule, which should 
have an adequate safety margin, is that the number of stages 
should be no greater than the number of row elements of the 
matrix divided by 2. Hence with a 9 x 9 matrix, the number 
of stages per bit should not exceed 4. 

3. The number of stages chosen is important. If, for example, 
the optimum number of bits to drill an interval is three, 
then, the number of stages should be a multiple of three. 
This accounts for the strange curves obtained when plotting 
the cost against the number of stages. The effect of in- 
correctly choosing the number of stages produces a large 
error, but this error does decrease as the total number of 


stages increases. 


The results of the optimization of the Mobil Oi] data were quite 


encouraging, and it is believed that no difficulties should be encountered 
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in using the method. Although, as previously stated, the number of 
stages was important, it was found to be quite easy to select the number 


correctly, two or three trials being sufficient. 


E. Comparison of the Methods 


It can be seen from the results summarized in Table 7 that 
as the complexity of the method used increases, the amount of savings 
increases. The rig cost was $200 per hour so that the saving between 
using the best and the worst method is approximately $20,000. The 
charge for use of the IBM 360/67 at the University of Alberta is $6 
per minute which is low by commercial standards. However, even assuming 
a figure of $10 per minute, the cost of a run using the most time-consuming 
method is $50, while that for the simplest is about $1. It can thus be 
seen that the saving that should be realized by use of the most. compli- 
cated method is well justified. 

The results show that the methods which have a practical use are 
numbers 2, 3 and 6 in Table 7. These are the modified "point" method, 
the "interval" method (in its simplest form) and the "multi-interval" 
method. 

It is suggested that the modified "point" optimization procedure 
be used when there is little data available. It is very fast, and it 
would be possible to run the program several times, at little cost, using 
different intervals and parameters over the likely range in order to get 


an approximate optimum procedure. As more knowledge is gained about the 
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Table 7 


Comparison of Different Optimization Techniques 


a ge ee ee ee eee eee 


Method Optimization Cost Saving over Percentage Computing 
No. Type (Rig hours) Method 1 Saving Time 
eee renee Mee eee Se a= ed ee ke ae eg tle ee ee 
] Pont | 1366 - 0.10 min 
2 Modified 1337 29 25 0.10 min 
“point. 

3 "Interval" 1303 63 4.6 0.16 min 

4 “Equal incre- 1303 63 4.6 2.5) min 
ment interval" 

5 "Interval by 1304 62 4.5 2.30 min 
dynamic pro- 
gramming" 

6 "Mul ti- 1265 10] Vhs 4.88 min 
interval" 


geology, and more bit runs are made, the "interval" method becomes more 
appropriate, and finally, probably after three or four wells have been 
drilled, the "multi-interval" optimization should be used. 

It is important to realize that, even if the predicted per- 
formance does not match the field performance very closely, this does 
not mean that the procedure is useless. It will always be true that 
the procedure can act as a guide towards better operating procedures. 


One other use of the programs can be made. This is to examine 
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the effects of increasing the maximum weight and rotary speed. The 
programs can be run with very high maximum limits in order to see what 

the cost saving might be. This type of analysis is useful, for example, 
when hole deviation is a problem and the decision on the use of stabilizers 
is to be made. A good economic decision can be made by determining 
whether the cost of renting the stabilizersis recovered because of the 


faster drilling rates. 
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V. Conclusions and Recommendations 


RE 


The object of this investigation was to examine the progress 
made to date in developing a satisfactory model of the drilling process 
and how optimization procedures could be used in conjunction with the 
model in order to reduce drilling costs. 

The mathematical models available at this time have the dis- 
advantage of including parameters which vary with operating conditions, 
so that changes in bit type or fluid properties alter the parameters. 
Thus any experimentation interferes with the optimization procedures. 

However, assuming constant down-hole conditions are maintained 
from well to well, the following conclusions can be drawn: 

1. Three useful methods have been developed for optimizing the 
drilling of a well. These are: 
(i) "Point" optimization modified in order to consider 
groups of intervals. 
(ii) "Interval" optimization with the bits drilling equal 
parts of the interval. 

(iii) "Multi-interval" optimization. 

2. The methods listed above are in order of increasing complexity 
and along with that complexity goes the need for better data. 

Hence the method used should be chosen according to the amount 


of data available. 
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3. The stage distribution for the "multi-interval" optimization 
has to be carefully selected. In order to do this, several 
trials will probably be necessary. 

4. The cost of running the programs is small compared with the 
predicted saving. | 

5. The programs can be used to investigate the effects of raising 
the weight and speed limits. If a large cost saving results, 
this may spur investigation into how the limits might be 
raised. 

6. Any new development should be incorporated in the optimization 


methods that have been derived. 


It is suggested that the methods developed are as sophisticated 
as needed with the models developed thus far. It is therefore recom- 
mended that future work be concentrated on developing more adequate 
mathematical models of the drilling process. 

The author believes that the methods developed, if used with 
care, can act as a guide toward good drilling procedures and that the 


cost saving should be considerable. 
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Nomenclature 


formation abrasiveness factor. 

bearing life, hours. 

vector of normalized bearing wear values at stage 

i in a process or an element thereof. 

fraction of bearing life expended. 

final fraction of bearing life expended. 

change in normalized bearing wear over a section. 
formation drillability factor. 

total cost, expressed in rig hours, of drilling an 
interval with n bits. 

total cost of replacing bit and drilling one stage . 
total cost, expressed in rig hours per foot, of a bit run. 
normalized tooth wear. 

final normalized tooth wear. 

vector of normalized tooth wear values at stage i of a 
process, or an element thereof. 

distance drilled by bit, feet. 

total distance drilled by bit, feet. 

length of stage i of a process. 


h 


value of the 3° footage during a bit run. 


length of a section. 
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bit or hole diameter, inches. 

number of elements in vector D.. 

minimum number of bits to drill an interval. 

number of elements in vector B.. 

rotary speed, r.p.m. 

vector of possible rotary speeds at stage i of a process, 
Or an element thereof. 

matrix of best rotary speeds at stage i of a process, 
or an element thereof. 

value of the cu rotary speed during a bit run. 

Cost of replacing bit at stage’ i. 

ariilang: fiuidst accor. 

rotating time, hours. 

total rotating time, hours. 

rotating time for stage i. 


th 


rotating time during the j” interval of a bit run. 


rotating time to drill a section: 


0.309375 D> + 3D“ + D. 


4 


0.113328 D* - 0.412768 D° é 


+1229230 De +e O19 7eeDr 
bit weight, 1000 Ib. 

equivalent 7.88" bit weight = 7.88 W/H. 

vector of possible weights at stage i of a process, or 
an element thereof. 

matrix of best weights at stage i of a process, or an 


element thereof. 
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r(N) 


value of the Aa. weight during a bit run. 


0.9281250° +6 D+ 1. 

bit cost, dollars. 

fixed cost of a bit run, expressed in rig hours. 

rig cost, dollars per hour. 

vector of possible footages drilled at stage i, or an 
element thereof. 

vector of best footages drilled at stage i, or an element 
thereof. 

depth. 

depth at beginning of stage i of a process. 

depth of the bottom of an interval. 

depth at the top of an interval. 


N + 4.438 x 107° N°. 


N 
weight exponent in drilling rate equation. 

1 in most cases but 0.6 for very soft formations. 
21340./(1 + 0.03W)9-4 

1359.1 - 714.19 10g, 9W . 


number of stages in a process. 


exponent of a(D) in drilling rate equation. 


1.0 for flat-crested wear, 0.5 for self-sharpening wear and 


Os0 form button bits. 
exp(-100/N2) N°°428 + 0.2 N(1 - exp(-100/N*) ) 

een Oy hell) 
in the case of hard formations, and exp(-100/N") N + 


0.5 N(1 - exp(-100/N°)) in the case of soft formations. 
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time to change bit and alter the number of drill collars, 
hours. 

trip rate, time to go one foot (both ways) during trip. 
vector of footages remaining to be drilled at stage i of 
a process, or an element thereof. 

length of an interval, feet. 


distance drilled by ae bit of a series, feet. 
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Derivation of Footage Distribution 
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Appendix A 
Derivation of Footage Distribution 


From equation (III-D-14) 


n 
P 2 
C. hy (apta,Z,+a,z5) + t (e. (n+1) : (Ge + e)/2) 
n 
ce A (n+1-i)z, Pee C Ca (eee (1) 


For a given value of n, certain of the above terms do not vary. 


Hence defining: 


C= t Ce. (n+1) - (e.te,)/2) ton Ch/c. on te (2) 


xs 
krr15 


n 
Gee a) (aytayz,+a,24) + t. 


(nt1-i)z, + C (3) 
=] i 


It is desired to minimize Cy with respect to Z;, subject 


to the condition: 
)z. = y = constant (4) 


Using the method of Lagrangian Multipliers (23) a new cost 


function is constructed: 
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where L is the Lagrangian multiplier. Differentiating equation (5) 


with respect to z5: 
a, + 2aoz, + t(n+1-7) +L=0 (6) 


Adding the n equations obtained in this manner gives: 


n n 
na, + 2a, egret Y (nt+1-i) + nb = 0 (7) 


; (n+1-i1) = n(nt1) - n(nt1)/2 
ey tual) (8) 


n 
and) 2, = y by definition. 
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Hence equation(7)may be modified to give: 


a, + 2a,y/n + t.(nt1)/2 ile=20 (9) 
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Hence solving for L gives: 


L=- (a, ~ 2ay/n + t(n+1)/2) (10) 


Substituting in equation (6) and simplifying: 


Zaz, - 2aoy/n + t ( (nti) /2-4) = 0 CT) 
Hence solving for z. gives: 
= See(oq 
z= y/n + ta, (21 - (nt1)) (12) 
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Appendix B 
Program and Results of 1(W) Fit 
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C MAIN 
MAIN 


THES PROGRAM IS AN ADAPTATION OF AN 1eBsMe APPLICATION 
PROGRAM CALLED 1130 DECLINE CURVE ANALYSIS PROGRAM 
(1130=MP-03X)e THE PROGRAM CAN BE USED TO APPROXIMATE 
DATA BY EITHER AN EXPONENTIAL OR HYPERBOLIC CURVE. 

THE ESSENTIAL DIFFERENCE BETWEEN THIS PROGRAM AND THE 
TeBsMs PROGRAM IS IN THE SEARCH TECHNIQUE. THIS 
PROGRAM USES THE FIBONACCI SEARCH METHOD WHILE THE 
ORIGINAL USED AN INTERVAL HALVING PROCESS» FOR DETAILS 
OF THE METHOD THE MANUAL ISSUED BY Is«BeMse CAN BE 
CONSULTED. 


DIMENSION P(100) sT(100) »xX(100) ¥(100) sV(100) s@Q@(2100) 
DIMENSION IFIB(40) 
L PORMAT(4012) 
2 FORMAT(8F10.3) 
600 FORMAT(6(/) s2l1Xs'BEST FIT CURVE FOR HYPERBOLIC! 
Pe eepee. wives 7) 
Ol FORMAT (6(/) sZixs' BEST FIT CURVE FOR EXPONENTIAL! 
Le DECI NE? s/) 
602 FORMAT (24Xs!STANDARD raed eae ee # pS oF Oe ay) 
BUS BORMAT(2Z7Xs'VALUE OF Cl IS. "sFlOse2s//527X 
er VeUCe (OF G2 LS her lose s7 } 
GUS PORMAT(27Xs' VALUE OF C3 IS 'sFlO0.237) 
605 FORMAT(17Xs'!X=VALUE ACTUAL COMPUTED i 
SUPER PERENCE Spe oX es  YVALUEl es SoX se "YeVALUE 67) 
606 FORMAT UR XsF6els6XsFi«eOs5XaF7sO0s6XsF7e0) 
607 iat 
608 FORMAT('1t 


THE INPUT DATA IS AS FOLLOWS 
NUMBR=NUMBER OF DATA POINTS 


TOPTISCHOICE OF CURVE 
L=EXPONENTIAL 


22BOTH 
BS=HYPERBOLIC 


IOPTZ=DATA SPACING 

O=EQUAL SPACING 

L=UNEQUAL SPACING 
TOPT3S=NUMBER OF ITERATIONS 

20 ITERATIONS ARE SUGGESTEDs 


TET THE DATA ARE EQUALLY SPACED ONLY THE LIMITS OF THE 
X=VALUES NEED BE GIVENe OTHERWISEs ALL THE POINTS 
MUST BE GIVENe THE Y-VALUES AREs OF COURSEs READ IN 
IN BOTH CASESe 
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C MAIN Pee ay 


READ(5 91)NUMSR»IOPT1»sIOPT2s10PT3 
IF(TOPT2eNEeQ) GO Lea) 3 
READ(592) TINsTOUT 
DELT=TOUT=TIN 
NUMBR1=NUMBR=] 
DELT=DELT/FLOAT(NUMBR1) 
DO 4 T=lsNUMBR 
TUL)=DELT#FLOAT( I=1)4+TIN 
ne rons 
Reaeivoee) CF CRS 
READ(552) (PUI 
SG a04 
FN=FLOAT(NUMBR) 
DO 6 [=l»sNUMBR 
TUTIGALCG (PED) } 
SP=SP+Y (1) 
GOUTO (7075100) TOPTI 
REXP= 


SPT=06 

DO 8 I=l1»sNUMBR 

Sr=STti typ) 
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Rmaere ese] =ST*ST 
Sp 


SGleaxPi(SST*SP=ST*SPT)/RRI) 
Bgatm ist PN*SPT<ST*SP)/RR1) 


Do ¢ =liaNUMBR 

=P{T }—CLl*EXP(=T(1)/C 
REXP=REXP+Q*Q 
S2ene (2009100) sITOPTI 
A2=0+05 
J=1 
AA2Z=O6 
AA3=06 
RR2=1eEt50 
RR3=leE+50 

RHYP=O6 
ST=06 
SST=06 
SPT=06 
DO 102 IT=l1»sNUMBR 
KCI )SALOG(1e+A2*T(I)) 
SES ree tT} 
SSTH=SST+X(1)*X(1) 
SPT=SPT+X(1T)*Y¥(1) 
Q=FN*SST=ST#ST 
Al=(SST*SP=$ST*SPT)/Q 
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S MAIN ess (CONT'D) 


A3=(FN¥SPT=ST#SP)/Q 

DO 103 IT=lsNUMBR 
Q=Y(T)—-Al=$A3*X(1) 
RHYPH=RHYP+Q*Q 

GO TO (10555009200) sy 
RRI=RR2 

RR2=RR3 

RRZ=RHYP 

AAL=AA2 

AA2Z=AA3 

AA3=A2 

TF (RR2“-RR3)500581045104 
A2Z=A2+A2 

Se 1O 101 

Pet J@NEwl) GO TO 502 
IFIB(1)=1 
TPIBtzZh=h 

DO 501 1=3rI0PT3 
IFI8(I)=IFIB( I<-1)+ 
FIB=FPEOAT(IFIB(IO 
J=Z 

JJ=] 

GO TO (59025913592 55935594) suJ 
RANGE=AA3=—AA1 
DEL=RANGE*FLOAT(IFIS(IOPT3<2))/FIB 
A2=AA1+DEL 
X1l=A2 

JJ=2 

GO FO 10} 
FXE=RHYP 
A2=AA3—~DEL 
Y1=A2 

JU=: 

GO TO.1L01 
PYI=REYP 
=[TOPT3=3 
KKK=IOPT3—2 
BOUSZ0  if=1 
LL=KKK=I1 
DELZERANGE*FLOAT(IFIBI(LL))/FIB 
PW wt Gel ) DEL=DELXO.6 9 

Pelee isGlerF¥1) GO 70 5223 
AA3=Y1 

Y1=xX1l 

FY1l=FX1 

KL=AALSDEL 

A2=X1l 

JIJ=4 

eo. 7O. 202 

FX1=RHYP 
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c MAIN ese(CONT!D) 


merTe S20 

AA1=X1 

Al=Y1 

FRIEFYI 

YIL=AA3=DEL 

A2=Y1l 

JJ=5 

SO FOC L102 

FYL=RHYP 

CONTINUE 

PPUPXLeGTseFY1l) GO TO 536 
A2=X1 

J=3 

Geto (LOL 

A2=Yl 
J=3 
GO-TO 10 
oes ee 
ITF (RHYP 
RR1=Os 
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DO 205 I=l»sNUMBR 

KOUNT=KOUNT +1 

CicouUN TRL ES i) “GO. TO 206 

WRITE (69607) 

WRITE(69605) 

KOUNT=10 

KOUNT=10 

WRITE(69606) T(I)sP( IT) sVI1) sQQ(T) 
CONTINUE 

STOP 

WRITE(69608) 

WRITE(62600) 

CL=EXP(AL) 

eae Li e/ oS 

C2=C3/Az2 

RR1=O6 

DO 207 I=l»sNUMBR 

WEL eCcl# (Let BHT (1) /C21** (ele) 7Ca) 
@Ot Ty =eP(ijeVil) 


7 Suu THa dite J 
» 457% ae 


- yy’ 1g 


7 : i L 
he Los at & ah 
To Lex -$06+ffEh OF Go aes 
e al be i ave aT de W993 | S08 


ae Gath TOs 
(etn 4)3T1F 

; : -)4660 6) SER 
Night 205 ~~ 

G4 +pOxarisettay 


| ey 7 (ii Ve S20) Se eee 
ie kd ‘AT OOnL PAAR 
. if 4) Lo aban 
- + ro) (Soeeeraiae 
ved d= (SO aeainen 4s 
7 i 


Oded, ” | 
Deas, a i 
viet, Soe — 
% 4 . its JCAaTI 

a: BD ie Ly vp 


ae)" Sow ee ee) 


MAIN eee (CONT'D) 


. 207 RR1= =RR1+QQ(1)*QQ(1) 
- RRIS(RRI/FN) #¥005 
WRITE (63602) RR1 
WRITE (63603) Clsc2 
WRITE (69604) C3 
WRITE (69605) 
KOUNT=19 
DO 208 I=1sNUMBR 
KOUNT=KOUNT +1 
IF (KOUNTeLE«57) GO TO 209 
WRITE (69607) 
WRITE (69605) 
KOUNT=10 
209 WRITE(65606) Té(I)sP(I) »Vil)sQart) 
208 CONTINUE 
STOP 
END 
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BEST FIT CURVE FOR HYPERBOLIC DECLINE 


STANDARD DEVIATION IS STelb 

VALUEDOF Cl 7S 2, 238687 

VALUE” OF C2 15 L0e20 

VALUE’ OF C3 15 Oe3l 

X=VALUE ACTUAL COMPUTED | DIFFERENCE 
YTVALUE YeVALUE 

15.40 5246 60956 14D. 
1620 58406 Sul Oee Love 
Ex eO 54406 53446 966 
1820 50806 5014s 566 
1940 47506 4/106 406 
2000) 44396 4430 Ge 
2160 41706 41726 zee 
2340 2680.6 Sills =a 
2400 324706 3206s “368 
2500 32/70% S3Le a “46 
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27020 29106 2I13B6 bal oo | 
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X=VALUE ACTUAL 

Y=VALUE 

5340 5953's 

54.0 918-6 

5509 884. 

Dae. % 8536 

5720 823% 

5840 T94e 

5920 7666 

60.0 739% 

6160 T14e 

6220 6396 

63460 6658 

6440 6426 

6540 6206 

66.0 5996 

a 2 6720 5786 
: 66%0 Botte 
69.0 5384 

7040 520 
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COMPUTED 
Y=VALUE 


9360 
9026 
8706 
B39. 
BOG. 
7T8le 
754. 
T29 
7046 
681. 
656. 
6376 
616. 
5966 
Sis 
2095 
5426 
5256 
509% 
494. 
4796 
4656 
4516 


Appendix C 


Program and Results of Integral Fit 


2 xfbneqgA 
$t3 [sipodnl to 2dfues bins mapovd 


oO 


CV EY OVO OF CYC OY GY FY CCA CY CVO MO ay 


Cc MAIN 
MAIN 


THIS PROGRAM WAS WRITTEN IN ORDER TO FIND A POLYNOMIAL 
APPROXIMATION FOR THE INTEGRAL OF THE SQUARE ROOT OF 
DULLA(D) (EVALUATED BETWEEN D=0 AND D=D). THE STANDARD 
LEAST SQUARES TECHNIQUE IS USED AND PAE Ce | Oe Tete 

FIT OF POLYNOMIALS FROM DEGREE 2 THROUGH 5 WAS TESTED. 
THE ‘SOLUTION OF THE SIMULTANEOUS EQUATIONSs FROM WHICH 
THE CONSTANTS OF THE POLYNOMIALS ARE FOUNDs WAS 
ACHIEVED USING THE GAUSS=SEIDEL TECHNIQUE. 


THERE IS NO DATA INPUT 
INTERNAL VARIABLES 


Di = VECTOR OF D VALUES WITH WHICH THE VARIOUS VALUES 
OF THE INTEGRAL WERE EVALUATED AND FITTED. 

DINICJ) = VALUE OF THE INTEGRAL CORRESPONDING TO Ditty) 

D2 =~ VECTOR OF D VALUES WITH WHICH THE ACCURACY OF THE 
Pit Who TESTED. 

DINZ{(J) = VALUE OF THE INTEGRAL CORRESPONDING TO D2tu) 


DIMENSION 01(30)%D2(30) sDIN1{(30) sDIh he sA(10) 
GO0 FORMAT('L' s6(/)s31lxXs! Nala: OF FIT 18S 'ehes7} 
601 LOE Salen OF Alslis! [TS ur2gaya7) 
602 FORMAT ( a ACTUAL COMPUTED ' 


Ls DIFFEREN he/ 29K oe Y=VALUE eS Aaa ye V ee te) 
FORMAT {/ 124% 4° STANDARD Pegi tie LS ele or} 


603 FORY erm 10s 3X eF Le 7 pak FL oye) 
GiO FORMAT('1! 


=O 
DINI(1)=0Q+6 
DIN2(1)=0+6 
DO J=2221 
Pete el yal) eke 


in ate gL tw a 

ata) =~ SIM POOL th) 
DIN2(J)=SIM1(D2(0)) 

See Ea 

DO 3 N=2Z95 

Cai L BEST(DIN NL»sD1ls GsNoA) 

NMAX=N+1 

WRITE(6 9600) 

DO 4 J=l»NMAX 

WRITE(62601) UsAlJ) 

SUM=0 6 

WRITE(6 9602) 
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és MAIN 
DO 5 Ja=ls2]1 
ANS=A (N41) 
DO 6 K=lsN 
L=N+ l= 
ANS=ANS*D2(U)4A(L) 


Y=DIN2(J)=ANS 
SUM=SUM+Y#Y 


eee (CONT tD) 


WRITE16 5603) 02(0) sDIN2(U) sANSsY 


CONTINUE 
SUM=(SUM/216) #65 
WRITE(6 9604) Sum 
WRITE(69610) 
CONTINUE 

STOP 

END 
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FUNCTION SIM1 
FUNCTION SIML{DULL) 


THIS FUNCTION SUBPROGRAM INTEGRATES THE SQUARE ROOT OF 
DULLA(D) BETWEEN 0 AND D USING SIMPSON'S RULE. 


DULLAH(D)=(160+(6604+0928125*D)*D) ¥#005 

PEC DULL Jl s1s2 

Y N2=2* TF IX (DULL /0*«005) 

H=DULL/FLOAT(N2) 

NU=N2/2—1 

S=1e+4e*DULLAH(H) 

DO 3 T=lsNU 

S=2e#DULLAH(H*FLOAT(2*1) )+4e*DULLAH(H*¥FLOAT(2*I4+1) )+S 
SIML=(S+DULLAH(DULL) )*H/3. 

GO TO 4 

SIM1=O06 
RETURN 

END 
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40 


60 
30 


80 
70 
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SUBROUTINE 3EST 
SUBROUTINE BEST(YPOINsPOINsusNoA) 


THIS SUBPROGRAM CALCULATESs FROM THE INPUT DATA» THE 
ELEMENTS OF THE MATRIX AA AND THE COLUMN VECTOR ¢€ IN 
ORDER THAT THE EQUATION (AA) (A)=C MAY BE SOLVED FOR 
THE COLUMN VECTOR As 


DIMENSION YPOIN( 30) sPOIN(30) sA(10)5SUM(30) sSUMY(30) 
DIMENSION €C(10)sAA(10910) »¥(30) 

M=e2*N+1 

NMAX=N+1 


DO 60 L=1s™ 
SUM(LY=SUM(L)+Y(L) 
SUMY GLI SSUMY (LI+FYCLI*XYPOIN(K) 
Y(LtL)ISY(L)*POIN(K) 
CONTINUE 

CONTINUE 

DO 70 K=l»NMAX 

DO 80KK=1 »NMAX 
Lak+KK—] 
AA(KsKK)=SUM(L) 
CONTINUE 

DO 90 Kl sNMAX 
CiKI=SUMY (K) 

CALL GASE(AsAAsCsaNMAX) 
RETURN 

END 
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SUBROUTINE GASE 


SUBROUTINE GASE(XsAsBoN) 


x 
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THIS SUBPROGRAM CALCULATES THE SOLUTION TO N 
SIMULTANEOUS EQUATIONS BY THE GAUSS=<SEIDEL TECHNIQUE. 


DIMENSION A(10910)sB(10) 9X10) 
DO 50 IT=loN 

SO X(T) 20. 

99 BIG=0+ 
DO 100 I=29N 

102 SUM=0.« 


oe CiveGel) GO TO 165 
LAST=I-1 
DO 106 J=l»sLAST 
106 SUM=SUM+A(T9U)*X(U) 
IFA lweQeN? GO TO 103 
ooo IND PLS) 41 
DO 107 J=INITLsN 
107 SUM=SUM+A(I9U)*xX (U) 
OS TEMP=(BV1)-SUM)/ACT oT) 
IFCABSC(TEMP=X(1))/TEMP) «GTeBIG) 
LBLG=ABSt (TEMP=X(1))/TEMP) 
LOO X(1I)=TEMP 
DEVS iGelLtecek@07) RETURN 
GO TO 99 
END 
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A@VALUE 


020 

Oe1l00 
O+«200 
00300 
02400 
Q0«500 
04600 
02700 
0+800 
06900 
La000 
1.6lo0 
12200 
16300 
16400 
16500 
16600 
i waene) 
18800 
12900 
24000 


ORDER (OF Fit (he «4 


VALUE OF Al 
VALUE OF A2 
VALUE OF A3 
VALUE OF A4 


VALUE OF A5 


ACTUAL 
eV ALOE 


0.0 

O.1l138904 
022523586 
064122117 
025914192 
Oe (885516 
120025415 
12325583 
144779329 
1273281144 
2¢0126419 
263011255 
2206032276 
209186506 
Be2471418 
365884542 
329423857 
4.3087692 
466874447 
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Program and Results of the Constants Calculations 
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THE FUNCTION OF THIS PROGRAM IS TO TAKE DATA FROM 

A WELL PREVIOUSLY DRILLED IN THE AREA AND CALCULATE 
THE CONSTANTS WHICH ARE REQUIRED FOR THE GALLE AND 
WOODS EQUATIONS AS DESCRIBED IN SECTION III-B 


THE INPUT VARIABLES ARE AS FOLLOWS 


N = NUMBER OF STAGES IN BIT RUN 
WWtS) = WEIGHT USED AT STAGE J 
Re va = SPEED USED “AT STAGE "J 
beuriuy = DISTANCE DRILLED AT STAGE 
Petit = TIME TO<pRIEE "STAGE 2 
SDEPTH = DEPTH AT START OF BIT RUN 
TDEPTH = DEPTH AT END OF BIT RUN 
DF - FINAL TOOTH WEAR 
BXF =~ FINAL BEARING Wi 
TE™= FOTAL ROTATING TIME 
Be =O 7OT AL DISTANCE DRILLED 
FORM = FORMATION RESPONSE 
Gx iF CSO 7 
be TF HARD 
FORMK = EXPONENT GN WEIGHT IN DRILLING RATE EQUATION 
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DULBAD(D he] (1 e0+(30+63093 /5*D)*D)#D 
DULLAH(D)=(160+(6s0+e928125*D)*D) #X005 
SIML(D) = (1601977946 10292304+4( 61133279*D—04127675)*D) 
¥H#¥D ) 4D 
KOUNT=10 
WRITE (69600 
WRITE(69601) 
20 FORMAT(12s8Xsl4F5ele/sl6FS5ele/slOF5el) 
21 FORMAT(2F10e3912F501) 
99 READ(5920) Ne (WW( J) »RR(U) sDELF( J) sDELT(U) su=lsN) 
fFt{NeEQeO) STOP 
TERR=0 
AF=O+« 
CON1=06 
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Cc MAIN o486 (CONT'D) 


CONZ=O0s 
CON3=0. 
FEET=06 
READ(5s21) SDEPTHsTDEPTHsDF s8XF 9 TF» FF »FORMsFORMK sHOLE 
%—5fP 
IP=IFIX(P/0e9)41 
DO 50 J=l5N 
WETIGHT=Ww( uy) 
REVS=RR(uU) 
A=EXP(=1l00e/REVS*#2 ) 
Per ORM—Os Ss) > Lelse 
ARR=A*REVS##O6 7540 e5*REVS#(1e=A 
Ge 1) 3 
es ARR=A*REVS##0 642840 e2*REVS*(1e—A) 
3 WBAR=7488#*WEIGHT/HOLE 
EMWATE=1359«el=/14e19*ALOGIO(WBAR) 
ROTI=REVS#*(letsO00004348*REVS##2) 
XLOAD=21338e87/(leteOZLIS7*WRAR) #¥#3.42258 
X1L=REVS/XLOAD 
XZ2=ROTI/EMWATE 
X3=(WBAR#X*#FORMK ) ARR 
PEE Mere. +PELEt J} 
CONS=CONS+DELT(U)*X1 
PPL IPebQeS) GO TO 30 
CONL=CONIT+X2*DELT (u) 
CONZ=CON2Z+X2*#DELF(I)/X3 
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30 CON2Z=CON2+DELF(U)/(X2*¥DELT(U)) 
50 CONTINUE! 
DF S=CON3/BXF 
GO TO (60%60570)sIP 
60 AF=CONL/DULLAI (DF) 
IF(IPeEQe1) CF=CON2/(AF¥SIM1(DF)) 
IF (IPeEQe2) CF=CON2/(AF*DF) 
GO TO 80 
70 CF=CON2/FLOATIN) 
O IF C(ABS(TDEPTH=(SDEPTH+FEET))sGTele0) IERR=1 
IF(IERRsEQel) WRITE(69603) 
IF (FORMeLT e045) TuUK=1 
TFUFORMeGTsOe>) [aK=2 
WRITE (6:602) SDEPTHs TDEPTHsAF 9CF sDFSsRESP(IUK) » FORM 
* yWEAR(IP) . 
KOUNT=KOUNT+2 
FE CKOUNT«LEsS2) GO TO 90 
WRITE(69600) 
WRITE (69601) 
KOUNT=10 
90 GO TO 99 
600 FORMAT('1's5(/)) 
601 FORMAT(14X9! INTERVAL AF CF DFS SPEED 
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* WEIGHT! s 
‘ WEAR! 9/945Xs'RESPONSE RESPONSE TYPE's/) 
602 FORMAT(LIXsF 600s '=' oFGeOe 1X sF6e2s1XsFGeds 1X sF 5029 3X 9AG 
*96X9F4e295X 
Le Ae ) 
COS FORMAT (//s22Xs'POSSIBLE ERROR's/s20Xs'DEPTHS 
* INCONSISTENT! ,/7}) 
END 
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EXAMPLE OF 


INTERVAL 
42800«-13000. 


T2800 e= 130006 
12000+«¢-12320. 


THE PRINT=OUT FROM THE CONSTANTS PROGRAM 
AF CF S SPEED WEIGHT WEAR 
RESPONSE RESPONSE TYPE 

Da f2 O4«0471 2003 HAIR D le ROL, ae 

Os U Je 0500 ZeVU HARD 1«0U a 

4e82 060504 1699 HARD LsQ So 


Appendix E 


Mobil O11 Test Data 


3 xfbnaqaA 


6360 sesT FtO frdom 


Sb eelzZe 
(INCHES) 


oe 750 
AF a 
5200 020450 
AF Cr 
4e50 0.40400 
AF ar 
4000 0560350 
AF core 
4050 040400 
ae 5 Sa 
(INCHES) 
827/50 


AF CF 


2000 040400 


ae St2E 
(INCHES) 
80750 


GROUP NOs id 
BIT COST TRIP RATE CHANGE TIME WEAR 
(HOURS ) (HOURS) (HOURS) TYPE 
1e860 04650/1000FT Par pon @) S$ 
INTERVAL LG4006 = 2920. 
> SPEED WEIGHT CONSTRAINTS 
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RESPONSE RESPONSE WMT WMAX RMIN  RMAX 
180 SOFT 1aU0 306 bUe 4Ue 1406 
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S SP BED WE LGHT CONSTRAINTS 
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INTERVAL 6810. = 8820. 


AF CF & SPEED WE TGHT CONSTRAINTS 
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GROUP NOe 4 


BIT SIZE BIT COST TRIP RATE CHANGE TIME WEAR 
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AF CF S SPEED WE IGHT CONSTRAINTS 
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RESPONSE RESPONSE WMIN 
2600 HARD 1200 306 


INTERVAL Lagiée. = 


5 SPEED WEIGHT 
RESPONSE RESPONSE WMIN 
200 HARD Lee OD BU « 


CONSTRAINTS 


WMAX 
456 


CONSTRAINTS 


WIhiAA 
{De 


CONSTRAINTS 


WIMAX 


4b5e 


CONSTRAINTS 


WMAA 
45s 


i [Mj N 


356 


aw 
RMI 


356 


RMIN 


354 


RMIN 
356 


RMAK 
71Vs6 


RIMAX 
tee 


RMAX 
T0e 


RIAA 
Tue 


Appendix F 


Program and Results of "Point" Optimization 
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C MAIN 
MAIN 


THIS PROGRAM OPTIMISES THE DRILLING PROCEDURE USING 
THE METHOD OF 4 "POINT?! OPTIMIZATIONe 


DIMENSION HOW(4) 
DIMENSION RESP(2) sWEAR(3) 
COMMON AF sCF sDFSsFORMsFORMK sPe ERR 
COMMON HOLE *BITHR»TRIPT »sFOOTIN» FEET sDOUT sBOUTsREPTIM 
COMMON D N»eBINs TP sXBAR sDELTAT sDELTAFsFIBsIFI8(30) 
COMMON INDICLIsINDIC2Z,sINDIC3sAINCi1 sAINC2 
COMMON XLIMLO(2) sXLIMHI (2) sCORRI2) 
COMMON WBESTsRBEST sDETBES »DEFBES sDBEST »8BEST 
COMMON FOOTLsTIMELsIOPTS 
DATA WEAR/ SS UstPrctgttTet 7 
DATA RESP/'SOFT!' s'HARDI/ 
DATA HOW/! i Pa ED eee ig: joe 

Ll FORMAT(8F10.3) 

2 FORMAT(4012) 

3 FORMAT(SF10e356F5el1) 
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4g! (HOURS) (HOURS } (HOURS } TYPE! afi 

GUS BORMATLTSSxXs'P RO €C EDU R-Efs/7 3 

6002 FORMAT(2XsF6e0s! “=! oF GeO lXsFGe2slXsF6e4s1XsF5e2 93% 
LAG a SX aF4e2s5XsFBe0e BX oF oe s2X oF 400s 2XaF4e0s3XoF beds 
LoOXsF6e3 92XeF4eZ2s5!/LOCOFT! s4XsFSe2s6XsA2s///) 

6020 FORMAT(27Xs'BIT NOs INTERVAL WELGHT  -SPEED t 
Lg¥'ROTs® TIME TRIP Time -EQUIVs GOST WEAR 4/s87Xs 
2'(RIG HOURS) tee 

S005 FORMATI(29XsT2s2XaFlelat=—!sF7els2XsFSe0s4XaF4e053X 
LFV e2e4X oF Tela SXsF7e2s5Xallse2xXvils//} 


6007 FORMAT(39Xs'FORMATION ROTATING TIME = 'oeFTe2s 
i" RIG HOURS's//s 
1 39Xs'FORMATION TRIP TIME = 'gF7e2ys 
L* RIG HOURS! o/7/% 
1 29Xe'FORMATION RIG TIME = 'sF[e2s 
1! RIG HOURS!'s//s 
2 39Xs'FORMATION COST = 'sFTle2s 
1* RIG HOURS’ s// 4 
2 329Xe' FORMATION BIT COUNT = '53XsI125//) 
6008 FORMAT(39Xs' CUMULATIVE ROTATING TIME = "sFfa2y 
1! RIG HOURS'»s//s 
Bi ZOXs "CUMULATIVE TRIP TIME = "gr T7eZs 
1! RIG HOURS! »+//s 
1 29X%e' CUMULATIVE: RIG TIME = loFT7e2s 
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C MAIN eee (CONT'D) 
Ti ORIG HOURS! s//y 
2 Z9Xv' CUMULATIVE CosT = ar pee e 
LPERIG"HOURS! »//s 
2 39Xs' CUMULATIVE BIT COUNT = Te 3Xaten/7)} 


THE INPUT VARIABLES ARE AS FOLLOWS 

ERR = STEP SIZE AT WHICH SEARCH 1S STOPPED = TRY Oel 

AINC] = WEIGHT INCREMENT FOR CALCULATING DERIVATIVE 
TRY O85 


AINC2Z = SPEED INCREMENT FOR CALCULATING THE DERIVATIVE 


CORR(1) = WEIGHTING FACTOR = TRY le 
CORR(2) = WEIGHTING FACTOR = TRY a. 
NEST =" NOs OF SECTIONS IN GROUP 
NOSE? = SWITCH — ZERC UNTIL THE LAST DATA SET WHEN A 
NON=ZERO NUMEER IS NEEDED 
Witwee = (MOLE SP7F 
Siicsie= BIT cCOsT 
REGGS| = RIG COST 
TRIPT = TRIP RATE = TIME PER RETURN FOOT 
Ree rime CONNECTION COST 
P = WEAR TYPE 
Celif SS 
Ce one 4 ore 
Er ie Lic 
SUEPTH = DEPTH AT TOP OF SECTION 
TDEPTH =-DERTH AT BOTTOM OF SECTION 
BR t=eCONSTANT IN DULLING EQUATION 
CF = CONSTANT IN THE DRILLING RATE EQUATION 
DFS = CONSTANT IN THE BEARING WEAR EQUATION 
FORM = RESPONSE TO SPEEDs Qe IF SOFTs ls IF HARD 
FORMK = EXPONENT ON WEIGHT 
ALIMLO( ) = LOWER LIMITS ON WEIGHT AND SPEED 
MLIMHT( ) =~ UPPER LIMITS ON WEIGHT AND SPEED 


IMPORTANT INTERNAL VARIABLES ARE AS FOLLOWS 


DIN = INITIAL TOOTH WEAR 
DOUT = FINAL TOOTH WEAR 

BIN = INITIAL BEARING WEAR 

BOUT = FINAL BEARING WEAR 

1CBIT = CUMULATIVE NUMBER OF BITS 

CTROT = CUMULATIVE ROTATING TIME 

GTRIP.=- CUMULATIVE TRIP TIME 

CEQ = CUMULATIVE COST 

STROT = ROTATING TIME FOR SECTIONS IN A GROUP 
STRIP = TRIP TIME FOR SECTIONS IN A GROUP 
STEQ = COST OF SECTIONS IN A GROUP 
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C MAIN eee (CONT'D) 


ICBIT=0 
CTROT=0. 


94) ERRSATNCIsAINC2s (CORR(U) sUuel92) 
PR PBL) =] 
=) 


READ S62) 1OPT& 

Do eks. JS 5 TOPTS 

IFIB(J)=IFIB(U=1)+IFIB(y-2) 

PIB=FLOAT¢CIFIB(IOPT4) } 

READ(592) NINTsNUSET 

STROT=0+6 

STRIP=06 

Suicenar he 

DIN=0 «6 

BIN=06 

INDIC3=0 

NBIT=0 

FOOTL=A06 

TIME1=O6 

DO 10 JJ=l»eNINT 

READ(S5s1) HOLEsBITCST»sRIGCSTsTRIPTsREPTIMsP 
IP=IFIX(P/009)+1 

XBAR=7e98/HOLE 

BITHR=BITCST/RIGCST 

TRIPHRETRIPT#1000.6 

READ(593) SDEPTHsTDEPTH»sAF »sCF  »DFES sFORMsFORMKs (XLIMLO 
Lisi Je lw is (XLIMHI (2) ved 2) 

IJK=IFIX(1e¢6+FOR™M) 

WRITE(6 6000) 

WRITE (696001) 

WRITE (696002) SDEPTHsTDEPTHsAF sCF sDFSsRESP( TUK) sFORMKs 
LXLIMLO(1) sXLIMHI(1) sALIMLO(2) sXLIMHI(2) sHOLE®SBITHRs 
LTRIPHRsREPTIMesWEAR(IF) 


WRITE (696003) . 
WRITE(6 96020) Pr 
FEET=SDEPTH 


CALL ASCENT 
FTEMP=FEET+DEFBES 

IF (FTEMP*GTsTDEPTH) GO TO 20 

IF(INDIC3eEQel) TRIP=(TRIPT*FTEMP+REPTIM) /2« 

IF (INDIC3eEQ.60) TRIP=TRIPT*(FEET+FTEMP) /2e+REPTIM 
TEQ=DETBES+TRIP+BITHR 

STROT=STROT+DETBES 

STRIPSSTRIP+TRIP 

STEQ=STEQ+TEO 
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C MAIN seal CONT'D) 


ID=IFIX(8e*DBEST+0, 4 

IB ety esau 

SRIG=STROT+STRIP 

IF CINDIC3« ie NBIT=NBIT+41 

WRITE(6 96006) NBITsFEETsFTEMP sWBEST»RBEST»DETBES» 
1TRIPsTEQsID v1B 
FEET=FTEMP 

Da =O s 

BIN=06 

INDIC3=0 

FOOT1=0.6 

TIME1=0+% 

io VY 21 

FOOTIN=TDEPTH=FEET 

INDIC2=1 

CALL TEST(WBESTsRBEST»DCOST) 

DBEST=DOUT 

BBEST=BOUT 

DETBES=DELTAT 

FTEMP=TDEPTH 

TRIP=Ols 

IF(INDIC3eEQe0) TRIPS(TRIPT*FEET+REPTIM) /26 

IF(JJsEQeNINT) TRIP=TRIP+(TRIPT*#TDEPTH+REPTIM) /20 
TEQ=DETBES+TRIP+BITHR 

STROT=STROT+DETBES 

STRIP=STRIP+TRIP 

STEQ=STEQ#TEO 

ID=IFIX(8s6*DBEST+0-5) 

IB=IFIX"4 8se*BBEST+0.5) 

SRIG=STROT+STRIP 

IFC INDIC3ZceEQe0) NBIT=NBIT#1 
WRITE(626006) NBITs»FEETsFTEMP»WBEST »sRBESTsDETBES>s 
LTRIPsTEQsIDs IB 

EMD LC 3 at 

DIN=DBE ST 

BIN=BBEST 

POOTL=FOOTLS+FOOTIN 

ime i= | LMELSDETBES 

CONTINUE 

WRITE (696007) STROT®sSTRIPsSRIGsSTEQ®NBIT 
LCBIT=ICBIT+NBi7T 

CTROTH=CTROT+STROT 

CTRIP=H=CTRIP+STRIP 

CEG=CEQ+STEQ 

CRIG=CRIG+SRIG 

WRITE(656008) CTROTsCTRIPsCRIGsCEQsICBIT 
IFUNUSET»sEQs0) GO TO 99 

STOP 

END 
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SUBROUTINE ASCEN 
SUBROUTINE ASCENT 
THIS SUBROUTINE USES A GRADIENT SEARCH METHOD TO GET 
TRE OPTIMUM SPEED AND WEIGHT. NORMALLY THE SEARCH 
TERMINATES AFTER THE STEP SIZE BECOMES SMALLER THAN ERR 
BUT IF THE SEARCH REACHES A CORNER OF THE PERMITTED 
AREA IT IS TERMINATED FOR THAT REASON. THE SUBROUTINE 
ALSO CONDUCTS A SEARCH ON THE FINAL WEAR CONDITION 
IF THE BIT IS FOUND TO BE RATE CONTROLLED. THE START= 
ING VALUES FOR THE SEARCH ARE THE UPPER LIM l1ITS ON THE 
WEIGHT AND SPEED 
INTERNAL VARIABLES ARE AS FOLLOWS 
SS = STEP SIZE 
PBEST = BEST VALUE OF FUNCTION FOUND AT ANY STAGE 
BAS PUNGTION VALUE 
Nook! <2 DERIVATIVE WeRets WEILGHT 
MELTS Ye PERIVATIVE WaReTe SPEED 
Rl’ ) = BASE VALUES FROM WHICH STEPS ARE MADE 
Ke{ ) = POINTS AT WHICH FUNCTION IS EVALUATED 
DEIBES ~ BEST VALUE OF ROTATING TIME 
DEPSES<— —BEST.VALUE OF FOOTAGE 
DBEST =~ VALUE OF FINAL TOOTH WEAR ASSOCIATED WITH THE 
BEST WEIGHT ANG SPEED 
BBEST = VALUE OF FINAL BEARING WEAR ASSOCIATED WITH 
THE~BEST WEIGHT AND SPEED 


DIMENSION AL(2) 9X2 
LOGICAL AGB God 


Ve Oeics 


COMMON AF sCFsDFSsFORMsFORMK sPsERR 
COMMON HOLEsBITHR»sTRIPTsFOOTIN»s FEET sDOUT »BOUT »sREPTIM 
COMMON DINs8 [Ne TPsXBAR sDELTAT sDEL TAF sFIBsIFIB{(30) 
SOMMON INDICITs INDIC2sINDICSsAINCIsAINCZ 
GOMMON XLIMLO(2)sXLIMHI (2) sCORR(2) 
COMMON WBESTsRBESTsDETBESsDEFBES»DBEST sBBEST 
COMMON FOOTIsTIMELs IOST4 
INDIC2=0 
INDIC1= 
Li LVSEAL IMAI C1) 
X1L(2)=XLIMAI (2) 
SS=ls 
LL=0 
KL=0 
CALE *TESTURICL) sX1(2) sFBEST) 
COLO PERT VEX Olds KD C2) sXL IL ePBESTs 
ieee Geka rai. E-10 ANDsABS(XL1(2))ebLTelsE=10) 


1 GO TO 100 


GO J=ls2 
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ON 


30 


100 


SUBROUTINE 


bee TS a ae )=SS* 
PF EX2 61 eGTeXLI 
re tke Cul le Eeet 
eae Teo hing (2) 
IF (FXeLTeF BEST) 
[IF (SSeLTseERR) 
SS*8S/2 6 

ik Sd bh 

GO TO 5 
FREST=FX 
KLEKL+1 
Peat weGT sO) 
SS=SS*2.e 

DO FF Je 1s2 
X2(4)=X1(U)-SS 


GC 


Canon veat (K2 019 


frPoPA@LESPBEST) GO TO 6 

SS=S5/2 « 

DG <2 Jel s2 

ALASTEAL OI)—-SS*XLI(U)*CORR IU) 

ee ALIMAT eT) Adtoyvexi iM 

om et eLTeXL IM ions All I)=XLIM 

A=  bretand kad DA crores Ee ge 

BSABS CXL CZI-XLIMHI(2) bel Teel 

| Day bight GO TO 100 
S=ABS(AL(2)=XLIMLO(2))eLTeel 

TE ( AeANDe entice aes G 

Ce AS (el be KL TMLOC Ley elt. od 

IF(DeANDeB) GO TO 100 

SS=2e*SS/2e**KL 

LL=0 

KL =O 

GO TO 9 

etait ol ww Oe 1.) GO te aia 

Sean weet Ad bleh be er BEST? 

WREST=X1(1) 

RBEST=K 142) 

TSEST=FBREST 

DETBES=DELTAT 

DErBES=DELT AF 

DBEST=DOUT 

BBEST=BOUT 

PPtINDICLeEQs2) RETURN 

DIFF=DOUT=DIN 

IF (DIFF el Tse0206) RETURN 


COU T=DOUT=0605 
INDIC1=1 
KNT=0 


ASCENT 


ese | 


cee *# CORR 


MAT ( 


( a) 
MLO(O))} Ae Jii= 
9X2(2) 9FX) 

GY TO 6 


So ae oe 


TO 30 


SAL LU IPECORR ( J) 
IF(X2( 4) eGTeXLIMHI (uU)) 
> i ee i Nag On a | 


X2(4U)=XLI 


9X2(2)9FX) 


) X2(0)=xXLIN 


CONT'D) 
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SUBROUTINE ASCENT see(CONT! 


Oo TO. 10 
KNT=KNT +1 

IF CKNTeNEo 1) oO to 192 
IF(FRESTs«GEeTBEST) 
RANGE=DBEST=DIN 

AL=DIN 

A3=DBEST 


DEL=RANGE*FLOAT(IFIB(IOPT4=2))/F1IB 


DOUT=A1+DEL 

KX L=DOUT 

rH ie eS 

SOM ry 10 

Canto imslOsiliglies is issu 
PXIl=FBEST 
DOUT=A3-DEL 
KY1=DOUT 

Pe 4 

Git 20 
PL BE SF 

VK K=ITOPT4=3 
KKK=ITOPT4=2 

DO 120 I=1leMKK 
LL=kKKK=] 
DEL=RANGE#FLOAT(IF 
Pet LLekQal) DEL=DE 
Por adel eh hl i GO 
A3=xY1 

XY1=2XxXl 

PYL=FX1/ 
XX1LSA1+DEL 
DOUT=XX1 

JJ=3 

GO 710 10 

PRI=F BEST 

GO To 120 

Al=xXxX1 

AX1=XY1 

FX1=FY1 
AYASAA—-DEL 
DOUT=xXY1 

JJ=4 

GO 70 10 

PY TeF BREST 
CONTINUE 

INDIC1=2 

TFC FX1leGTeFY1) GO TO 136 
DOUT=XX1 

GO TO 10 

DOUT=XY1 

Go TO 10 
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SUBROUTINE ASCENT ++e(CONTID) 
—_ - | 
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eer Lignan ai 
SUBROUTINE TE STI WEIGHT 9 f Ne Vee es 1) 


THIS SUBPROGRAM USES THE GIVEN VALUE OF WEIGHT AND 
Seem PLUS THE INI PIAL WEAR CONDITIONS PRESENT AT THE 
TIME TO DO ONE OF TWO THI) ae CITHER fT CALCOLATES, THe 
PINAL WEAR CONDITIONS AND THE FOOTAGE: DRIULEDs THE TIME 
TAKEN AND THE COST PER FOOT) OR» [IF THE SWITCHES ARE 
oe) IN THE OTHER POSITION IT CALCULATES THE TIME,» 

COST AND FINAL WEAR CONDITIONS FOR A GIVEN FOOTAGE 


COMMON AP sCFsDFSsFORMsFORMK sPsERR 

COMMON HOLE sBITHRsTRIPT »sFOOTINs FEET sDOUTs BOUT» REPTIM 
COMMON DINsBINsIPsXBARs DELTAT sDELTAF sF IB» LFIB(30) 
COMMON INDICLsINDIC2sINDIC3sAINC1sAINC2 

COMMON XLIMLO(2) oXLINHI(2) sCORR(2} 

COMMON WBEST»RBEST»DETBESsDEFBES »DBEST »BBEST 
COMMON FOOTIsTIME1sIOPTS4 

DULLA(D)= 12604+(6604+.928125*D) #0 
DULLAT(D)=(1.5,0+(3e+.309375%D)#D)*D 
DULLAH(D)=(1.04+(6004+6928125%*D) #D) ¥¥#0.5 
SIM1(0)=(16019779+(16292304+(21133279%*D-064127675) 
LOD Pe) #D 

IP=IFIX(P/069)+1 

IF (AF eLTelsE-05) AF=les 

A=EXP (=100/REVS#*2) 

IF (FORM=O051 19192 

ARR=A#REVS##0 «7540 e5*REVS#(1e=A) 

GO To 3! 

ARR=A*REVS##06428+00e2*#REVS#(1Le—A) 
WBAR=7e88*WEIGHT/HOLE 
EMWATE=1359e¢1=714e19*ALOG1O(WBAR) 
ROTISREVS*(leteO0004348*REVS#*2) 
KLOAD=21338687/(1lete03118 (*WBAR ) ¥*¥#36.2258 
FISAF*EMWATE/ROTI 
F2Z=FL*¥CF* (WBAR*#FORMK ) ARR 
PO=FISREVS/DFS/XLOAD 
IF ( INDIC2e¢EQel) GO 7O Lid 
GO2TO (18518328) siP 
PFC EINDIC Ls NEO) GO: Gao di 
DOUT=l1l. 

BOUT=BIN+(DULLAI (DOUT)=DULLAI(DIN) )*F2 

IF (BOUT eLEel) GO TO) 5 

BOUT=1l-. 

C=DULLAI(IDIN) #(Lle@BIM) /F3 

KE le 

DOUT=X=(DULLAI(X)=-C)/ DULLA(X) 

TE CARS( (DOUT=X)/DOUT) aL Tel»sE=05) GO TO 5 
Yi wiow 7 

GO TO 4 
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SUBROUTINE TEST 


ieee 2) 


eee (CONT'D) 


DELTAT=F1*(DULLAI (DOUT)=DULLAI(DIN)) 


SO TO (51552528) le 
DELTAF=F2*(DOUT=DIN) 

GO TO 48 

DELTAF= F2%*(SIM1(DOUT)=SIM1(DIN)) 
FOOT2=FOOTI1+DELTAF 


DCOST= (B81 THR+REPTIM+TRIPT* (FEET+(DELTAF- =FOOTL)/20 4 


IF DELTAT+TIMEL) FOOT? 


RETURN 

CONTINUE 

mo TO (19529539) .,]p 

X=DIN+0 62 

C=SIM1(DIN)+FOOTIN/F2 
DOUT=X=(SIM1(X)=C) /DU NC 

IP {ABS ( (DOUT=xX)/DOUT) «eL Tslee=05 
=DOUT 

Go FO. °9 


DIFF=DULLAI (DOQUT)=<DULLAI{(DIN} 
BOUTH=BIN+F3*#DIFF 
DELTAT=FI1*DIFF 

cos Vis ih aes et dea 

DCOST= 
1S Ls FOOT 2 
RETURN 
FIl=DFS*XLOAD/REVS 
BOUT=POOTIN/(FLEC 
DELTAT=FLl4#(BOUT<$68 
DOUT=O06 

GO TO 4¢& 
DOUT=FOOTIN/F24DIN 
DDIFFEDULLAI (DOUT)=$DULLAI (DIN) 
BOUT=F3*DDIFF 
DELTAT=FL*¥DDIFF 
Oo 2442 
Peli tard (FSA lis 
DELTAF=F2/F3*(les 
DOUT=0.« 

GO TO 48 

END 
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RESULTS OF SIMPLE 
- CONSTANTS IMPLIES 


BIT NOs 


wn e 


BIT NO» 


OPO Fe 


BIT NOs 


uw SW Ph re 


BIT NOs 


wm SW NM be 


BIT NOQe 
1 


"POINT": METHOD. 


A NEW GROUP. 


GROUP NO.w 1 
INTERVAL 
1SG0s. = 2847, 
LGale = 2263. 
22630 — 2700. 
27006 +. 2920+ 


GROUP NOe 2 


INTERVAL 
gaeOe | = Nese Soe 
DO Sta er Siena 
SOT [MONS SY hae 


GROUP NOe 3 


INTERVAL 
37708 = 4086-6 
4086¢ = 4405. 
44056 =~ 4730. 
47308 = 5058s 
505868 = 5070. 


GROUP NOe 4, 


INTERVAL 
50700 = 5476s 
24766. > 2884.6 
5B84e = 629Te 


6297s = 67176 
Me vies 67706 


GROUP NOe 5 
INTERVAL 


GTIR Hh GBLOs 


A CHANGE JIN 
WEIGHT SPEED 
606 1506 
60s L506 
606 148.6 
600 146.5 
WEIGHT SPEED 
606 1406 
60 1368-6 
bUe vA Cia 
WEIGHT SPEED 
596 1326 
596 1326 
59° 1296 
596 1256 
596 L27e 
WEIGHT SPEED 
60 1326 
596 132. 
60's 1306 
596 1296 
59 «8 129.6 
WETGHT -SPEED 
50. 


60+ 


GROUP NO. 6 
BIT NOs INTERVAL WEIGHT SPEED : 
1 668101 = .71500 BB. 110. fa 
2 TL50a =~ (7491. 5B. 110. 
3 T4921 = 18346 5B« 110 
4 78346 = B1l80s 5768 110» 
Ls S180. = B527. bys LOG 6 
6 Soe 7s = BB2U0. 57. 1096 
GROUP NO. 7 
BIT NO. INTERVAL WEIGHT SPEED 
1 88206 >= BB6U. 506 46 
i 
GROUP NO. 28 
BIT NOs INTERVAL WEIGHT SPEED . 
1 S860e-= 9020% 586 100 
GROUP NOs 9 
BIT NOs INTERVAL WEIGHT SPEED 
ii 9020.6 -_ 9060.6 50-6 466 
} GROL?® NOs 10 
BIT NQe INTERVAL WEIGHT SPEED 
: 
L SO CO. = 9270 556 Lo 
GROUP NOs 11 7 : 
BIT NOs INTERVAL WEIGHT SPEED 7 
if 92706 = 99699 60. 90° i 
2 9969s (e° 20080s 160K 90» 
_ 7 
GROUP NOs 12 a 
ve 7 : 
BIT NOs INTERVAL WEIGHT. SPEED 


: ae 
100806 = 104014 456 70% 
104016 - LO7Z21e) «4be> (> 970% 
LO721« =: LCP RR Sade amas 


=) 4 ay i? af) 


BIT NOe 


wm & WN Ee 


Sit NOs 


wn 


CUMULATIVE 
CUMULATIVE 
CUMULATIVE 
CUMULATIVE 
CUMULATIVE 


ROTA 


GROUP NOs 


LO7506 
10869. 
109886 
as 
112266 


13 
INTERVAL WEIGHT 
ee! LO&69.6 45 6 
= 10988. 456 
= VELOT « OS 
— L1l2266 son Ve 
oe a SDae 46 
NOs 14 


GROLIP 


INTERVAL WEIGH 


Li3Z136 a 11680. 45 6 
GROUP NOe 15 

INTERVAL WEIGHT 
116806 = 11899. 45 6 
L1i899.« po 12084. 45 « 
GROUP NOe 16 

INTERVAL WEIGHT 
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ICBIT=0 
CTROT=0.% 
CTRIP=06 
CRIG=0. 
GEQ=0% 


READ(Ss2)IWlsIR1IsIOPT1»,IOPT4s1IOPTSs10PT6 

STROT=0% 

FOR INPUT DETAILS SEE PROGRAM ON 'RPOINT? OPT? 

ADDITIONAL INPUT IS AS FOLLOWS 

MUSET Lo Ir NOT ZERO. ANOTHER DATA SET ExXiSis 

NUSETZ ~ IF NOT ZERO ANOTHER DATA SET EXISTS 
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CONTINUE 


+RRANGE*XFL OAT (J=1) /FLOAT(IR=1) 

CONTINUE 
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CONTINUE 
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TID=IFIX(8e*DTEMP+0<5) 

TIB=ITFIX(8e#BTEMP+025) 

STROT=FLOAT(MBITI*#TTEMP 

STRIP=FLOAT (MBIT) *( (TOEPTH+SDEPTH) #TRIPT/2«e+REPT IM) 

SRIG=STROT+STRIP 

STEQ=SRIG+FLOAT (MBIT) *8ITHR 
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WRITE(6s6001) 
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: CTROTSCTROT+STROT 
CTRIP=CTRIP+STRIP 
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STOP 
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m= «2001 CONTINUE 
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AA2Z=O6 
AA3=06 
RRI=1leE+50 
RR2Z=]4eE+50 
AzZ=56 
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90 CONTINUE 
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L+FLOAT(MBIT) *#(REPTIM+81THR) 
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RRL=RR2 
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RRZ=COST 
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IF(RR2=-RR3) 50091049104 
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GO TO 93 
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RANGE=AA3=AA] 
DEL=RANGEXFLOAT(IFIS(IOPT4=2) )/FIB 
A2=AA1+DEL 
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JJ=2 
GO TO: 93 

FX1=COST 

AZ=AA3=DEL 

Y1=A2 
JJ=3 

GO 70.93 
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KKK= TOP T4<2 

DO 520 IT=lseMKK 

LLEKKK=I[ 

DEL RANGeer LOA (CT EIRWUiL) /EIB 
SCL weGiwl) DELS DELO. 9 
Dies OS Parra © iGO TO: S23 
AA3=Y1 

Y1l=xX1 

ie ae 

XL=AA1+DEL 
A2=X1 
JJ=4 
SO Oi 9 
FX1L=COS 
GO) 1G: 3 
AA1=X1 
Al=Y1 
FX1=FY1 
YL=AA3—DEL 
A2=Y1 
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Sn O. 93 

RY L=COST 
CONTINUE 
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WWClsl)=WWwl29u) 

RR(lsU)=RR(2s8¥U) 

TIMECLs J) =TIME (2s) 
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CONTINUE 
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CONTINUE 
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WRITE(6 36008) CTROTsCTRIPsCRIGsCEQsTICBIT 
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Vee rOURMATI? L's6(/)) 
END 


'eF4Gelsl0Xs' COST 


eee (CONT'D) 
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SUBROUTINE SEARC 


SUBROUTINE SEARCH (FOOTINS WREST sRBESTsTRESTs 
LDSEST+BBEST) 


THIS SUBROUTINE CALCULATES THE BEST WETGHT AND SPEED 
BY A GRID SEARCH METHOD. THE ASSOCIATED COST AND THE 
FINAL WEAR CONDITIONS ARE ALSO FOUND. 


C 
€ 
a 
c 


COMMON AF sCF sDFS »FORMsFORMK 

COMMON TDEPTH sSDEPTHsHOLE sRIGCST »sBITCSTsTRIPTsREPTIM 
COMMON NOLsNOMAX 5IW912 

COMMON KLO(2) »XH1 (2) sERRo W550) sR(1400) »P 
TBEST=1.E+10 

ITEMP=IW 

DO 13 K=lsIR 

IND=0 

IJ=ITEMP 

TFCITEMP2EQ<e0) GO TO 14 

DO 12 JdFlelyu 

LFREINDsE@el) GO.TO-13 

J=TUtl=ou 

CALL TEST(W(U) »R(K) sFOOTINS FX sDX 5 Bx) 
IF (FXeGTeleE+08) GO TO 15 

IND=] 

fr wpaecee BESh) GO TOckz2 

TBEST=FX 

RBEST=R(K) 

WBEST=Wil J) 

DBEST=D'x 

BREST=8xX 

GOL io k2 

ITEMP=J-1 

CONTINUE 

CONTINUE 

RETURN 

END 


SUBROUTINE TEST1 
SUBROUTINE TESTI (WEIGHT sREVSsFOOTIN) 


C THIS SUBROUTINE CALCULATES THE FOOTAGE POSSIBLE USING A 
C GIVEN WEIGHT AND SPEED. 


COMMON AF sCFsDFS»FORMsFORMK 
COMMON TDEPTH»SDEPTHsHOLE»sRIGCSTsBITCSTsTRIPT»REPTIM 
COMMON NOLsNOMAX» IW SIR 
COMMON XLO(2) sXHI(2) s9ERReW(550) s9R(1400) 5P 

DULLA(DI= 1460+(600+.928125#D) *D 
OU LAT 4D feet Le O+( Bete ZQOIATSXD) HD) #D 
DULLAH(D)=(160+(660+e928125*D)¥D) ¥#065 
SIM1(D)5(16019779+(10292304+(+61133279%D=04127675) 
1¥*D)*D)*D 

IPSIFIX(P/0.9)+1 

LP CAR sel TeleE=-05) AF=l. 

ASE Xe (=LOO06e/REVS**2) 

IF(FORM=065) Isls2 


rl ARR=A*REVS**O 6 75+0e5*REVS*( 16 =A) 
GO TO 3 
2 ARRS=A*REVS#*0 oe 42840 e2%*RPEVS#(Le=A) 


3 WBAR=/7s688#*WEIGHT/HOLE 

| EMWATE=135961=714619*ALOG1O(WBAR) 

3 ROTISREVS#(LleteO0004348*REVS¥##2) 
XLOAD=21338e87/(loetseO3 LIA 7#WRAR) ¥¥3 62258 
FL=AF*EMWATE/ROTI 
F2=F1*CF* (WRAR**FORMK)*#ARR 
F3=F1¥REVS/DFS/XLOAD 
GO TO (20520330) »IP 

| 20 BOUT=F3*DULLAI(1.) 


LPreouUrelealeO)- GO TC 10 
X=]. 

| & DMAX=XK=(FS*DULLAT(X) lel /(FSZ*DULLA(X)) 
DFC ABS( (OMAX=X) /DMAX) eLEeleE=05) GO TO 60 
X=DMAX 


Go TO 6 

10 DMAX=1e 

GO TO (4940230) s5IP 

4 FOOTIN=F2*SIMI1(OMAX) 

RETURN 

40 FOOTIN=F2*DMAX 
RETURN 

30 FOOTIN=F2/F3 

RETURN 

END 
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SUBROUTINE TEST 
SUBROUTINE TESTU JEITGHT sREVS sFOOTINs TEFF sDOUT «BOUT ) 


THIS SUBROUTINE CALCULATES THE TIME TO DRILL A GIVEN 
FOOTAGE USING THE GIVEN WEIGHT AND SPEED.e THE FINAL 
WEAR CONDITIONS ARE ALSO CALCULATED. 


oy 'ON AP s CFs DFS sFORMsFORMK 
CON Hee TDEPTH»SDEPTHsHOLE RIGCSTsBITCSTsTRIPTsREPTIM 
COMMON NOlsNOMAXsIWsIR 
OMMON XLOC2)eXH1 (2) sERROW (550) 9R(1400) sP 
DULLA(D)= 120+(600+6928125*D) xD 
DULLAT ( (D)=(160+(36+0299375*D) *D) ¥D 
DU peer 160+(60046928125%*D)#D) ¥¥0.65 
SIM1(D)=( 160197794 (10292304+( 61133279%*D=—s4127675) 
(SAD 1340 
IP=IFIX(P/069)+1 
A=EXP (=1l000e/REVS*#2) 
TFC FORM=025) lsls2 


i ARR=A#REVS##0 © 15+0e5#REVS#(1Lle=A) 
GO fO 3 
2 ARRFA*REVS#XO oe 42840 6 2*REVSH(Le™A) 


3 WBAR=7e88*WEIGHT/HOLE 
EMWATE=1359e1l=—714e19*ALOG1OI(WBAR) 
ROTISREVS*¥(1 ete QC U0U4G34R2#REYV SHR ) 
XLOAD=213386¢87/(1le+s031187*WBRAR) *¥#362258 
GO 710 (20u20930)+41P 

20 FI=AF*#EMWATE/ROT I 
F2=FL*EOF* (WRAR#*FORMK ) #ARR 
FS=FI*REVS/DFS/XLOAD 
GO TO (40950930) 51P 

40 C=FOOTIN/F2 
X20 

4 DOUTH=X—(SIM1(X)—-C) /DULLAH(X) 

IF(ABS( (DOUT=X) /DOUT) «eLTelsE-05) GO TO 5 
X=DOUT 
GO TO 4 

SB Tr( GOUT s*eGTele) GO TO 10 

DDIFF=DULLAI(DOUT) 
BOUT=F3A*DDIFF 
IF(BOUT.GTe1l.) GO TO 10 
terrer L*ODIFF 

S970. Lit 

10 TEFF=leE+10 

bi RETURN 

BQ DOUTSFOOTIN/F2 

FipouTsGTel.0) GO TO 10 
DDITFFSDULLAI(DOUT? 
BOUTH=F3*DDIFF 
Te(BOUTs*GTeleQ) GO TC 10 
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SUBROUTINE TEST ese (CONT'D) 


TEFFSFL*¥DDIFF 

RETURN 

F1l=DFS*XLOAD/REVS . 
BOUTSFOOTIN/( F1¥CF*(WBAR¥*FORMK) ARR) 
TF(BOUT+GT»1l.0) GO TO 10 

TEFF=F1l*BOUT 

RETURN 

END 


o* ieZ 


G MAIN 
C MAIN 
C THE PROGRAM USES THE DYNAMIC PROGRAMMING APPROACH 
G DESCRIBED IN SECTION III=D TO ALLOCATE THE DISTANCE 
2 TO BE DRIELED BY EACH BIT IN AN OPTIMAL MANNE! 
C IN MANY RESPECTS THIS seibotae [S SIMILAR TOoTHE «SIMPLE 
C INTERVAL PROGRAM AND SO ONLY THE Ua rete ee ey CW LE ie ek 
ei INDICATED 
& 
DIMENSION FINT(40) sWINT(40) sRINT (40) sTROT(40) sDINT 
1(40) eBINT(40) sTTRIP(40) sTEQ(40) sID(40) 5 IB(40) > 
LRESP(2) sWEAR(3) sFOUT(41) 
COMMON AF 9CF sDFSsFORMsFORMK sBITHR 
COMMON TDEPTHs SDEPTHsHOLEsRIGCSTsBITCSTsTRIPTsREPTIM 
COMMON Sige at ie coneeee tee 
COMMON XLO(2) sXHI(2) sERRsRFsFIN( 100) sTIN(100) sw(550) 
COMMON R{14OO) sTIME (40950) 9 (40950) 
DATA WEAR/'SS's'FC fe 
DATA Sekt a os 
lL FORMAT(8F10.3) 
2 FORMAT(4012) 
3 FORMAT (5F100e3 56F 501) 
GO000 FORMAT ('1L's///s52Xs'SUMMARY OF INPUT DATA's/) 
6001 FORMATI5Xs!' INTERVAL AF c S SPEED!» 
aie WEIGHT CONSTRAINTS CT Sie err" 
2%' COST,, TRIP RATE CHANGE TIME WEAR's/937Ks 
S3"RESPONSE RESPONSE WMIN WMAX RMIN RMAX (TROHES.))" 
a4o! (HOURS) (HOURS) (HOURS) TYPE ey 


clos MORMA TN SS Xe KRG CLUE  DeGeR CE Sor 7 99 

6002 FORMAT( 2X sF6e08! “=! pF SeOSlLXsFOeZslXsF6e4sl1XsF5e0253Xs 
LAS SSX aF4e295XaFBZeOsBXoFZeO0s2XaF oes 2XoF4e0s3xXoF Ge Bs 
LoXeFGeB92X9Poe2s!'/LOOOFT ! s4XsF5e296X9A28///) 

6020 FORMAT(27Xs'BIT NOs INTERVAL WEIGHT SPEED f 
Te tROTs TIME TRIP FIME EQUIVe (COST WEARS s/ S87R 5 
2 RVG HOURS, tT 9B? a7) 

6006 FORMAT (29Xsl2Z292XaFlelotm=—! sF7ale2krsFae0s4X oF 4e093Xs 
LP7e2svGX oF Te2sSKsFlezeSXelle2Xelia//) 


6007 FORMAT(39Xs'FORMATION ROTATING TIME = '»sF7e2, 
1' RIG HOURS's//» 
1 39Xs"'FORMATION TRIP TIME = ‘ek 7eZs 
1! RIG HOURS! »//» 
i 39Xs'FORMATION RIG TIME a ees 
1! RIG HOURS's//9 
2 39Xs'FORMATION COST = 'sFT7e2ys 
1’ RIG HOURS'»s//>» 
2 39Xs'FORMATION BIT COUNT = 's3Mal2e//) 


6008 FORMAT(39Xs' CUMULATIVE ROTATING TIME = 'tsFT7e2s 
1’ RIG HOURS's//s 
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a MAIN eee (CONT'D) 
Fy 39Xs' CUMULATIVE TRIP TIME = Ver ire s 
L©° RIG HOURS! s//s 
2 39Xs'CUMULATIVE RIG TIME = '5Flae2s 
TT * RIG HOURS!" y/ 7s 
2 39Xs'! CUMULATIVE COST = Rey wes 
1" RIG HOURS! s//5 : 
2 S9OXs! CUMULATIVE BIT COUNT = Ns SKa 1259/7) 


NEW INPUT VARIABLES 


OF INTERPOLATION TO BE USED WHE! 
Or 2 .2N feelin VSED WHEN FINDING. COST 


NOl = THE NU R OF CIFFEREN FOOTAGES FOR WHICH A 
ROTAT I TIME Is CALCULATED. 
NOMAX = THE NUMBER OF STATES AT EACH STAGE. 
NEWSET = NON ZERO VALUE RESULTS IN READING ANOTHER 
DATA SET 


ICBIT=0 

CTROT=O0 6 

CTRIP=0.6 

CRIG=06 

CEQ=068 

READ(592) IWlsIR1sIDEG1sIDEG2sNO1 »sNOMAX 


STEQ=0 o/ 

READ(Se1) MHOLEsBITCSTsRIGCSTsTRIPT 
READ(5s3) SDEPTHsTDEPTHsAF sCF sDFS» 
LIXLO( J) aJSlo2) s( AHI (0) evel s2) 


a ae eae 
FO RM »sFORPIK » 


INTERNAL VARIABLES 


PINC ) = FOOTAGE VALUES FOR WHICH ROTATING TIMES ARE 
FOUND 

TING 9) = ROTATING TIMES ASSOCIATED WITH FINI } 

Pe taie 9) OF IDERP TRS At SPAGE wu 

TIMEt Ss) = VALUES OF OPTIMUM COSTS ASSOCIATED WITH 

ae ae) 

WINT( UO) = OPTIMUM WEIGHT FOR THE STAGE J 

RINT(w) = GPTIMUM SPEED FOR THE STAGE ug 

TROT( J) = ROTATING TIME FOR STAGE ov 

TTRIP(J) = TRIP TIME ASSOCIATED WITH STAGE J 

TEQ(s)- = COST ASSOCIATED WITH STAGE J 

FOUT(U) = DEPTH AT TOP OF STAGE J 


IP=IFIX(P/069)+1 
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CG“ 
e MAIN eae (CONT TD) 


BITHRSBITCST/RIGCST 
TRIPHR=TRIOT*1OOO. 
LF UFORMsGTs0.5) ITUK=2 
PERFORM sk Tw 0051 I JK= 
WRANGE= ANTI I=XLOC1) 
RRANGE=XHI(2)=xXLO(2) 

I 


W=IW1L*I1F 1 X(WRANGE +005) +1 
IR=IR1*I FIXURRANGE+005) +1 
TFCUIWsEQsel) Iwe2 
[IF(IRsEQel) IR=2 
DO 16 J=lsIw 


16 W(J)=XLO ) 
DO 1? J=lsIR 
17 R(J)=XLO(2)+RRANGE*FLOAT (J=1) /FLOAT(IR=1) 
RANGE= TDOERTH =SD SEPTH 
GALL FESTIU(XLO( 1) »XLO(2) sFOOTIN) 
FOOTINSFOOTIN#®4999 
IF CRANGEsGTeFOOTIN) FMAX=FOOTIN 
IF CRANGEeLEsFOOTIN) FMAX | 
CALL TESTIMXAHD Cl leXHI (2) sFOOTING 
IFS FOOTINeGTaFMAX) FOOTINEF 
FMINSFOOTIN/2.6 
RF=FMAX“FMIN 
DO 10° J=lsNOl 
9 PINGS) SFMINtRFHPLOAT ( J=1)/7FLOAT Ui ei) 
O CONTINUE 
DO ll J=leNOl 
CAL SEARCH(CFING J) sWRESTsRBEST sTINis) sDOUT sBOUT) 
11 CONTINUE 
MINBITS=IFIX(RANGE/SFMAX) +1 
TP C(MINBITseGis40) GO TO 1lil 
FOPTHlseE+20 
DO 50 J=1940 
FSTART=TIDEPTH=FLOAT (WI) *(FIN(NOL)=3.4) 
PPL FSTARTsLisaSDEPTH) #*START=SDErPTF 
PENDS (DEPTHRFLOATIO) PUR IN ( 2) Foe] 
RANGEL=FEND=FSTART 
DO 70 K=lsNOMAX 
TO FlosK)=FSTART+RANGEL*FLOAT (K=1) ZFLOAT (NOMAX=1 ) 
MeiwecQed) GO TO 71 
CALL FIND( Jel) 
TFC(MINBITeGTev) GO TO 72 
IF( TIME (J91)eGEeTOPT) GO TO L000 
TOPT= Laan a 
72 DO 73 Ke2sNOMAX 
jae Sacle kh ae 
GO TO 50 
71 DO 74 K=1lsNOMAX 
DELF=TDEPTH=Filsk) 
Tee Time ULV) = XNPOLZ (DELFI+B I THRSREP PING tie ee DEP TH) 
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fk MAIN eee (CONTID) 


LETRIPT 2 
SE UMINBITTeEQel} TOPT=TIME (2,2) 
20 CONTINUE 
Cogn 27] 
1000 LBIT=J=-1 
FNOW=SDEPTH 
FOUT(1)=SDEPTH 
ey 40.00 “K=1 shat 
KKSLEIT+1L=<K 
LEVKKReEGQal) GO TO 3000 
CALL FINDI(LBEITs KK sFNOWsFINT(K) } 
GO TO 5000 
3000 FINT(K)=TDEPTH=FNOW 
BO00 CALL SEARCH(FINT(K) sWINT(K) sRINT(K) sTROTI(K) sDINT(K) 
*eBINT(K)) 
ID(KISIFIX(8e* DINT(K) +065) 
PR(K) = IP IX(8e* BINT K +4005) 
TTRIP(K) S(FNOW+FINT(K) /20)*TRIPT+REPTIM 
TEQ(KISTTRIP(K)4+TROT(K)+RITHR 
STROT=STROT#TROT(K) 
STRIP=STRIP+TTRIP(K) 
STEQ=STEQ+TEO(K) 


POUT CKTLISFOUTIK I +RINT (CK) 
4000 CONTINUE 
CUROTECTROM+S TROT 
CTRIPHCTRIP+STRIP 
CEQ=CEQ+STEQ 
TRIPT=TRIPT#1OO0« 
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(Z)sXKHT (2) sHOLESBITHRSTRIPTSs REPT IM 


6005 WRITE(6536006) UsFOUT(U) sFOUT( U4l) sWINT(u) sRINT (CU) 
SeTROT( I) eo TTRIP() 
ete a) RO Le LB) 
SRIG=STROT+STRIP 
WRITE(6 96007) STROTsSTRIPsSRIGsSTEQsLBIT 
CGH IT=H=TCBIT+L Bit 
CRIGHECTROTFCTRIP 
WRITE (626008) CTROTsCTRIPsCRIGsCEQs ICBIT 
haan ye is Re 
Pont WRITE (6s, 112) 
1112 FORMAT (! ERROR =RROR ERROR Vs/y! MAX WO 
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FUNCTION XNPOL1 
FUNCTION XNPOLL(DELF) 


THIS FUNCTION SUBPROGRAM INTERPOLATES THE APPROPRIATE 
ROTATING TIME FOR A GIVEN FOOTAGE USING LAGRANGIA 
INTERPOLA TIONe 


OMMON AF sCFsDFS sFORMsFORMK sBITHR 

COMMON TDEPTs SDEPT MesHOLE sRIGCSTsBITCSTsTRIPT sREPTIM 
COMMON NOlsNOMAX sIWsIRsPsIPsIDEG1sIDEG 
COMMON XLO(2) sXHT(2) sERRORFE sFIN(100) sTIN(100) swt550) 
COMMON R¢1400) sTIME(40550) sF (40550) 
TP CCDEL Pel TePINGD ))eORe (DELFseGTeFIN(INOL))) GO TO 3 
K=IDEGL+1 


De OP areLeeP LN (Kiel GO TO: '2 
K=ak+] 

TFP(KeGTeNOl) GO TO 3 

Se VG) i 

Jl=K=-IDEG1 

J2=K 

XNPO =O-6 

Oe i a ie 

XL=Ale 


COV 2 Kelas? 

De tank Gak) GOTO B2 
Aie=hL eC BELPHFING KY US CE EN Col OAKS 
vedi 

XNPOLL=XNPOLI14+TIN(U)*XL 


CONTINUE 
RETURN 
ANPOLL=leE+LO 

RETUR N| 

ENT 
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FUNCTION XMPOL2 
FUNCTION XNPOL2(FPREVsNSTAGE) 


THE FUNCTION SUBPROGRAM INTERPOLATES THE COST 
ASSOCIATED WITH THE GIVEN VALUE OF THE DEPTH. 


COMMON AF » CF sDFS sFORM»FORMK sBITHR 
COMMON TDEPTHs SDEPTH sHOLE sRIGCSTsBITCSTsTRIPTsREPTIM 

COMMON NOLSsNOMAX sIWsIRsPolPsIDEG1sIDEG2 

COMMON XLO(2) 9XH1I(2) s9ERR»RFoFIN(1O0) sTIN( 100) 9W(550) 

COMMON R(1400) » TIME (40950) sF(40550) 

JeNSTAGE= 

LECFPREVelLTsePlssii) GO To 3 

K®IDEG2+1 

Pret PPREVALESFUUsK)) GO TO Z 

K=K+1] 

TF (KseGT eNOMAX) GO 7O°3 

SO. TO 1 

w#iL=K=IDEG2 

J2=K 

XNPOL2=06 

PO 1g Kevlaw2 

XL=le 

OO" a2 abe is v2 

lPCLeEQ@eK) GO TO 12 

AueRee tC PPREVSP (Jeb) ) SCR Le Ke R Coy Ed) 

CONTINUE. 

ANPOL2Z=TIME(UsK) *XL+#XNPOL2 

CONTINUE 

RETURN 

ANPOLZ=leE+10 

RETURN 


END 
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SUBROUTINE FIND 
SUBROUTINE FIND(NSTAGEsNO) 


THIS SUBROUTINE CALCULATES THE OPTIMUM COST FOR A 
PARTICULAR VALUE OF THE STATE VARTABLE (F(NSTAGEsNO))« 


COMMON AF sCF sDFSsFORMsFORMK sBITHR 
COMMON TDEPTHs SDEPTH sHOLE »RIGCST+BITCST »TRIPT sREPTIN 
COMMON NOLSNOMAX sIWsIRsPsIPsIDEG1s1DEG2 
COMMON ALO(2) sXH1(2) sERRe RFs FIN( 100) sTIN( 100) sw(550) 
COMMON R(1400)sTIME(40550) sF (40550 
SS —3.6 
TIMECNSTAGEs NO) =1 eE+10 
DELF=(TDEPTH=F(NSTAGEsNO) ) /FLOAT(NSTAGE)=1e 
DELF1l=DELF 
999 R=XNPOLL(DELF)+(F(NSTAGEsNO)+DELF/2¢)*TRIPT+BITHR+ 
IXNPOL2( (F(NSTAGE »NO)+DELF) sNSTAGE)+REPTIM 
IFC TIMEXCNSTAGEsNO) «GT ex) Gov hO:- 7 
oo =S5/3¢ 
TF (CABS(SS) «LTeERR) RETURN 
DELEaDELF L+5S 
Bev TO 99 
L TIME NSTAGE »s NO) =X 
RPeiE d= DeELr 
SS=SStZ » 
DELFa=DELFI+tSS 
GO TO 999 
END 
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SUBROUTINE FIND1 


SUBROUTINE FINDI(LBITsNSTAGE sFNOWsDELF1) 


C THE SUBROUTINE CALCULATES THE OPTIMUM DISTANCE TO BE 
c DRILLED AT A PARTICULAR STAGE IN THE PROCESS.« 
C 


COMMON AF sCFsDFSsFORMsFORMKsBITHR 
COMMON TOEPTHs SDEPTH sHOLEsRIGCSTsBITCSTsTRIPT sREPTIM 
COMMON NO1LsNOMAX sIWsIRoPsIPsIDEG1sIDEG2 
COMMON XLO(2) »XH1I(2) s9ERR»oRFsFIN(100) sTIN(100) sw(550) 
COMMON R(1400) s TIME (40950) 9F (40550) 
SS=~36 
TIM=1lsE+10 
DELF=(TDEPTH=SDEPTH)/FLOAT(LBIT) = he 
DECFLEDELF 
999 X=XNPOL1(DELF)+( FNOW+DELF/26)*TRIPT+BITHR+XNPOL2 | 
lL ((FNOW+DELF) sNSTAGE)+REPTIM 
PFATIMeGt er) “GO°( TO] 
S5=—"SS/3's 


DELF=DELFI+SS 
GO TO 999 
1 TIM=xX 
SS=SS%*2 0 
b. PELPISOELF 
VELreDELrFi+ss 
SOLO see 
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SUBROUTINE SEARC 


SUBROUTINE SEARCH(FOOTIN»s WREST sRBEST»TBEST> 
LOBEST»BBEST) 
COMMON AF sCFsDFSsFORMsFORMKsBITHR 
COMMON TDEPTHsSDEPTHsHOLE sRIGCSTsBITCSTsTRIPTsREPTIM 
COMMON NOLsNOMAXsIWsIRsPsIPsIDEG1sIDEG2 
COMMON XLO(2) »XH1(2) sERRe RF oFIN(100)-9TIN( 100) sw(550) 
OMMON R(1400)sTIME (40550) sF (40950) 
TBEST=1.E+10 
ITEMP=IW 
DO 13 K=lsIR 
IND=0 
TJ=ITEMP 
Dr OPTEMP eek O a0) GOTO. 14 
Be oie wweleig 
LFtTNDwEG«l) GO TG 13 
J2TUtleud 
CALE TEST (WK gar LK 


)sFOOTINS FXsDXs 8X) 
LF(FX«eGTeleEt+O8) GO T hae 


Ral a 
IND=1 
PPMEARCESTBESS} 160 TO.12 
TBEST=FX 


RBEST=R(K 
WREST=W(u) 
DBEST=Dx 
BBEST=BxX 
Go 7O 12 
ITEMP= 2-1 
CONTINUE 
CONTINUE 
RETURN 


END 
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SUBROUTINE TESTI 


SUBROUTINE TESTI(WEIGHT sREVSsFOOTIN) 
COMMON AF sCF sDFS sFORMsFORMK sBITHR 
COMMON TOEPTHs SDEPTHsHOLEsRIGCSTsBITCSTsTRIPT»REPTIM 
COMMON NOL SNOMAX sIWeIResPsIPsIDEG1sIDEG2 
COMMON XLO(2) sXH1 (2) tERReRFsFIN( 100) sTIN( 100) swt550) 
COMMON R(1400) sTIME(40550) sF (40950) 
DULLA(D)= 1.0+(660+.928125%D)#D 
DULLAT(D)=(1604+(30+0299375#D) *D) *D 
DULLAH(D)=(1.0+(640+.928125%*D)*D)**0.5 
SIM1(D0)=(16019779+(14292304+( 61133279*D=.4127675) 
1*#D)*D)*D 
IP=ITFIX(P/069) +] 
TPC AF el Tele E=C5) AF=le 
A=EXP(“100e/REVS*¥#2) 
LFNPORM=005) lal #2 
ARR=A*REVS##O 0 7540 e5*%REVS*(1Le—A) 
SOO" Tad 
ARR=A*REVS#*#0 642840 e2*#REVS#(1le—A) 
WBAR=7«e688*WEIGHT/HOLE 
EMWATE=1359e¢1-714e19%ALOG1O(WBAR) 
ROTISREVS# (Lette 00004348 #REVSH#2) 
KLOAD=Z21338e87/(1LeteOB1187*WBRAR) ##362258 
FI=AF#EMWATE/ROT I 
F2=FI1*CF* (WBAR*¥*FORMK ) *ARR 
FS=FL*REVS/DFS/XLOAD 
GO TO (20820930) s1P 
BOUT=F3*DULLAT (ls) 
Pe sBOUlebEsl sO) GO TO 20 
K=1ls 
DMAX=X=(F3*DULLAI(X)=1e)/(F3*DULLA(X)) 
IF (ABS( (DMAX=X) /DMAX) eLEsleE=05) GO TO 60 
X=DMAX 
GO TO 6 
DMAX=1 + 


GO.TO (4840330) s5I1P 
FOOTINZFZ*S IML (DMAX) 


RETURN 
FOOTIN=F2Z*DMAX 
RETURN 
FOOTIN=F2/F3 
RETURN 

END 
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SUBROUTINE TEST 


SUBROUTINE TEST( WEIGHT sREVSsFOOTIN »TEFF sDOUT »80UT) 
COMMON AFsCFsDFSsFORMsFORMK sBITHR 

COMMON TDEPTHs SDEPTHsHOLEsRIGCSTsBITCSTosTR RIPTsREPTIM 
COMMON NOLsNOMAX sIWsIRePsIPsIDEG1 sI DEG2 
COMMON XLO(2) XH 62) vERR ORE sFING OO} eh int 100} ew { S50) 
COMMON R(1400)sTIME (40950) sF (40350) 

DULLA(D)= 140+(600+.928125%D)#*D 
DULLAT(D)=(1e04(30+6309375%D)¥D)¥*D 
DULLAH(D) =(14041640+. 928125*D)#D) ##0.5 
SEML(D)=(1e019779+( 102925044 ( «1153279%D=.41276 75) 
#D14D1*D 
[P= FIX(P/069)+] 
ASEXP(-100e/R REVS*#2) 
IFC FORM=005) lsls2 

ARRE=AXREVS#*#O 6 7540 e5*#REVS#(Le=A) 
GO-TO 3 

ARREAXREVS##O 042840 e2*REVS#(1le—A) 
WBAR=7e658*WEIGHT/HOLE 
EMWATE=1359e1=-714.19*ALOG1LO(WBAR) 
ROTI=REVS#( Leto 00004324 8#REVSHHD) 
XLOAD=21338e87/(1lseteOZ1LLISTHWBRAR) ¥¥#3.2258 
GO TO (20320530) sIP 
PLE=AF*EMWATE/ROTI 
F2=F1l*CF* ( WBARX*FORMK ) #ARR 
F3=FI*REVS/DFS/XLOAD 


RO 


GO TO (405350330) sIP 

C=FOOTINAR2 

K=e5 

DOUTH=XKSiSIML(X)}—-C) /DLLLAH CX) 

1G ABS I (DC UT—xX) SDOUT sai Tedeb=057 GOo 1G S 
X=DOUT 

GO TO 4 


Let DOUTsGTele) GO TO 10 
DDIFFSDULLAI (DOUT) 
BOUT=F3*DDIFF 

IP (BOUT «GTels) GO TO 10 
TEPFEPLeYeDDIFF 

GO JO il 

TEFReL eb? 10 

RETUR 

DOUFSFOOTIN/Fe2 

CEtpOVT eGleiseO)? GO TO 20 
DDIFF=DULLAI(DOUT) 
BOUT=F3*DDIFF 

TE ( BOUT seGTeled) wo TO oO 
TEFF=FL*DDIFF 

RETURN 

FI=DFS*XLOAD/REVS 
BOUT=FOOTIN/( F1*CF*(WBAR**FORMK ) ARR) 
Pe CROUTeGTele0) GO TO 10 
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RESULTS FOR BASIC ‘INTERVAL! METHOD WHERE 
DRILL EQUAL PARTS OF THE INTERVAL. THERE 
GROUP FOR EVERY CHANGE IN THE CONSTANTS« 


GROUP NOs 1 


BIT NO INTERVAL WEIGHT SPEED 
Ag 14008 = 1907. 596 LOA 
2 1907e cae 24136 59 « 1336 
2 24136 = 29206 598 1336 
GROUP NOs 2 
BIT NOs INFERVAL WEIGHT SPEED 
d 292068 =~ 33454 596 1286 
2 3345 ath ay: O-* 59 « 128.6 
GROPP Os 3 
BIT NOs INTERVAL WEIGH SPEED 
a 3770+ = 40956 59» LS Os 
Zé 4095@ = 44206 596 ae 
2 44206 te Lf Le Sie 59s L306 
& “ 474568 °—= 5070» 59 130.6 
GROUP NOe 4 
BIT’ NOe INTERVAL WEIGHT SPEED 
Ih 5070 « baal 54956 59s 128. 
A 54956 = 59206 59's LZee 
3 59206 oe 63456 BOs Leb 
4 63456 ae 67706 59 « 128.6 
GROUP NOs 5 
BIT NOe INTERVAL WEIGHT SPEED 
1 6/106 = 68106 506 606 
| GROUP NOs 6 
7 
7 BIT NOs INTERVAL WETLGHT SPEED 
1 68106 Roms, (Lada 596 L086 
2 T1455 = T4805 59 108-6 
3 74808 = T7ElS5s 59 108. 
& 78156 = 8150 « 59-6 1086 
e) E1504 =~ 84855 596 L086 
6 B4B5<« = 88208 55s 108s 
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GROUP NOe 7 


BIT NOs INTERVAL WEIGHT SPEED 
1 Ba20 6 mi 88606. 50 606 


GROUP NOe 8 


BIT NOs INTERVAL WE LGHT SPEED 

1 oS608e = 90205 60+ LOOe 
GROUP NOs. 9 

BIT NOe INTERVAL WEIGHT SPEED 

1 90206 = I060 DOs 606 
GROUP NOs 10 

BIT NOs INTERVAL WEIGHT SPEED 

i 90606 = 92705 60 90.6 


GROUP NOs lil 


. BIT NOe INTERVAL WETGHT -“SPEED 
Hk 92/09 —~ LO0SOs sea S9e 


GROUP NOe le 


Bil NOt INTERVAL WELGRT «SPEED 
1 LOO8Os. — 104154 456 65-6 
2 10415 end 107506 45» 656 

GROUP VO Wes) 

BIT NOs INTERVAL WEIGHT SPEED 
i ONC ene Loe LO8636 456 Gk e 
cs LO975 es = TL0CGe 43 44. @ 
& 110866 =~ RE PATON 454 Abe 
5 1120008 =~ 1131308 456 “deg 

GROUP NOs 14 
BIT NOe INTERVAL WEIGHT SPEED 


i PLS ae <= £1680. 456 10° 
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BIT NOs 


ho 


CUMULATIVE 
CUMULATIVE 
CUMULATIVE 
CUMULATIVE 
CUMULATIVE 


GROUP NOs 15 


INTERVAL WEIGHT 
116800 = 11882. 45.6 
118820" = 12084. 45.6 

GROUP NOs 16 

INTERVAL WETGHT 
120848 = 123636 436 
123636 oJ 12641. 456 
12641. kn L2920« 456 


GROUP NOs 17 
INTERVAL 


1292065 = 


120686 som 


13068. 43 
1321600 74a 


He) 


Y el m1 Ais ¢ 
GROUP NOs 18 


INTERVAL WEIGHT 
13216. lsd LoS Go 456 
133638 = 13510 456 
ROTATING TIME = 806666 
TRIP TIME = 2862674 
RiG TIME = 1091640 
COS s = 1503444 
BIT. COUNT = 42 


SPEED 


SrPLep 
686 
636 
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RESULTS FOR THE EQUAL INCREMENT VARIATION 


PINTERVAL! 


Wo Mh 


to PO 


BIT, NOs 


& Ko 


BIT NO« 


Ow fF Wh 


OPTIMIZATION PROCEDURE. 


GROUP NOe 1] 


INTERVAL WETGHT 


14008 = 190C. 60. 

19C0Os = 24060 59 6 

24066 = 29200 59 6 
GROUP NOs 2 
INTERVAL WE LGHT 


29206 bs) 


33454 596 
33450 = 9 


Ot taglines acs 


CEAIID AIC A 
‘ tf 4 
GROUP iwUe 3 


INTERVAL WEIGHT 
3/106 os GO94.4 59s 
40946 =— 4alGe 59 ~@ 
44180 = 4744. 59. 
Le eae = 50706 59s 

GROUP NOe 4 

INTERVAL WEIGHT 

5070.« — 5484 6 596 

4845 —= 5905. 59 «6 

5905e¢ = 6334s 29° 

33 | tae 67706 59s 
GROUP NOs 5 

INTERVAL WEIGHT 
67702 = 68106 30s 


GROUP NOe 6 


INTERVAL WETGHT 
68lQ. = 7138s 60% 
7T13d50 = 746958 596 
T4690 —~ 78038 596 
T8036 = B1396 59. 
813958 ~= 84785 58s 
8478 « Me! 85206 586 
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1306 
129. 
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110.6 
LU96 
1086 
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BIT NO« 


FORMATION 


BI NOs 


wm Ew he 


GROUP 


NO 


INTERVAL 


GROUP 


NOe 


INTERVAL 


$2606 


GROUP 


9020 


NOs 


INTERVAL 


GROUP 


9060-6 


NO» 


INTERVAL 


INTERVAL 


100806 
104136 


GROUP 


INTERVAL 


Vv 
RY 
~~ 


10 
108626 
109756 
L1IO87« 
L1IZ200~8 


INTERVAL © 


Lisle 


116806 


P 
10 


WE TGRT 


WE TIGHT 
45 
4 Oe 
456 
45 
45» 

14 

WETGHT 


45 » 
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GROUP NOe 15 
BIT NOs INTERVAL WEIGHT SPEED 
Va 11882. = 120846 45 6 106 
GROUP NOe 16 
bit “NOs INTERVAL WEIGHT SPEED 
ty 120845 = 1236256 45. 68s 
2 123626 2 L2Z641+6 45 6 686 
2 12641le = 129206 45. 68 
GROUP NOs 17 
BIT NOs INTERVAL VETGHT SPEED 
L 129206 = 13067.» 46% 60 
fi: L3206/. aa 132164 Dy 5 Die 
GROUP NOs 18 
BIT NOe INTERVAL WEIGHT SPEED 
= i 32 LG. ec 13363 45 6 106 
133636 = LS pO 4 De 106 
CUMULATIVE ROTATING TIME = 808637 RIG HOURS 
CUMULATIVE TRIP TIME = 282669 RIG HOURS 
CUMULATIVE RIG TIME = 1091206 RIG HOURS 
CUMULATIVE COST = 13034610 RIG» HOURS 
CUMULATIVE Bia COUNT = 42 
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RESULTS FOR THE DYNAMIC PROGRAMMING VARTATION OF 
PINTERVAL! 
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BIT NOs 
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BIT NOe 
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OPTIMIZATION 


GROUP NOs 
INTERVAL 
1400. = 1894.6 
18945 = 2395. 
2a956 be cx 29206 

ROUP NOe 


INTERVAL 


) = 33426 
c S77 O% 


1 ee am 4087s 
408/58 ™ 4410 
44106 = 4731 
47316 = 5070s 


INTERVAL 


50706 oe 


54/736 
54730 = 59136 
59136 = 63376 
633/656 = Gf roe 
GROUP NOs 
INTERVAL 
67708 —@ 6810.4 
GROUP NOe 


INTER VAL 
681046 = 71426 
TL4&a2e ie 74756 
T475¢ = 78096 
78098 = 8143s 
81430 —- 6480. 
8480e8 = BB2U.6 


PROCEDURE s 


1 

wE TIGHT SPEED 
606 “laa% 
606 L326 
596 1296 
2 

WELGHT SPEED 
59s L296 
296 L27s« 


WEIGHT SPEED 


DI e Lee 
60. 128. 
29s LS ls 


y 


WETGHT  - SPEED 


60. L131. 
290 124.6 
aoe 128 
59 6 1266 

D 
WEIGHT SPEEDO 
5.06 606 

6 


WEIGHT SPEED 


D9 « L109 
53 6 LOG. 
59 1096 
606 1066 
596 1086 
28» 109-6 


G=34 


-_ *) ie 
“werd 


BIT NOe 


4 

e 
RIT NOs 

1 

2 

3 

4 

S 
SIT NOe 


GROUP NO« 7 
INTERVAL WEIGHT 
86c040 = 886054 50.6 
GROUP NOs 8B 
INTERVAL WE TGHT 
E8606 = 90206 606 
GROUP NOs 2 
INTERVAL WEIGHT 
902056 =~ 90604 50s 
GROUP NOs 109 
INTERVAL EIGHT 
90606 ~ 92706 60 
GROUP NOe Tl 
INTERVAL WEIGHT 


1B 0e 52s 


INTERVAL 


30 10405 
10405» = 10750 


INTERVAL 


LOTSOn, = LOSes 
108638 = 10976 
10976+ = 11068¢ 
LLOBBs. = 1L2o00 
LL ZOU ee Se 


GROUP NOe 


INTERVAL 


peas = Lee e 


12 
WEIGHT 
» 45.6 
« 45e 
te 
WEIGHT 
eo 6456 
0 45e 
© 456 
» 456 
e 456 
14 
WEIGHT 
0 6456 


SPEED 
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GROUP NOe 15 


BIT NOe INTERVAL WETGHT SPEED 


2 l1l88le = 12084. 45 5 waer 


GROUP NOs 16 


BIT NOs INTER VAL WEIGHT SPEED 
1 2084s = 123625 458 586 
2 1236256 = 12640. 45+ 686 
3 1264056 =~ 129205 45.6 67s 


GROUP NOe 17 


BIT NO’ INTERVAL WEIGHT SPEED 
1 1292048 = 13068. 456 59-6 
2 130686 ™ 132168 456 59 
GROUP NOs 18 
BIT NOs INTERVAL WELORI SPEED 
1 Ll3s2l6ée = 1336246 45% T06¢ 
=- 2 LSe62 0) SLs SiO 45 6 106 


808e«80 RIG HOURS 

282067 RIG HOURS 

1091648 RIG HOURS 

1303651 RIG HOURS 
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CUMULAT- EVE ROTATING TIME 
CUMULATIVE TRIP TIME 
CUMULATIVE RIG TIME 
CUMULATIVE COST 
CUMULATIVE BIT COUNT 
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Program and Results of "Multi-interval" Optimization 
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C MAIN 
MAIN 
THE PROGRAM IS DESIGNED TO OPTIMIZE THE DRILLING OF 
AD abate USING THE DYNAMIC PROGRAMMING METHOD OF SECTION 
Tee 
DIMENSION XLO(2) sXH7T(2) sRESP(2) sMAPING(10) sWEAR(3) 
COMMON AF(150) sCF(150) sBFSt150) 
COMMON FORM(150) sFORMK (150) sDEPTH(151) 
COMMON XLIMLO(15092) sXLIMHI (15092) sNOS1(151) 
GOMMON WIN(150)sRIN(150) sTOTALT( ine et bpracbetat ng 
COMMON XNDIC(C1ISs15)swliSs15)sR015815)sDE(15) »BB(15) 
COMMON ERR» ATi CLeAINC2yREPTINyCORR(2) TH IPT »P 
COMMON HOLE sBITCSTsRIGCSTsFOOTIN IT sBOUT sLAGsDELT 
COMMON ITPsXBAR 
hg TOPT3s1OPT4sIC 2TS5s10PTS 
DATA RESP /*SORT! s “HARD! 7 
DATA fEURA' 854s ee Ug tee hs 
zh essai ety 
2 RMAT(4012 
3 FORVAT(SF10«3 96°32) 
6004 RMAT{'LT'’s4GXe"P O L IT -C ¥¥as/7) 
6005 en aan Matar INTERVA WEIGHT SPEED ROTe! 
Bee a Me PRP TAME OSs WEAR! s/s72X%¢s' (HOURS) T 8B! 
28e7,) 


6006 FORMAT(LEX aI 293X sF6e09'=! sFGn0slXaF4e0e3X oF4e053Xs 
SP ieexegrnslise2rxal less) 

600] FORMAT ¢25X>s'SUMMARY'! sl Xa FGeOs'=!' sFG6e0e1S5XsF le2a4Xs 
LF feZ2rZz2XsFfa2s2XeTls2XsTla/ 7} 

2002 FORMAT(S7Xs "BIT CHANGED AT ‘sFGe0s" FEET "e977 ) 

60G0 FORMAT('1's55X5' SUMMARY OF INPUT DATA'’»s///) 

6008 FORMAT(3Xs!' INTERVAL AF Cr 5 SPEED! 
L4Xs'* WEIGHT CONSTRAINTS SII CSLZE LBA -COSt 
Let TRIP RATE! 
2es2Xs'CHANGE TIME YG XG D WEAR's7935Xe'RESPONSE! o 1X9 


3'RESPONSE WMIN WMAX RMIN RMAX (INCHES) (HOURS)! 
493Xs! (HOURS) (HOURS) TN Oe hig ay 


6009 FORMAT(1XsF6e Ost='sFGeOslXsF6e29lXsFGae4alxoF5e2s 
LsBXsAGsSXoF4e2o3XoFoeOelxeFGe0olXsF4e0slXsF4e093Xs 
2F60394K oF Ge da2XaF4e2et/lLOOOFT 's4X oF Sets4XoT2elxXel2s 
31XsIl1l93XsA2s/) 


6003 FORMAT(28Xs'CUMULATIVE ROTATING TIME = 's FTe2s 
1' RIG HOURS's//s 
; 28Xs! CUMULATIVE TRIP TIME = ly Feds 
1! RIG HOURS! s//s 
1 28X%s'CUMULATIVE RIG TIME ee ee 
1! RIG HOURS! »s//s 
2 28Xs'CUMULATIVE COST SIGE PREZ 
1' RIG HOURS'»s//s 
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RK MAIN eee (CONT!D) 

4 BX» CUMULATIVE BIT COUNT Sty boa he) 
6010 FORMAT(28X9'6 ROUP ROTATING TIME = Se Fheg s 

1! RIG HOURS's//s 

LA,* 28Xs'GROUP TRIP TIME 2's Fie2s 

1’ RIG HOURS'»s//s 

2 28X%'GROUP RIG TIME ed aie ees 

LRG OHOUREIS 77% 

3 28Xs'GROUP COST = Np Pras 

L* RIG HOURS! 3/75 

et 28X»'!GROUP Bit ae iT ee Wt a 
BOZO FORMAT (814 9 BX g tH 55H 4696 Ee EE TE ESE SE IE FE IE HE HE SE FEE FE 9 HE 96 9b HEE dE 


hai « *Ty1lOXs WARD 1 IG! 11x! HV gg / aig! Bl y2Q2X gtk! yy 
24X9'* DEPTHS ARE NOT IN SEQUENCE *'5/s4Xs 
3 et ee ae CH HIE Te HEH FE TE TE HE FE GEE 46 EET) 

GO2Z21 FORMAT (81! 3X 9 EHH EHH HEHEHE HNEKRHHHHKREHHET g/g 
LAPP VELSLESKLERROR' JIN 5 ay fay Onl em ras Meee 
» i r ‘ cr FN Me *< 
2'* NO OF INTERVALS EXCEEDS 150 *'9/54X51%!529K5 tx! 
Bs/o4Xs!#'y1lOXe! RUN Pete liey a ops 
4 V3 FE FE SE IE SE SEE SE 9b FE HE SE S638 SE SE TE GE GE SE ES HICH EH NY 


IPUT IS SIMILAR IN MANY “RESPECTS TO. THE- PREV Rous 
GRAMSs. AND. SO ONLY. THE DIFFERENCES WILL BE NOTED 


INPUT DATA 


LAG = DEGREE OF INTERPOLATION IN SUBROUTINE INPOL 
I IF -NON ZERO THE RESULTS OF THE 'BACKWARDS! 
CALCULATIONS ARE PRINTED OUTs 
10PT2 + IF NON ZERO A MAP OF THE SEARCH AREA IS GIVEN 
TOPT3 = FIRST DULLNESS ELEMENT FOR WHICH A MAP 1S 
GIVEN 
IOPT4 = NUMBER -OF DULLNESS ELEMENTS SEPARATING EACH 
MAP EVALUATIONs 
IOPTS - NUMBER OF THE FIRST BEARING WEAR ELEMENT FOR 
WHICH A MAP If PRODUCED. 
1OPT6 — NUMBER OF BEARING WEAR ELEMENTS SEPARATING 
EACH MAP EVALUATIONs 
MAPNO = THE NUMBER OF DISTINCT STAGES FOR WHICH A MAP 
1S PRODUCED. 
MAPING = VECTOR WHOSE ELEMENTS GIVE THE STAGES AT 
WHICH A MAP WILL BE PRODUCED. 


©) 
Do 

4 

[= 
f 


INO = THE NUMBER OF STAGES FOR THE FORMATION 

INOS1 — THE NUMBER OF DULLNESS VALUES AT EACH STAGEs 

INOS2 - THE NUMBER OF BEARINGS VALUES AT EACH STAGE. 

IND1 - IF ZERO THERE IS ANOTHER FORMATION. IN THE GROUP 

IND2 - IF NON ZERO THE STATE VARIABLE MATRIX VARIES« 

IND3 < IF NON ZERO THE CUMULATIVE COST FIGURES ARE NOT 
PRINTED 
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C MAIN ees (CONT'D) 
c 

CEQ=06 

CROT=05 

CTRIP=QO6 

KKOUNT=0O 


9999 READ(592) LAG» IOPT1»1OPT2»10PT3s10PT4s1OPT5s10PT6> 
IMAPNO s(MAPING(J) »J=1sMAPINO) 
READ(591) ERR*AINC1sAINC2»(CORR(J) sv=192) 

999 ILEV=0 


Hira Os 

STRIP=06 

READ(5s1) HOLEsBITCSTsRIGCSTsTRIP Ts REP TI! sP 
99 READ(592) INOsINOS1»INOS2sIND1sIND2sIND3 


DEEVLSbREV+tL 
ILEV=ILEV+INO 
CPADGEVevEel SG) GO.Ta 98 
WRITE(696021) 

STOP 


98 READ(5s3) eee ae er ee sAF ORK 2 
Pe iac Weck ocd We (ABS(TEMP=$SDEPTH) eGTeOs.5) } 


LWRITE(G 96020) 
a PWOTSETVEPTH=RSDEPTA) /FLOAT ICING) 

BO TO veLeeVis TLEV 
DEPTH(JPSSDEPTHSFLOATCIFIX(FLOATCUF=LEEVI)*#FO0T+065) ) 
AF (GU) =AA 

CF (UY) =ACF 

DFS(U)=ADFS 

ALEMBEObs LIHRXLOt]) 

XLIMLO(J92)=XLO{ 2) 

XLIMHI(J91)=XHI(1) 

XLIMHIT(U¥s2)=XHI (2) 

FORM( J) =AFORM 
FORMK (J) =AFORM 
NOS1L ( ws WO Sl 
NOS2(J)=INOS2 

10 CONTINUE 

TEMP=TDEPTH 
IFCIND2sEQe0) GO 


IF IND2 1S NO RO THE SIZES OF THE STATE VARIABLE 
GRID ARE READ IN HERE 


AAD 


READ(592) (NOSL(4U) sJ=ILEVIsILEV) s(NOS2(4) sJFILEV1» 
a he 
20 TFUIND1L«EQs60) GO TO 99 
ILEVMI=STLEV 
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C MAIN eoee(CON 


ILEVSILEV+1 
NOSI(TILEV)=9 
NOS2(ILEV)=9 
DEPTH(ITLEV)=TDEPTH 
DO 38 J=1s9 
DO 37 K=159 
COTATI EV io K Y= e 
CONTINUE 
CONTINUE 
BITHR=BITCST/RIGCST 
TRIPHR=1000e#TRIPT 
KOUNT=1 
OO IZ21 KKFlsITLEVM1 
IF CKOUNTseNEsl) GO TO 31 
WRITE(6 96 
WRITE(696008) 
KOUNT=6 
ce k KOUNTeLTe61) GO 10743 
KOUNT=1 
(30 Oo 
KOUNT=KOUNT +2 


IP=IFIX(P/009)4+1 
XBAR=/e88/HOLE 

TFCFORMIKK )eGTeO05) ITUK=2 
IFUFORM(KKieLTe#O005) LUK#1 
WRITE(696009) DEPTHI(KK) »sDEPTH(K 
DFS(KK) sRESP(IUK) sFORM(KK) »XLI 
pALIMEQTK Kod) sXLIMHT (KK 92) 9 HO LE 
SBNOSL(KK) *sNOS2(KK) sLAG»wEARCIP) 
CONTINUE 

DOUT=0 « 


DOULZO K=LsILEVMi 
Se=TLEVMISLl=K 

KOUNT=KOUNT #1 
FOOTIN=DEPTH( JJ+1)=-DEPTH( Uy) 
eee 

ies Sree) ka TO Ls 

XNOS= FLOATINOSI (Jp <1. 

DO 15 J=lsKL 

DET QYY=FLOAT(OQ=17 7 ANOS 
DE(1)=Q6« 

IF(NOS2Z (Ud sEQel} GOTO as 
KMSNOS2 (Ue) 

XNOS= REGhT INOS aise ket 


DO 16 Je1sK™M 
BB J)=FLOAT(J-1)/XNOS 
BB(1)=0« 


fetvudreQel) GOTO 120 
ChuL LEVEL NS) 
TFLTOPTI«eNEsO) CALL RESULT (OW) 
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C MAIN eee (CONT'D) 


fF (TOPT2Ze«EQe0) GO To 126 
BO 150. MM=] sMAPNO 

150 [F(susEQ@eMAPING(MM)) CALL MAPStuu 

120 CONTINUE 
CALL TEST2(XLIMLO( 191) sXLIMLO(192) sTEFF slslsTANS#1) 
CALL PIND(XLIMLO( 191) sXLIMLO(1 92) slelels TERE) 
WIN(L)=Ww( 151) 
RIN(CL=RCLs1) 
KOUNT=0 
DO 130 JguH=leTLEVvMl 
FOOTIN=DEPTH(Jut])=DEPTH( UU) 
IF ( JUsEQel) GO TO 8 
Cac leo 2 (Ke IMO ed 
IFYLANSsEQs1) GO TO. Li 
CALL FINDCUXLIMLO( Juel) sXLIMLO(J092) sleds JusTEFF) 
REPLACSTIMEO+TRIPT*DEPTH(Uu)+REPTIM+BITCST/RIGCST 
PP REP ACL Tat OPAL (Geis Bde Go Te 1) 
CALE TEST 2(Welal) sR 0 Pel) sTERF ei vis IANS eat) 
TOTRTFTOTRT+DELT 
IDES=IFIX(DOUT%8s+0.5) 
IBB=IFIX(BOUT#S8.+0.5) 


PeXLIMLOC(Uds2) STEFF s leis lANSsdl) 


NSEC=NSEC+1 
WRITE (696006) NSECs CEPTH( JJ) sDEPTH( uutl) sWllalhs 
IR( lel) sDELTsIDEsI68 
- TE PE ORT OAT 0 Jule ach, wD) 
DE (1) =DOCUT 
BB(1)=BOUT 
Let sch Qe TLEVMIT ) -GO> TO. 21 


2 Ew 
11 TTRIP=TRIPT#DEPTH( UU) /2e+REPTIM/2e+TTRIP 
TEQ=TIRIP+TOTRI+BITCST/RIGCST 
CEQ=CEQ+TEO 
CTRIP=CTRIP+TTRIP 
CROT=TOTRT+CROT 
\ lel edad 1) FSTART» DEPTH( UU) sTOTRTsTTIRIPsTEQs 
11D FsI8B 
WRITE(696002) DEPTH(Jd) 
SEQ=SEQ+TEO 
STRIP=STRIP+TTRIP 
SROT=SROT+TOTRT 
CRIG=CROT+CTRIP 
SRIG=SROT+STRIP 7 . 
WRITE(6:6010) SROT»STRIP»SRIGsSEO sKOUNT - 
IF (IND3eEQe0) WRITE(696003) CROTsCTRIP»sCRIGsCEQsKKOUNT 
| TF (ddeEQeILEV) GO TO 130 
8 TTRIPHTRIPT#HDEPTH(JU)/2e+REPTIM/26 
FSTART= DEPTH( JJ) 
KOUNT=KOUNT#L1 
KKOUNT=KKOUNT+1 
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C MAIN eae (CONT'D) 


NSEC =) 
VRITE(6 36004) 
WRITE (696005) 


Wi lell=WIN( JU) 
R(lyLIFERIN( JU) 
DE(1)=06 
BB(1)=O05 
CALL TEST2(Wilsl)sR¢1l41)sTEFFslslsIANSsJv) 
TOTRT=DELT 
IDE=TFIX(DOUT#@e+005) 
PBBETFIX(E sat iret inde ) 
WRITE(696006) NSECs DEPTH (UU) sDEPTH( Uutl) sWllel)s 
R(lsl)sDELT»sIDE»IBB 
PIMOS) OT AL. THU eg les 7) 
DE(1)=DOUT 
RR(1)=R0UT 
Lei oeeGe lr CEYVMT ) GO ‘TO 2a 
CONTINUE 
READY v2) USETLsNUSET2 
DF Gh IUSETleNE«0) GO TO 999 
TF (NUSET2~NEs0) GO TO 9999 
STOP 
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CY RY CIsey 


re 


ANS) 


SUBROUTINE RESUL 
SUBROUTINE RESULT au) 
THIS PRINTS OUT A STAGE OF THE PROCESS SHOWING THE 


OPTIMAL PROCEDURE FOR THAT STAGE FOR EACH OF THE STATE 
VARTABLE GRID POINTS. 


COMMON iis fe cea DF Si1309 
COMMON FORM{150) sFOR! ES DerPiht 152) 


° 
OMMON XLI* NLO(15042) »XLIPHT(15002)» IOS D1 157) 
COMMON WIN(1L50) sRIN(L50) s TOTALT(151.15515) sNOS2(151) 
OMMON XNDIC (15515) 9W(15015)sR(15915) sDEC15) MeeBt led 
COMMON ERR AINCIsAINC2 sREPTIMsCORR( 2) sTRIPT »P 
COMMON HOLE sSITCSTsRIGCST sFOOTINs DOUTs» BOUT sLAGsDELT 
Ale heat DEPTH (uu) *FOOTIN 
PORMATC 1 '958Xs't INITIAL CONDITIONS's/s51Xs'DERTH = * 
LeFGsOv2Xoa INTERVAL) = '5F5.087) 
FORMAT (8 La A /3 
FORMAR(': TEETH) 
END=(NOS2 (UU) =1) /541 
KLENOSL (UU) 
KMeNOS2 (JU) 


TFAKMeEQei)° GO TO lil 
KOUN 


WRITE(6 9102 (BB ld) seelakM) 


BO. Sy we lak 

a Ni A 

TFC KOUNTseLEe60) GO TO 105 

WRITE(6 ees 

neck (BBi tL) aL. H=lskmM) 
RIiITE(Sea109) 

KOUNT= Lb 

WRITE(6s200) (WlU9K) sK=l9KM) 

WRITE(6Gs201) (RC Usk) 9K=Ll5K™) 

WRITE(69202) DE(U) »(TOTALT (yuus09K) sK=19KM) 

WRITE(6 3203 (XNDIC(UsK) s9K=1l9K™M) 

CONTINUE 

GO7i TOt7 

WRITE(G 9112) (BB(U) svsl eK) 

WRITE(69109) 

DOW 143) Vale Re 

KOUNT= KOUNT+9 

‘old TeLEs60) GO TO 115 
WRITE (69300) 

Beemre te 12) (BB(L) sk =l9KM) 

WRITE (69109) 

KOUNT=15 

WRITE(69210) (WlusK) 9K=l9KM) 
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SUBROUTINE Ea || o0oe CCONT ' D) 


WRITE(6s271) (RiUsK) sK=lskKM) 
| SE ihe 2) DEC) s (TOTALT (uu rusk VekKe=l eM) 
WRITE ( 5213) (XNOLC US cK kK =o KM) 
pid relies gid 
GO Ta: 7 
6 WRITE(6s122) (BARU) suelskmM) 

WRITE(69109) 

DO 123 J=bsKL 

K OUR Whe Ne hata 

IF (KO JUNTeLE«60) GO TO 125 

WRITE (69300). 


WRITE(6 3122) (SB(E) sl =hekM) 


ra) 
er 
iw 


(WC JekK) sKelsKM) 
(Ri Usk)» Peete! 


WR : DE CH) etTOTALT (alist) sK=H 1 kM) 
WRITE (69223 (XNDIC(UeK) sel 9 KM) 

veo CONTINUE 

T CONTINUE 
3 RETURN 

Zi IND=(NOS1{J0)=2) (oe OA 
GO-TO (l2vleeL4?) sINb 

PA OWRIFES 631227) (DEVO) sJelekKhL ) 
WRITE(Gs220) (WlUsl)sJU=1sKL) 
WRITE(G9221) (Rluel) suelskKL) 
WRITE(6 9227). (TOTALT(JUausl) a0=l sKL) 
WRITE LGe223) CAND ERClUs is J= lyk) 
RETURN 

Ce WeRiTetoell 7?) (DEL) sJ= Lak) 
WRITE(69210) (Wiel) sJ=lsKL) 
WRITE (652212) (Ridel)svelskKL) 
WRITE(69217) (TOTALT(Udeuel) steal sKL) 
WRITE(69213) (XNDIC(J91) sU=1sKL) 
RETURN 

14 WRITE(63107) (DE(U) sJ=l9KL) 
WRITE(6:200) (Wllsl) Je oe 
WRITE(62201) (RiUsl) sUelsKL 
WRITE(6s207) (TOTAL Ee eae 
WRITE(62203) (Xi ORE (ded awa! 
RETURN 


107 FORMAT(62Xs'DULLNESS! s/s1l2Xe15(3XsFS503)) 
117 FORMAT(62Xs!DULLNESS's/s11Xs10(7X oF D6 Ce. 
iar FORMAT (62Xs'DULLNESS!9/s SXsS(1LIXeF 560 ae 


ZO07 FORMAT ( hie Ness nae Lee ee 
Z2L7 FORMAT ( Js9Xs' COST! os 2XeFGe22 D(4XsFBe2)} 
227 FORMAT ( Js9Xs' COST! s BX sFBe2e4(lLbXsFBe2)) 


102 FORMAT(62X»9'BEARINGS! »/s12X915(3X9F503)) 
112 FORMAT(62X»! BEARINGS! »/slIXslOl(7XsF5e3)) 
122 FORMAT(62X%s'BEARINGS' s/s SxXsS5119X9F 503) ) 


4 ' % 
Geis 


Me un 


= 


7 - 
; Ch 
es | el a epaeed a 

a be ed T% he oe a. : 
es «lene, Lawl ECR ait 


= 


elshe (LGN) ee 


ms \ VOI 
Fa ab tw 7 

os f= - hehe Nausea VaPteld 
Leoke eat a 
__ Pitas 


cis {Ae J bis Lr 
{ MicteleeGh Si eet ey) 
: t\ 1 (00 92 R71 e8 
i= (> a1 H S30 oPRFre 
BYP TIA 
Rear ata 
rast. 
L=e Li-(uaple wie 
“wie { efLeBlt + 
{ lacy ti) i 4 SA ed 
(8¢ Legetieuw )) (0Reee 
aleageleri' i {ic aeg 
= ators fue y eS Td i Ss aia ‘ 
eitives Devi dt *) (¢ShekI8 autos 
: SUT 


, de Deve (1 90% (Xi cee ee EE 
-pwefeket bev) wi (Riis elaine 

(jie bee tLeut 9) “ERS of)9 Pw 

~. eie(deGeuGITuAToOLh is ‘ea ar gts 


plegel Pen yaic cry tae tive) 3te 


( )Oe Tele ty} 20 (POReeeaT, 
(jsélvietice?) (egal arar 
(jveimuatLeL li 
( we foie {le Gece Pha Top 
- ; prs Ana aa 


SUBROUTINE RESULT ses(CONT!'D) 


FORMAT(//97Xs'WEIGHT! 91XsF6elsl4(2X%sFGel) 
FORMAT (//s7Xs'"WEIGHT! s2XsFSe2s 9(4XsFBe2) 
FORMAT (//97Xs'WEIGHT!98XsF8e294(16XsF8e2) 
FORMAT (/s8X9'SPEED! s1XsF60el914(2XsF60l)) 
FORMAT (/s8Xs'!SPEED! 52X sFBe2s G(4X%eFBe2) ) 
FORMAT(/:s8Xs!SPEED! 98XsFBe2s4(16EXsFEe2) ) 
ORMAT(FGe35/99X st COST! as 1X sF60el914(2X%sF60el)) 
FORMAT (CF G039/59X 6 "COST! 2X sFBe2s 904K FB e 2) ) 
FORMAT (F6035/99X s'COST' s8XsFSe2s4(1EXsF8e2)) 
FORMAT(/ s7Xs!POLICY! s4xXsAls14(7X%sAL) ) 

FORMAT (7s 7Xs'POLICY!s6XsAls9(11XsAl)) 

FORMAT (/s7Xs'POLICY! 912X%9A194(23X%sA1) } 

END 
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SUBROUTINE LEVEL 


SUBROUTINE LEVEL (uy) 


THIS SUBROUTINE CALCULATES THE OPTIMAL POLICY FOR 
EACH OF THE STATE VARIABLE GRID POINTS FOR A GIVEN 
STAGE. 


DIMENSION RESP(2) 


COMMON AF (150) sCF(15¢) »DFS(150) 
COMMON FORM(150) sFORN« (150) sDEPTH( 151) 
COMMON XLIMLO(150%2)sXLIMHI (15092) sNOS1(151) 


COMMON WEN( 150) sRIN(150) sTOTALT (151915515) sNOS2(151) 
OMMON ANDIC( 15915) ow (Los le eRALS4e1S) spe tis) sehen 2) 
COMMON ERR*AINCLsAINC2ZsREPTIMsCORR(2)sTRIPT oP 
COMMON RS BE TEHE VES ELST 9 UR i SOM ee beer »DELT 
DATA RESP/IK'  sIRI/ 

BO de Le ks. 15 

Poe ie od 5 


XNDIC(I sL)=RESP(1) 

CONTINUE 

NUSIML=NOS TA Oy = 1 
FI(NOSLMIsEQse0) NOSIM1=1 

NOS2M1=NO0S2 (UU) = I 

PR(INOS2M1*eEQe0) NOS2Z2M1l= 

CA te S tne’ XLIMLO( Jul) XLIMLO(Uu92) 9TEFFol9lsTANS»uu) 
Cae PENDIXLIMLOU JI s1) sXLiIMLOl Qiak ll sielewurPerr. 
REPLACE GtALT (ele Lie TRIPT ADEPT AA OBL ICS iy RiGte. 

L+REPTIM 

DO 32 IelsNOSIMI 

DO 4 L=lsNOS2M1 

Te CCL EEGs) eANDa (LeEGel)} GO TO 4 

CALL TEST2Z(XLIMLO( Josl ls XLIMbLO( Use de CERF £2 Pom ranSee 
Pex DANS«EQal)) GO 7O\ 5 

CHILL FIND(XLIMLOS Jol) XL IMLOlUIs2) 9 Teles ERR) 
DEX TOPRALT( Woe lebleGTeREPLAC) GO TO 3 

CONTINUE 

GO TO 2 

tet LweeQal? GO ime) 6 

DO ¢* K=LsNOS2ZM1 

TOTALT(JUsTsK)=REPLAC 

WiIsK)=W(1lsl) 

Ri IsK)=R( lel) 

XNDICCIsK)=RESP(2) 

CONTINUE 

CONTINUE 

Goo TO) 10 

DO 8 J=I»NOSI1M1 


DO 9 K=L»sNOS2M1 
TOTAL T (Jee K DEREPLAC 
WiwesK)ewllsl) 
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SUBROUTINE LEVEL 


R(UsK)=E=R( 191) 
XNDIC(UsK)=RESP(2) 
CONTINUE 

CONTINUE 

KLENOS2 (uu) 

KM=NOSI] (uJ) 

TF tKMeEQs1) GO TO 13 
OGL ake Le Rie 
TOTALT(UJsKMsK)=REPLAC 
W(KMsK)SW(lel) 
R(KMsK)FR (191) 
XNDIC(KMsK)=#=RESP(2) 
CONTINUE 

LPvaLeeGed) “GO TO 14 
DO 12 K=lsNOSI1M1 
TOTAL T (Ud skKsKL)=REPLAC 
W(KeKL)=wllsl) 
R(KsKL)JFR( lel) 
XNDIC(KsKL)=RESP (2) 
CONTINUE 
WIN( UU) =WClel) 
RIN(JJ)=ER( 191) 

RETURN 


leat Dans 
END 


eee (CONT'D) 
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THATS *SUBROUTINE PRODUCES A MAP OF TH 
WHEN CALLED. 


SUBROUTINE MAPS 


SUBROUTINE MAPS( JU) 


DIMENSTON WT(16)sRV(26)sTIME( 
DI} IENSION ae ake: 26 } 
COMMON ot 50) sDFS(150 
COMMON FORM( 150 ae R 9D 
: (] 


COMMON ALIMLO(150 Vs XAL IMAL (015092) sNOS1T¢1 
COMMON WIN(150) »RIN(150) sTOTALT(151915815 
COMMON XNDIC(15515) sW(15915) »R(15515) DET 
COMMON ERR»sAINC] sAINC2sREPTIMsCORR(2)sTRI 
COMMON HOLE sBITCSTsRIGCST sFOOTIN»s DOUT s BOUT 


3, 


COMMON we sXABAR 

COMMON LOPT3sIOPT4s IOP T5s10PT6 

RANGE L= epee ithe LD) XE TMEO CI a ) 
RANGE Z=KLIMHI(JJs2)=xXLIMLO( vus2) 
IRI=ITFIA(RANGE2Z+0.5)+1 
IRL=MINO(IR1s26) 

wih 1 K=1s16 

WICK) =XLIMLO(J491)+FLOAT(K=1) *RANGELSI 
DO 2 K=lsIR1 


—_ 


= 
) »sNOS 
LSiae 
Py 
l 9 


9? 


LAGsDELT 


RV(K)=XLIMLO(JJs2) +FLOAT(K=1)*RANGE2/FLOAT (IRI=1) 


NOS i LeNOST( J =1 


IF C(NOSIMIeEQeO) NOSIM1=1 
PPXCNOST AI) sbGel ) LOPT3S=1 
NOS2MLSNOS2(Ju)=1 
IF(NOS2M1leEQe0) NOS2M1l=1 
PEt NOS2tJa)ieEGa dl). TOPT Ss =i 
DO 100 1=IOPT3sNOSIM1sTIOPT4 
DO 101 L=ITOPTSs NOS2M1siOPT6é 


DO 103 K=l1sIR1 

DO 102 J=1516 

CALL TEST2Z(WT(U) @RVIK) sTEIME(UsK) sislalaA 
DT (UsK)=D0OUT 

Bites h) =SOUT 

TFC TIME (Usk) eLTeleE+06) GO TO 102 
TIME( UsK)=00 

DT(JsK)=26 

BTiJsK)=2Z.s 

CONTINUE 

IND=0 

DO 99 IJK=1916 

ee ea ce, IND=IND+1 
TEL INDseNEal6) GO TO 103 


IR1=«K 
GO TO 300 
ic ON a if NUE 
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200 


204 
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SUBKOUTINE MAPS wee (CONT'D) 


WRITE(69200) DEPTH( uu) sFOOTINSDE(L) sBB{L) 
PORMAT(C'1L'sS8Xs"INITIAL CONDITIONS's//s32X%s'DEPTH = ! 


LeFGe092Xs!' INTERVAL = "esF5e052Xs' DULLNESS =° 'o F604 62% 
2¢'BEARINGS = "oF6e45/) 


WRITE (63204) (WT(U) su=1916) 

FORMAT (64X9"WEIGHT! s//s20Ks16(1X9F6e2) 9/93X%9'SPEED') 
KOUNT#8 

DO 104 K=lsIR1 

KOUNT=KOUNT +4 

TFISOUNT*LE*60) GO. .70 105 


KOUNT=4 
oe 
FORMAT(!'1'* 


Pearce scons (WIS) sy=1916) 


Weide V6a20)) (TIME Lusk i rueleie) 
WRITE (69202) RV(K) 9 ( (DT CUsK) svelel6) 
WRITE(G6 9203) (BT lUsK) sdJ=1516} 
CONTINUE 

CONTINUE 

CONTINUE 
PORMATCLSK 9 COST 's3%5 160 1XsF 602) } 


FORMAT (2XsF6e293Xs'DULLNESS!' s1xXs16(2X5F5e3)) 
FORMAT(1LIXs' BEARINGS! »s1Xs16(2X5F 563) 5/) 
RETURN 

END 
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THIS SUBPROGRAM FINDS THE TOTAL MINIMUM COS 
IT EXAMINES BOTH THE INTERNAL MINIMUM AND THE ED 
CONDITION AND TAKES THE LOWER VALUE. 


H= 14 


SUBROUTINE FIND 


SUBROUTINE FIND{WATEINSREVINe TsLouue TIN! 


LOGICAL AsBsCaDsEvF 


DIMENSION X12) 9X2 (2) axL2 (2) 

COMMON AF { Oi sCF(150) sBFS(150) 

COMMON FORM(150) sFORMK(150) 9b SEP TH (151) 

COMMON XLIMLO(1503s2)sXLIW APtLe0 2 ysNOSL ELL) 

COMMON WIN(1L50)sRIN(150) s TOTALT (151915915) oN 24 2on4 
GOMMON XNDIC(15915)eW(15Ss15 JeR( 15915) 9DE Fioy eerie 
COMMON riddled \CLsAINC2ZsREPTIMsCORR(2)sTRIPT sp 
COMMON HOLEsBITCST sRISCSTsFOOTIN »DOUT o£ ACUT sLAGsDELT 
SS=l.s 

LL=90 

KL=0 

K2(1L)=WATEIN 

K2'2)=REVI! 

DO 1 J=els2 

pa Ge ee an ie 

FBESTSTIN 

Smid te ear dk St RRR Re ae alld dl J} 


TF (ABS (XL1(1) feb TeleE=10eANDeABS(XLI(2)}eLTelaE=10) 


DO 4 Jela2 

RE CQ ERL CY) =SSEXET (OV CORR CO) 

PE CRS tI eGTeXL IMA] (Jos) ) X2t as BRKEIMET (oa) 
TRUK e Ra hel Te ALIMLO Ve) AS) SREP MLO (ed? 
CALL TEST 2X!) eX2( 2) sFXstslLsTANSeue) 
IF(FXeLTeFREST) GO TO 6 

ITF ¢SSsLT*« ERR) GO TO 100 

SS=SS/2e 

Lobe ed 

GO TO 5 

FBEST=F x 

KL=KL41 

TF(LLeGTeO) GO TO 30 

SS=SS#2. 

Poe pedal ae 

RECO PEK it J )=SS% ‘XLII (U)* ORR (UU) 

TEYKS (CU hPGTeXLIMHINdds do) ) KZCOTERLENAT (Cs) 
Teun e to eal exe Eh bade age Keto ext IMEC 
CALL TEST2(X%261)5X2(2) sFXslsleTANSsuu) 
PDFCFX«eLTeFREST) GO TO 6 

SS=SS/26 

pO § JRlez 

Bot gyex ht byt =SS#KCL( #CORR( UW) 


avg yi 
SuGr Meuuad 
22 — 


oe 
> 


aha 
IugRets)s aie 
vis Seo Ri 


(gS eet) Le 
ae 
* eo) Tive tS erie 1\ 
i 7 i _ ~~ 4 aw = pa res ee 
oud oT 98 © 
Seleu @. 
wi ~s(. jtixne t2=(UT Peet sx 
a ( e*Leil Ljtwie »ferGvest 
5 4 she (usu ae Tin BTS 
{ Ri hsee re (ulema tae ie (I bd A eS 


3) (Tee 24 T sak 1? 
Lit. Oyen rhs r 


i focus 
. ipecaatasy 
ay i. 7 min.’ 


G& 
Vw 


ens 


Dg 


on 
<> 
© 
te 


1003 


1000 


SUBROUTINE FIND 
TP UXL (od eG Te XGIMMI (ove) ) 
LEI SLIW i OE TeXLIMLOt dae) 
AZABS(XL(1LI@XLIMHI (Quel) Jeb TeOel 
SEABS (XE (2) =XLIMHI(JJ92) ol TeOel 
IF(AsANDeB) GO TO 103 

CABS(XL(2)@XLIMLO(JJs2) ) eb Teds u, 
IF CAsAND eC) GO 70 160 


DSABS UXT CLI =XLIMLO( Jel) Jel Ts Oe l 


IF(DeANDsB) GO TO 100 
SS=2e%*5S/2 o*#*KL 

LL=0 

KL= 

GO TO 9 

CONTINUE 

CACK TEST2(XLIMA IT Cited ja XLIMAD C2 } 
DEC IMAXeL Ter BEST) GO TO. 102 
TFC TMAXeGTelseE+06) GO TO 102 
WtTslL)=X1(1) 

RilsbLiexl (2) 
TOTALT(JusTsL)=FR8EST 


RETURN 

WOLsLISXLIMHI(ld91) 
R(IsLI=XLIMHI(Uus2)} 
TOTALT(J0eTsbL)=TMAX 


RETURN 


RANGE L=(XLIMHT(J0el)=XLIMLO( Jos] ) 
LOtIwJe2 ] 


RANGE2=(XLI! 
XGRID=XLIMLO 
YGRID=RLIMLO(JUs2 
XTEMPSX1(1 

YTEMP =) 


PT el eae LG 
(JUeal)+RANGEL 
)+PANGE2 
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ITF(RANGEZeL1e00e13) GO 
a eee Bo GO 10 
TTEMP=TT 

KTEVP= XGRID 

See agli ite 
XGRID=XGRID+RANG 
YGRID= ceriee NG 
GO TO 1001 
XGRID=XGRID-RANGE 
¥GRID=YGRID=<“RANGE 
So 10 LOO] 


Kd bod = KL. PI 
AL(I)SXLIMLO ( 
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SUBROUTINE TEST2 


SUBROUTINE TEST2(WEICHT »REVS »TEFFoTslLsTANSsJJ) 


THIS SUBROUTINE CALCULATES THE TOTAL COST TO DATE IF 
A GIVEN WEIGHT IS USED AT THE CURRENT STAGE. 


COMMON AF (150) sCF (150) sDFS(150) 

COMMON FORM(150) sFORMK (150) sDEPTH(151) 

COMMON XLTIMLO(15092)sXLIMHI(15092) »NOS1(151) 5NOS2(151) 
COMMON WIN(150) sRIN(150) sTOTALT(151915915) 

COMMON XNDIC(15915) 9W(15515)sR(15s15) sDE(15) »BB(15) 
COMMON ERR AINC1 sAINC2 sREPTIMsCORR(2) sTRIPT sP 

COMMON HOLE sBITCSTsRIGCSTsFOOTINsDOUTsBOUTSLAG »sDELT 
COMMON IP»sXBAR 


DULLA(D)= 1e04(600+.6928125%D) *D 
DULLAT(D)=(61004(364+.309375%*D) #D) *D 
DULLAH(D)=(1e0+(600+.9281258D) *D) ¥¥0.5 
SIM1l(D)=(1.0197794(1.292304+( .1133279%0—-.4127675) 
¥D) #D ) %*D 
IANS=0 
Az ne apeee te Weeki 
LPC FORM( OU) =045) Lsls2 

eee A*REVS#H#O 6s (540 e5*#REVS*(Le=A) 


Mow 3 

ARR=A#REVS#¥0 0428400 28REVS#(1le—A) 
WP ie =) RAL 3, EIGHT 
EMWATE=1359e1l=—714e19#ALOG1O(WBAR) 
ROTI=REVS#( le + CO004348 #REVSERO) 
KLOAD=B1338e¢87/(LleteOZB1LIB7*WBAR ) ¥¥#3 02258 
GO TO (20920920) sIP 
Fl= ices aac er ROTI 
hem uae * (WBRARH*#FORMK (UU) )#ARR 
F3=FI*RE a ea Reais 
GO TO (40%50930) 5IP 
X=DE( 1) +062 
G=SIML(DE(T) }+FOOTINZAF 2 
OOUTHXK—=(SIMLIX) Cc) /DULLAH(X) 

[FLABS{( (DOUT=K) /DOUT) «LTeleE=05)} GO 10 5 


TF{(DOUTseGTels) GO TO 10 
BOITFFEDULLAT(DOUT)=—DULLAT(DE(T)) 
BOUT=F3*DDIFF+BB(L) 
TFP{BOUTeGTels) GO TO 10 

CALL INPOL(TPREVelsLeu) 
DELT=FL*DDIFF 

TEFFH=TPREV*DELT 

eo 70 11 

TEFF=1seE+10 

IANS=1 


es 
SUBROUTINE TEST2 eese(CONT'D) 


11 RETURN 

20 DOUT=FOOTIN/F2+DE(1) 
rEVOOUPSeGTsl6e0) GO. To" 15 
DDIFF=S=DULLAI (DOUT)=DULLAI(DE(1)) 
BOUT=F32*DDIFF+B8R(L) 
ITF(BOUTs«GTele0) GO 0 -Eo 
CALL INPOL(TPREVsTsLou¥u) 
DELT=FI*DDIFF 
TEFF=TPREV+DELT 
RETURN 

30 FI=DFS(J4)*XLOAD/REVS 
BOUT=FOOTIN/ (F1¥*¥CF (UU) * (WBAR*#*FORMK (JU) ) *#ARR)+BBI(L) 
Pe VooUT SS ha leur GOO po 
GALL INPORTPPREVesleiswal 
DELT=F1l* (BOUT=BB(L)) 
TEFFHRTPREV+DELT 
RETURN 
END 
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SUBROUTINE DERIV 
SUBROUTINE DERIV(WEIGHT»REVSsDTEFFolslLsTEFFlsuy) 


THIS SUBROUTINE CALCULATES THE VORMALIZED DERIVATIVES 
OF COST WeReTe WEIGHT AND SPEED. 


DIMENSION DTEFF(2) 

COMMON AF(150)9CF(150) sDFS(150) 

COMMON FORM(150) sFORMK(150) sDEPTH(151) 

COMMON XLIMLO615052) sXLIMHI (1502) sNOS1(151) 

COMMON WIN(150)sRIN(150) s TOTALT(151915815)sNOS2(151) 
COMMON XNDIC(15915) oW(15915)sR(15915)sDE(15) »68B8(15) 
COMMON ERRSAINCISAINC2 sREPTIMsCORR(2)sTRIPT oP 
COMMON HOLE »*BITCST»RIGCST»FOOTIN»>DOUTsBOUTsLAG sDELT 
X=WEIGHT™=AINC1 

CALL TEST2(XsREVSsTEFF2s1 sLa ITANSs JU) 

X=REVS=AINC2 


GALL TESTE ZIWELGHTSASTERFSS I sie TANSs JU) 

ET SeR OLE, LEFF LeTERF2) /ATNGE 
DIBFRE(ZISATEFFI@]=TEFF3S)/AING? 

[IFt(ABSUDOTEFPF (1) }elLTelwE=20sANDsABS(DTEFFI( 2) aLlTs 


lisE=20) GO TO 100 
DEV=(DTEFF (1) *#2+DTEFF (2) **2) + 
OrERE WE =eDIErR i) / DEV 
DTERF(2Z)=DTEFF (2) /DEV 

100 RETURN 
END 
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SUBROUTINE INPOL 


SUBROUTINE INPOL(TPREVsIsksJu) 


TRIS SUBROUTINE FINDS THE COST UP To THE PREVIOUS 
STAGE BY INTERPOLATION. 


COMMON AF(150)sCF(150) sDFS(150) 
COMMON FORM(150) »FORMK(15 O) sDEPTH(151) 
COMMON ALIMLO(150s2) sXLIMHI (15052) sNOS1(151) 

OMMON WIN(150)eRIN (150) sTOTALT( Vis lel? plo) eNOS 2a) 
COMMON XNDIC(15515) 9W(15915)9R(15015) sDE(15)9BB ee 
COMMON ERRsAINCLsAINC2 sREPTIM Is CORR | Z2)9TRIFPT 9P 
eats HOLE sBITCST sRIGCSTsFOOTINsDOUTsBOUTsLAG sDELT 

NOSL=FLOAT(NOS1(Uu+1)=1) 
eee: FLOAT(NCS2(JJ+1)=1) 
J=LAGTI 


Ke J 
LevPrGieivo)® GO TO 2 
TPS DOUT sLEe (FLOAT y= 


J= J+] 
Leis bOa NOS Lt gel) “GO Ta) 2 
So tO! 2 


PEXNOSZ\ IIL) eEGs lt) SO TO"3 
PPV BOUT estes (FLOAT (Ki) 7ANOS2) 7 -GO 70 


i) 


LPikebOaNOs2 (Jor) C2 106 


J2=5 * 
LE VPS GT ale 5) JUZaul 


Kak | 
PE U(NOS2( CUtl}eEOsl) K2=K1 
TPREV= Ue 


DO 10 J#K1sK2 

TEMP=0 6 

OO 21. ivsJLewe 

XLl=le 

TF (PsGlTeleS) GO TO 15 

OO Le hel id til 2 

IF(IuseEQsKy) GO TO 12 

XLL=XL1#( DOUT*XNOS1=<FLOAT (Kum) ) /FLUAT (Tu-ky) 
CONTINUE 

TEMP=TEMP+XLI*TOTALT(Vutlsluss) 

CONTINUE 

XL2=l-e 

DO 13 KL=Kl»K2 

LE fe EG eKL) GO TO. 13 
XL2=XL2*(BOUT#XNOS2—FLOAT (KL=1))/FLOAT (u™ KL) 


CONTINUE 
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GROUP 


14006 
LoZ ts 
16538 
1780+. 
1907« 
20336 
21606 
22876¢ 
24136 
25406 
2667 3 
21936 
29206 
30266 
37 338 
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34516 
3o 56's 
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BIT NOs SECTION 


12 
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BIT 


NOs 


FwWNr FWhw PE 


SECTION 


1 


SECTION 


EWN foOMe FwWNhHY FHM re FW PDP 


INTERVAL 
59208 = 60266 
6026.6 si 61336 
61336 ae 62396 
62396 e' G34 5.6 
63455 = 6451. 
645le = 6556. 
65588 = 6664. 
66646 = 6770s 


GROUP NOe 2 


INTERVAL 


677Cse = 6810.6 
GROUS NOs 3 
INTERVAL 
6BL008° = 6S Le 
Otte = COT bs 
CODe ee Wee 
Live = “(2120 
Veen = peiae 
T3136 = °7S13% 
THES = OTS Pas 
7514s = 7614. 
PEA Shel SA et ee 
(ae) eee ee 
oLo% = T1916. 
T9lés. = BOSS 
SULoe: = Bllve 
Siliet* Shed ies 
S2ive —~ —8313 
6318s bes: Sa18 4 
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GROUP NOe 4 
INTERVAL 


88204 = 8866 
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INTERVAL 
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INTERVAL 
902068 = 9060. 


COA 
GROUP N 


Ose f 


INTERVAL 


OnLA 
5060 6 = 


2 aa 
Sligo a = 
SA ee = 
9400 = 


Voy 
77 tee = 
9878s = 
99796 =, 


L052. = 
LU628 6, = 


Gus Gx eae 


LO¢eL = = 
108136 ee 
108446 ™ 
LOE f aa = 


109066 = 
LOO sae = 
LO969« = 
ELOUG a 
Loe. = 
110636 


110946 =) 


tee a 


91136 
SLo. § 
9218. 
92706 
9371.6 
94736 
9574.6 
96756 
9716s 
98786 
he Me 
LOUSUs 


LOi926 
LO303e 
104156 
LOQO63G6 
LO75U06 
LO7SB 1 
L1U844 
LOS i Se 
LU90Ge 
109386 
LlOoSs6éS« 
LLOOOs 
VEO be 
110636 
L1LU94e 
liane 
LPL Ts 


WELGHT 


606 


WEIGHT 


Dae 
60 
60e 
606 
oe 

Ove 
506 
568 
606 
K j 

Si @ 
606 
Pay es 

Cue 


if. 
Ww 


SPEEL 


1 


bod 


UUe 


EED 


WM hy eR Wb 
e 


opie z 
wliagel 
uu 
«95 \Gs 


baie 


t252 = Ope 2 Ey : 


te 


Ww 


uw 


CUMULATIVE 
CUMULATIVE 
CUMULATIVE 
CUMULATIVE 
CUMULATIVE 


& LL1S vs ee 
1 L1188. = 
Z LL 219s 
2 
4h 


Ve Sls = 


Ll282s 


LPLoots = 
11738. 
bat Soar = 
B= eB is 
LiQlle rs 
SR ine ies 
r202oses 
120846 a 
es Oia Pe 
Leecon = 
ings eee 


L22ese = 
125403 = 
LZ6168 = 
T2692 a 
12768. rz 

2844s. = 
Lee eo = 
Lieto 


WNMr ££ Wh FP FF Wh FE 


~ 
> 


Bi ee 
USN OLE Weep tan 
pars 
1 OS ik ee 

ave. vs 


Oo. -!1 0 
° 
i 


PFUNFEWNE FWONHE PUD E 
= 


Ww WhN be © 


3 
3 

3 Se = 
aha ae 
2352 = 
s45)e.= 


a ee ee ae ee 


ROTATING TIME 
TRIP TIME 

RIG TIME 

COST 

BIT COUNT 


INTERVAL 


112826 
pp Ne vee 5 


NO» ¥ 


INTERVAL 


INTERVAL 


11738. 
LPiso% 
LigS3< 
ak dala 
L1969.% 
120266 
120846 


rR pt PS Re 
o 

CE 

o 2 


Dh) BN) Ro RG PO DO PO PG 


=! o Os Ui & GW Wo Po 
or DW 
ine) 
es 


2 
a 


ee Ss 
> oO 
Or 
a a 3 2.40 5S al 2 a 
e . 2s e R 


i 


wc 


PO NO PR PO 
s © 
¢ 


en ae ee eo 
OW ~1h FO 


! fd 
i) Ww 
me Oo © 
s 


wi G ~) G ® 
e 2 


2a 


~ he 


- 
5 Ww Ww 


Wo Po 


WETGHT 


~W Po Ww 
o 


r 
~~ 
e 


wm f$wore Ww 
e 


ui 
cS 


Srrie 


106 


SPEED 


Oronoovuwdodcro sd 
e = © © e a e ° oe ° e Ss @ 


OO NN itn © 
a @ 


Fwy How SHOW HAO ww AS Sw 


* 


>O re 
e 


RIG HOURS 
RIG HOURS 
RIG HOURS 
RIG HOUKS 


: 
* 
2 
. 
° 
° 
> ? 
o 
* 
. . 
.* 
PP» 
sue 
aCi 
' < 
. :”* 
e ~~ *«.Fr 
Zt 
~ en 
: = 
; + 
4 , «a1 


" a 


cmeitictes 


; ehitte 
¥ i & LA Hs 


jaw A io 2 ee. a 


r 
* « eo > 


” 
Twi iT) 
— . rps 
— ] ' 


H 
‘ 
2 F 
— 
- 
ne Bre 


i 
* . ~ * vows 
on = 8 

ve & 


! 
a 
’ 
v 
a + 
- 


A 
«4° 
.¢ 
oat 
i) q 
&<* 
Fr 
eS « 
, 


th, Me & Mie Ow th 


— 
Tel L1T3a2. «0A one | 


tf 


7 "Ta 
ie! 


